tstLibrary.f90

./Core/dc/f/tests/tstLibrary.f90

program tstlibrary

#include "ScaleSTL/FortranTestMacros.h"

use origen_library_m
use origen_fakefactory_m
use origen_definitions_m

implicit none

type(origen_library)              :: lib
type(origen_libraryheader)        :: definition
type(origen_decaydata)            :: decay_data
type(origen_transitionstructure)  :: tstruct
type(origen_transitioncoeff)      :: tcoeff
type(origen_nuclideset)           :: nucset

type(origen_fakefactory)          :: ff
type(scalestl_vec_flt)            :: fvec
type(scalestl_vec_int)            :: ivec
real(C_FLOAT),allocatable         :: farr(:)
integer(C_INT),allocatable        :: iarr(:)

integer,parameter :: l=c_size_t
logical :: loaded

call lib % initialize()
loaded = ff % Library_random1(lib,333)
if( .not. loaded )stop 1

!
! Get different parts of the library.
!
call lib % getptr_definition(definition)
assert_c_assoc(definition)

!test the header
call nucset % initialize()
call definition % get_nuclide_set(nucset)
expect_eq(742,nucset % num_in_sublib(origen_sublib_1lt))
expect_eq(556,nucset % num_in_sublib(origen_sublib_2ac))
expect_eq(928,nucset % num_in_sublib(origen_sublib_3fp))

!test decay data
call lib % getptr_decay_data(decay_data)
assert_c_assoc(decay_data)
expect_near( 1.50642300e-3 , decay_data%decay_constants(1) , 1e-5)

!test transition structure
call lib % getptr_transition_structure(tstruct)
assert_c_assoc(tstruct)
call tstruct % getptr_num_parents_vec(ivec)
assert_c_assoc(ivec)
expect_eq(6,ivec%at(1_l))
expect_eq(2226_l,ivec%size())
!try with native arrays
call tstruct % get_num_parents(iarr)
expect_eq(6,iarr(1))
expect_eq(2226,size(iarr,1))

!test transition coeff
call lib % getptr_transition_coeff_at(tcoeff,1_l)
assert_c_assoc(tcoeff)
call tcoeff % getptr_matrix_vec(fvec)
assert_c_assoc(fvec)
expect_near(8.69777426e-03,fvec%at(1_l),1e-5)
!try with native arrays
call tcoeff % get_matrix(farr)
expect_near(8.69777426e-03,farr(1),1e-5)

call lib % destroy()

end program