crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 14:15:42 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Added ability to treat multiple interpolation regions in interpolate_tab1 functions.

Browse source
This commit is contained in:
Paul Romano 2011-11-12 10:37:57 -08:00
parent e3e35dec49
commit 1d73c6d509
2 changed files with 22 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

1
src/DEPENDENCIES
View file

@ -100,6 +100,7 @@ interpolation.o: endf_header.o
interpolation.o: error.o
interpolation.o: global.o
interpolation.o: search.o
interpolation.o: string.o
logging.o: constants.o
logging.o: global.o

25
src/interpolation.f90
View file

@ -5,6 +5,7 @@ module interpolation
use error, only: fatal_error
use global, only: message
use search, only: binary_search
use string, only: int_to_str
implicit none
@ -29,6 +30,7 @@ contains
real(8) :: y ! y(x)
integer :: i ! bin in which to interpolate
integer :: j ! index for interpolation region
integer :: loc_0 ! starting location
integer :: n_regions ! number of interpolation regions
integer :: n_points ! number of tabulated values
@ -81,8 +83,12 @@ contains
elseif (n_regions == 1) then
interp = data(loc_interp + 1)
elseif (n_regions > 1) then
message = "Multiple interpolation regions not yet supported."
call fatal_error()
do j = 1, n_regions
if (i < data(loc_breakpoints + j)) then
interp = data(loc_interp + j)
exit
end if
end do
end if
! handle special case of histogram interpolation
@ -111,6 +117,9 @@ contains
case (LOG_LOG)
r = (log(x) - log(x0))/(log(x1) - log(x0))
y = exp((1-r)*log(y0) + r*log(y1))
case default
message = "Unsupported interpolation scheme: " // int_to_str(interp)
call fatal_error()
end select
end function interpolate_tab1_array
@ -129,6 +138,7 @@ contains
real(8) :: y ! y(x)
integer :: i ! bin in which to interpolate
integer :: j ! index for interpolation region
integer :: n_regions ! number of interpolation regions
integer :: n_pairs ! number of tabulated values
integer :: interp ! ENDF interpolation scheme
@ -159,8 +169,12 @@ contains
elseif (n_regions == 1) then
interp = obj % int(1)
elseif (n_regions > 1) then
message = "Multiple interpolation regions not yet supported."
call fatal_error()
do j = 1, n_regions
if (i < obj % nbt(j)) then
interp = obj % int(j)
exit
end if
end do
end if
! handle special case of histogram interpolation
@ -189,6 +203,9 @@ contains
case (LOG_LOG)
r = (log(x) - log(x0))/(log(x1) - log(x0))
y = exp((1-r)*log(y0) + r*log(y1))
case default
message = "Unsupported interpolation scheme: " // int_to_str(interp)
call fatal_error()
end select
end function interpolate_tab1_object