diff options
author | miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2000-03-10 13:32:12 +0000 |
---|---|---|
committer | miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2000-03-10 13:32:12 +0000 |
commit | c9be9cbae8cc14223059293e26eeea8fdc84e095 (patch) | |
tree | 48869d3ccd4d9810831678c97c9ad40007d3cd42 /src | |
parent | 3a21b0e6628dbed7fadbaaff37ae815dbaf2e3c8 (diff) |
Extending elliptic interface to work with petsc, and adding some
parameters for boundary conditions.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@21 a678b1cf-93e1-4b43-a69d-d43939e66649
Diffstat (limited to 'src')
-rw-r--r-- | src/brilldata.F | 25 |
1 files changed, 22 insertions, 3 deletions
diff --git a/src/brilldata.F b/src/brilldata.F index 65614f7..1b1645d 100644 --- a/src/brilldata.F +++ b/src/brilldata.F @@ -44,7 +44,7 @@ c to find conformal factor. call CCTK_VarIndex (Mlin_index, "IDBrillData::Mlinear") call CCTK_VarIndex (Nsrc_index, "IDBrillData::Nsource") - AbsTol(1)= 0.0001 + AbsTol(1)= brill_thresh AbsTol(2)= -1 AbsTol(3)= -1 @@ -52,8 +52,27 @@ c to find conformal factor. RelTol(2)= -1 RelTol(3)= -1 - call Ell_LinConfMetricSolver(ierr,cctkGH,metpsi_index, - .field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"sor") + if (CCTK_EQUALS(brill_bound,"const")) then + call Ell_SetRealKey(ierr,brill_const_v0, + . "EllLinConfMetric::Bnd::Const::V0") + end if + + if (CCTK_EQUALS(brill_bound,"robin")) then + call Ell_SetIntKey(ierr,brill_robin_falloff, + . "EllLinConfMetric::Bnd::Robin::falloff") + call Ell_SetRealKey(ierr,brill_robin_inf, + . "EllLinConfMetric::Bnd::Robin::inf") + endif + + if (CCTK_EQUALS(brill_solver,"sor")) then + call Ell_LinConfMetricSolver(ierr,cctkGH,metpsi_index, + . field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"sor") + end if + + if (CCTK_EQUALS(brill_solver,"petsc)) then + call Ell_LinConfMetricSolver(ierr,cctkGH,metpsi_index, + . field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"petsc") + end if c Synchronize conformal factor. |