ETS  \$Id: Doxyfile 2162 2020-02-26 14:16:09Z g2dpc $
 All Classes Files Functions Variables Pages
check_equilibrium.f90
Go to the documentation of this file.
1 program check_equilibrium
2 !----------------------------------------------------------------------------
3 ! This small program compares the newly calculated equilibrium
4 ! file 'equilibrium.cpo' with the reference equilibrium file
5 ! 'equilibrium_opteron' in the current directory.
6 ! It uses the functions in the module error_analysis to find
7 ! differences between the two equilibria.
8 ! It helps you verify the quality of your equilibrium calculation and trace
9 ! bugs with the existing test cases.
10 !----------------------------------------------------------------------------
11 
12  use itm_types
13  use deallocate_structures
14  use read_structures
15 ! use diff_structures
16  use error_analysis
17  use euitm_schemas
18 
19  implicit none
20 
21  type (type_equilibrium), pointer :: reference(:), test(:)
22  character(len = 132) :: filename_reference_equilibrium
23  character(len = 132) :: filename_new_equilibrium
24  character(len = 132) :: path
25  integer(itm_i4), parameter :: iu6 = 6
26 
27 !-- define default behaviour without input arguments
28  path = './' ! run locally
29  filename_reference_equilibrium = 'equilibrium_opteron'
30  filename_new_equilibrium = 'equilibrium.cpo'
31 
32 !-- allocate empty CPOs
33  allocate(reference(1))
34  allocate(test(1))
35 
36 !-- read equilibrium CPO from file (reference equilibrium)
37  call set_read_verbosity(0) ! no verbose output
38  call open_read_file(12, trim(adjustl(path)) &
39  // trim(adjustl(filename_reference_equilibrium)))
40  call read_cpo(reference(1), 'equilibrium')
41  call close_read_file
42 
43 !-- read equilibrium CPO from file (test equilibrium)
44  call set_read_verbosity(0) ! no verbose output
45  call open_read_file(12, trim(adjustl(path)) &
46  // trim(adjustl(filename_new_equilibrium)))
47  call read_cpo(test(1), 'equilibrium')
48  call close_read_file
49 
50 !-- compare both equilibrium CPOs
51  call set_diff_verbosity(0) ! only summary output
52  call reset_diff_counter
53  call set_error_level(0._r8)
54  call diff_cpo(reference(1), test(1), 'equilibrium', are_identical_float)
55 ! write(iu6, *) 'total number of differences: ', get_diff_counter()
56  call diff_cpo(reference(1), test(1), 'equilibrium', average_relative_error_float)
57 ! write(iu6, *) 'final average relative error: ', get_error_level()
58  call diff_cpo(reference(1), test(1), 'equilibrium', maximum_relative_error_float)
59 ! write(iu6, *) 'final maximum relative error: ', get_error_level()
60  call diff_cpo(reference(1), test(1), 'equilibrium', chi_square_float)
61 ! write(iu6, *) 'final chi^2: ', get_error_level()
62 
63 !-- final deallocate
64  call deallocate_cpo(reference)
65  call deallocate_cpo(test)
66 
67 end program check_equilibrium
error_analysis::maximum_relative_error_float
character(len=132) function maximum_relative_error_float(diff_counter, error_level, reference_field_float, new_field_float, reference_field_array3dflt_type, new_field_array3dflt_type, reference_field_array4dflt_type, new_field_array4dflt_type, reference_field_array5dflt_type, new_field_array5dflt_type, reference_field_array6dflt_type, new_field_array6dflt_type, reference_field_array7dflt_type, new_field_array7dflt_type, reference_field_matflt_type, new_field_matflt_type, reference_field_vecflt_type, new_field_vecflt_type)
Definition: error_analysis.f90:738
error_analysis::are_identical_float
character(len=132) function are_identical_float(diff_counter, error_level, reference_field_float, new_field_float, reference_field_array3dflt_type, new_field_array3dflt_type, reference_field_array4dflt_type, new_field_array4dflt_type, reference_field_array5dflt_type, new_field_array5dflt_type, reference_field_array6dflt_type, new_field_array6dflt_type, reference_field_array7dflt_type, new_field_array7dflt_type, reference_field_matflt_type, new_field_matflt_type, reference_field_vecflt_type, new_field_vecflt_type)
Definition: error_analysis.f90:75
error_analysis::chi_square_float
character(len=132) function chi_square_float(diff_counter, error_level, reference_field_float, new_field_float, reference_field_array3dflt_type, new_field_array3dflt_type, reference_field_array4dflt_type, new_field_array4dflt_type, reference_field_array5dflt_type, new_field_array5dflt_type, reference_field_array6dflt_type, new_field_array6dflt_type, reference_field_array7dflt_type, new_field_array7dflt_type, reference_field_matflt_type, new_field_matflt_type, reference_field_vecflt_type, new_field_vecflt_type)
Definition: error_analysis.f90:1177
error_analysis::average_relative_error_float
character(len=132) function average_relative_error_float(diff_counter, error_level, reference_field_float, new_field_float, reference_field_array3dflt_type, new_field_array3dflt_type, reference_field_array4dflt_type, new_field_array4dflt_type, reference_field_array5dflt_type, new_field_array5dflt_type, reference_field_array6dflt_type, new_field_array6dflt_type, reference_field_array7dflt_type, new_field_array7dflt_type, reference_field_matflt_type, new_field_matflt_type, reference_field_vecflt_type, new_field_vecflt_type)
Definition: error_analysis.f90:309
error_analysis
This module is an example for the error_analysis module which shall hold the user defined error analy...
Definition: error_analysis.f90:19
check_equilibrium
program check_equilibrium
Definition: check_equilibrium.f90:1