00001
00002 include '../Include_file/include_modulefiles_peos.f90'
00003 include '../Module/grid_parameter_binary_excision.f90'
00004 include '../Module/grid_points_binary_excision.f90'
00005 include '../Module/grid_points_asymptotic_patch.f90'
00006 include '../Module/weight_midpoint_binary_excision.f90'
00007 include '../Module/radial_green_fn_helmholtz.f90'
00008 include '../Module/radial_green_fn_hrethadv.f90'
00009 include '../Module/radial_green_fn_hrethadv_homosol.f90'
00010 include '../Module/copy_array_3d.f90'
00011 include '../Module/copy_array_4d.f90'
00012 include '../Module/make_char2_array_2d.f90'
00013 include '../Module/make_char1_array_2d.f90'
00014 include '../Module/make_int_array_3d.f90'
00015
00016 include '../Include_file/include_interface_modulefiles_peos.f90'
00017 include '../Module_interface/interface_error_metric.f90'
00018 include '../Module_interface/interface_error_metric_type0.f90'
00019 include '../Module_interface/interface_grdr_gridpoint_type0.f90'
00020 include '../Module_interface/interface_grdr_gridpoint_type0_nosym.f90'
00021 include '../Module_interface/interface_sourceterm_helmholtz_solver_test.f90'
00022 include '../Module_interface/interface_sourceterm_exsurf_eqm_binary.f90'
00023 include '../Module_interface/interface_sourceterm_outsurf_eqm_binary.f90'
00024 include '../Module_interface/interface_sourceterm_surface_int.f90'
00025 include '../Module_interface/interface_sourceterm_surface_int_homosol.f90'
00026 include '../Module_interface/interface_sourceterm_insurf_asymptotic_patch.f90'
00027 include '../Module_interface/interface_helmholtz_solver_vol_int.f90'
00028 include '../Module_interface/interface_helmholtz_solver_surf_int.f90'
00029 include '../Module_interface/interface_helmholtz_solver_binary.f90'
00030 include '../Module_interface/interface_helmholtz_solver_binary_vol_int.f90'
00031 include '../Module_interface/interface_helmholtz_solver_binary_surf_int.f90'
00032 include '../Module_interface/interface_helmholtz_solver_outer_surf_int.f90'
00033 include '../Module_interface/interface_helmholtz_solver_asymptotic_patch_homosol.f90'
00034 include '../Module_interface/interface_copy_to_bsjy_and_sbsjy.f90'
00035 include '../Module_interface/interface_interpo_binary_to_asymptotic_patch.f90'
00036
00037 include '../Include_file/include_modulefiles_mpatch.f90'
00038
00039 include '../Include_file/include_subroutines_peos.f90'
00040 include '../Subroutine/IO_output_poisson_test_3D.f90'
00041
00042 include '../Subroutine/error_metric.f90'
00043 include '../Subroutine/error_metric_type0.f90'
00044 include '../Subroutine/printout_error_metric.f90'
00045 include '../Subroutine/sourceterm_helmholtz_solver_test.f90'
00046 include '../Subroutine/sourceterm_exsurf_eqm_binary.f90'
00047 include '../Subroutine/sourceterm_outsurf_eqm_binary.f90'
00048 include '../Subroutine/sourceterm_surface_int.f90'
00049 include '../Subroutine/sourceterm_surface_int_homosol.f90'
00050 include '../Subroutine/helmholtz_solver_vol_int.f90'
00051 include '../Subroutine/helmholtz_solver_surf_int.f90'
00052 include '../Subroutine/helmholtz_solver_binary.f90'
00053 include '../Subroutine/helmholtz_solver_binary_vol_int.f90'
00054 include '../Subroutine/helmholtz_solver_binary_surf_int.f90'
00055 include '../Subroutine/helmholtz_solver_outer_surf_int.f90'
00056 include '../Subroutine/helmholtz_solver_asymptotic_patch_homosol.f90'
00057 include '../Subroutine/grdr_gridpoint_type0.f90'
00058 include '../Subroutine/grdr_gridpoint_type0_nosym.f90'
00059 include '../Subroutine/test_source_helical_binary.f90'
00060 include '../Subroutine/allocate_poisson_solver_test.f90'
00061 include '../Subroutine/allocate_metric_and_matter_BHNS_test.f90'
00062 include '../Subroutine/test_analytic_solution.f90'
00063 include '../Subroutine/calc_radial_green_fn_hrethadv.f90'
00064 include '../Subroutine/calc_radial_green_fn_hrethadv_homosol.f90'
00065 include '../Subroutine/interpo_binary_to_asymptotic_patch.f90'
00066 include '../Subroutine/copy_to_bsjy_and_sbsjy.f90'
00067 include '../Subroutine/sphbess.f90'
00068 include '../Subroutine/sphbess_and_dx.f90'
00069 include '../Subroutine/sphbess_dx.f90'
00070 include '../Subroutine/bessjy.f90'
00071 include '../Subroutine/beschb.f90'
00072 include '../Include_file/include_subroutines_mpatch.f90'
00073 include '../Subroutine_mpatch/iteration_helmholtz_solver_binary_test_mpt.f90'
00074 include '../Subroutine_mpatch/IO_output_poisson_test_3D_mpt.f90'
00075 include '../Subroutine_mpatch/allocate_metric_and_matter_BHNS_test_mpt.f90'
00076 include '../Subroutine_mpatch/copy_metric_and_matter_BHNS_test_to_mpt.f90'
00077 include '../Subroutine_mpatch/copy_metric_and_matter_BHNS_test_from_mpt.f90'
00078 include '../Subroutine_mpatch/copy_to_mpatch_helmholtz_test.f90'
00079 include '../Subroutine_mpatch/copy_from_mpatch_helmholtz_test.f90'
00080 include '../Subroutine_mpatch/test_source_helical_binary_mpt.f90'
00081 include '../Subroutine_mpatch/sourceterm_insurf_asymptotic_patch.f90'
00082 include '../Function/chebev.f90'
00083
00084
00085
00086
00087 PROGRAM Main_helmholtz_binary_test_mpt
00088
00089 use phys_constant, only : nmpt
00090 use grid_parameter, only : indata_type, outdata_type, &
00091 & iter_max
00092 use grid_parameter_binary_excision
00093 use grid_points_binary_excision
00094 use grid_points_binary_excision_mpt
00095 use grid_points_asymptotic_patch
00096 use weight_midpoint_grav
00097 use weight_midpoint_fluid
00098 use weight_midpoint_binary_excision
00099 use def_vector_x
00100 use def_vector_phi
00101 use radial_green_fn_helmholtz
00102 use radial_green_fn_hrethadv
00103 use radial_green_fn_hrethadv_homosol
00104 implicit none
00105 integer :: iseq, iter_count, total_iteration, impt, impt_bin
00106
00107
00108 call allocate_grid_parameter_mpt
00109 call allocate_grid_parameter_binary_excision_mpt
00110 call read_parameter_mpt(1)
00111 call copy_grid_parameter_to_mpt(1)
00112 call read_parameter_binary_excision_mpt(1)
00113 call copy_grid_parameter_binary_excision_to_mpt(1)
00114
00115 call read_parameter_mpt(2)
00116 call copy_grid_parameter_to_mpt(2)
00117 call read_parameter_binary_excision_mpt(2)
00118 call copy_grid_parameter_binary_excision_to_mpt(2)
00119
00120 call read_parameter_mpt(3)
00121 call copy_grid_parameter_to_mpt(3)
00122 call read_parameter_binary_excision_mpt(3)
00123 call copy_grid_parameter_binary_excision_to_mpt(3)
00124
00125
00126 call set_allocate_size_mpt
00127
00128 call allocate_coordinate_patch_kit_grav_mpt
00129 call allocate_grid_points_binary_excision
00130 call allocate_weight_midpoint_grav
00131 call allocate_weight_midpoint_fluid
00132 call allocate_weight_midpoint_binary_excision
00133 call allocate_grid_points_asymptotic_patch
00134 call allocate_vector_x
00135 call allocate_vector_phi
00136 call allocate_radial_green_fn_helmholtz
00137 call allocate_radial_green_fn_hrethadv
00138 call allocate_radial_green_fn_hrethadv_homosol
00139 call allocate_metric_and_matter_BHNS_test
00140 call allocate_metric_and_matter_BHNS_test_mpt
00141
00142 call allocate_mpatch_all_test
00143
00144 do impt = 1, nmpt
00145 call copy_grid_parameter_from_mpt(impt)
00146 call copy_grid_parameter_binary_excision_from_mpt(impt)
00147 call coordinate_patch_kit_grav_mpt
00148 if (impt.eq.1) then
00149 call calc_parameter_binary_excision
00150 call calc_grid_points_binary_excision
00151 call calc_weight_midpoint_binary_excision
00152 call test_source_helical_binary_mpt(impt)
00153 end if
00154
00155 if (impt.eq.2) then
00156 call calc_parameter_binary_excision
00157 call calc_grid_points_binary_excision
00158 call calc_weight_midpoint_binary_excision
00159 call test_source_helical_binary_mpt(impt)
00160 end if
00161
00162 if (impt.eq.3) then
00163 impt_bin = 1
00164 call calc_parameter_binary_excision
00165 call calc_grid_points_asymptotic_patch(impt_bin,impt)
00166 end if
00167 call calc_vector_x_grav(0)
00168 call calc_vector_x_matter(0)
00169 call calc_vector_phi_grav(0)
00170 call calc_vector_phi_matter(0)
00171 call copy_to_mpatch_helmholtz_test(impt)
00172 call copy_metric_and_matter_BHNS_test_to_mpt(impt)
00173 end do
00174
00175 call iteration_helmholtz_solver_binary_test_mpt(iter_count)
00176 if (total_iteration.ge.iter_max) then
00177 write(6,*)' ** Solution did not converge **'
00178 end if
00179
00180 do impt = 1, nmpt
00181 call copy_from_mpatch_helmholtz_test(impt)
00182 call copy_metric_and_matter_BHNS_test_from_mpt(impt)
00183 if (outdata_type.eq.'3D') call IO_output_poisson_test_3D_mpt(impt)
00184 end do
00185
00186
00187
00188
00189
00190 END PROGRAM Main_helmholtz_binary_test_mpt
