diff options
Diffstat (limited to 'src/petsc_wrapper.c')
| -rw-r--r-- | src/petsc_wrapper.c | 14 |
1 files changed, 8 insertions, 6 deletions
diff --git a/src/petsc_wrapper.c b/src/petsc_wrapper.c index 1f25b3e..74d3512 100644 --- a/src/petsc_wrapper.c +++ b/src/petsc_wrapper.c @@ -7,17 +7,17 @@ static const char *rcsid = "$Header$"; CCTK_FILEVERSION(CactusElliptic_EllPETSc_petsc_wrapper_c) -/* The wrapper functions for the core PETSc solver, that performs the actual solve. - One wrapper function for each class, because these functions are registered - with a fixed set up arguments. Still, they all call petsc_confmetric_solver - as the core solver. */ +/* The wrapper functions for the core PETSc solver, that performs the + actual solve. One wrapper function for each class, because these + functions are registered with a fixed set up arguments. Still, they + all call petsc_confmetric_solver as the core solver. */ -/* Wrapper function for the class of ellitpic equations that needs a metric */ +/* Wrapper function for the class of elliptic equations that needs a metric */ void petsc_metric(cGH *GH, int *MetricI, int FieldIndex, int MIndex, int NIndex, CCTK_REAL *AbsTol, CCTK_REAL *RelTol) { - /* Arrays for the M/N size info */ +/* Arrays for the M/N size info */ void petsc_confmetric_solver(cGH *GH, int *MetricI, int MetricISize, int FieldIndex, int MIndex, int NIndex, @@ -25,7 +25,9 @@ void petsc_metric(cGH *GH, int *MetricI, int FieldIndex, if (GH->cctk_dim>3) + { CCTK_WARN(0,"This elliptic solver implementation does not do dimension>3!"); + } petsc_confmetric_solver(GH, MetricI, 6, FieldIndex, MIndex, NIndex, |