diff options
author | miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2000-03-10 14:05:04 +0000 |
---|---|---|
committer | miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2000-03-10 14:05:04 +0000 |
commit | 8247894d2a3055a83200eeea3a1639f91f026c4a (patch) | |
tree | 20e774b06dee25ea83e4eca6513b6a19acc5597f | |
parent | c9be9cbae8cc14223059293e26eeea8fdc84e095 (diff) |
Adding hook to BAM.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@22 a678b1cf-93e1-4b43-a69d-d43939e66649
-rw-r--r-- | src/brilldata.F | 33 |
1 files changed, 21 insertions, 12 deletions
diff --git a/src/brilldata.F b/src/brilldata.F index 1b1645d..4d5862a 100644 --- a/src/brilldata.F +++ b/src/brilldata.F @@ -29,20 +29,20 @@ c Set up background metric and coefficients for linear solve. call setupbrilldata3D(CCTK_FARGUMENTS) end if -c Call the Linear Elliptic solver interface -c to find conformal factor. +c Call the Linear Elliptic solver interface to find +c conformal factor. First some preliminaries. - call CCTK_VarIndex (metpsi_index(1), "einstein::gxx") - call CCTK_VarIndex (metpsi_index(2), "einstein::gxy") - call CCTK_VarIndex (metpsi_index(3), "einstein::gxz") - call CCTK_VarIndex (metpsi_index(4), "einstein::gyy") - call CCTK_VarIndex (metpsi_index(5), "einstein::gyz") - call CCTK_VarIndex (metpsi_index(6), "einstein::gzz") - call CCTK_VarIndex (metpsi_index(7), "einstein::psi") + call CCTK_VarIndex(metpsi_index(1),"einstein::gxx") + call CCTK_VarIndex(metpsi_index(2),"einstein::gxy") + call CCTK_VarIndex(metpsi_index(3),"einstein::gxz") + call CCTK_VarIndex(metpsi_index(4),"einstein::gyy") + call CCTK_VarIndex(metpsi_index(5),"einstein::gyz") + call CCTK_VarIndex(metpsi_index(6),"einstein::gzz") + call CCTK_VarIndex(metpsi_index(7),"einstein::psi") - call CCTK_VarIndex (field_index, "IDBrillData::brillpsi") - call CCTK_VarIndex (Mlin_index, "IDBrillData::Mlinear") - call CCTK_VarIndex (Nsrc_index, "IDBrillData::Nsource") + call CCTK_VarIndex(field_index,"IDBrillData::brillpsi") + call CCTK_VarIndex(Mlin_index,"IDBrillData::Mlinear") + call CCTK_VarIndex(Nsrc_index,"IDBrillData::Nsource") AbsTol(1)= brill_thresh AbsTol(2)= -1 @@ -52,6 +52,8 @@ c to find conformal factor. RelTol(2)= -1 RelTol(3)= -1 +c Boundaries. + if (CCTK_EQUALS(brill_bound,"const")) then call Ell_SetRealKey(ierr,brill_const_v0, . "EllLinConfMetric::Bnd::Const::V0") @@ -64,6 +66,8 @@ c to find conformal factor. . "EllLinConfMetric::Bnd::Robin::inf") endif +c Elliptic solver. + if (CCTK_EQUALS(brill_solver,"sor")) then call Ell_LinConfMetricSolver(ierr,cctkGH,metpsi_index, . field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"sor") @@ -74,6 +78,11 @@ c to find conformal factor. . field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"petsc") end if + if (CCTK_EQUALS(brill_solver,"bam")) then + call Ell_LinConfMetricSolver(ierr,cctkGH,metpsi_index, + . field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"bam") + end if + c Synchronize conformal factor. call CCTK_SyncGroup(cctkGH,"einstein::confac") |