From 8247894d2a3055a83200eeea3a1639f91f026c4a Mon Sep 17 00:00:00 2001
From: miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649>
Date: Fri, 10 Mar 2000 14:05:04 +0000
Subject: Adding hook to BAM.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@22 a678b1cf-93e1-4b43-a69d-d43939e66649
---
 src/brilldata.F | 33 +++++++++++++++++++++------------
 1 file 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")
-- 
cgit v1.2.3