WIP

docs/src/PETScPC/oseen.py.rst

@@ -139,8 +139,9 @@ We try a multiplicative preconditioner instead ::
          w[:] = 0
          dofs = BitArray(self.fes.ndof); dofs[:] = True
          smooth = KrylovSolver(self.a, dofs, p=smoother, 
-                              solverParameters={"ksp_type": "gmres", 
-                                                "ksp_max_it": 10, 
+                              solverParameters={"ksp_type": "fgmres", 
+                                                "ksp_max_it": 10,
+                                                "ksp_view": '',
                                                 "pc_type": "mat"})
          w += GMRes(self.a.mat, d, pre=smoother, x=w, maxsteps = 10, printrates=False)
          r = d.CreateVector()

ngsPETSc/ksp.py

@@ -6,7 +6,7 @@
 from ngsolve import la, BilinearForm, FESpace, BitArray, Projector
 from ngsPETSc import Matrix, VectorMapping, PETScPreconditioner
 
-def createFromBilinearForm(a, freeDofs, solverParameters, optionsPrefix):
+def createFromBilinearForm(a, freeDofs, solverParameters):
     """
     This function creates a PETSc matrix from an NGSolve bilinear form
     """
@@ -23,7 +23,7 @@ def createFromBilinearForm(a, freeDofs, solverParameters, optionsPrefix):
         mat = Matrix(a.mat, (dofs, freeDofs, None), solverParameters["mat_type"])
     return (a.mat, mat.mat)
 
-def createFromMatrix(a, freeDofs, solverParameters, optionsPrefix):
+def createFromMatrix(a, freeDofs, solverParameters):
     """
     This function creates a PETSc matrix from an NGSolve bilinear form
     """
@@ -39,7 +39,7 @@ def createFromMatrix(a, freeDofs, solverParameters, optionsPrefix):
         mat = Matrix(a, (dofs, freeDofs, None), solverParameters["mat_type"])
     return (a, mat.mat)
 
-def createFromPC(a, freeDofs, solverParameters, optionsPrefix):
+def createFromPC(a, freeDofs, solverParameters):
     class Wrap():
         def __init__(self, a, freeDofs):
             self.mapping = VectorMapping((a.dofs,freeDofs,{"bsize": 1}))
@@ -78,7 +78,7 @@ class KrylovSolver():
     :arg optionsPrefix: special solver options prefix for this specific Krylov solver
 
     """
-    def __init__(self, a, dofsDescr=None, p=None, solverParameters=None, optionsPrefix=None, nullspace=None):
+    def __init__(self, a, dofsDescr=None, p=None, solverParameters=None, optionsPrefix="", nullspace=None):
         # Grabbing parallel information
         if isinstance(dofsDescr, FESpace):
             freeDofs = dofsDescr.FreeDofs()
@@ -93,11 +93,11 @@ def __init__(self, a, dofsDescr=None, p=None, solverParameters=None, optionsPref
         #Construct operator        
         for key in parse:
             if isinstance(a, key):
-                ngsA, pscA = parse[key](a, freeDofs, solverParameters, optionsPrefix)
+                ngsA, pscA = parse[key](a, freeDofs, solverParameters)
         if p is not None:
             for key in parse:
                 if isinstance(p, key):
-                    ngsP, pscP = parse[key](p, freeDofs, solverParameters, optionsPrefix)
+                    ngsP, pscP = parse[key](p, freeDofs, solverParameters)
         else:
             ngsP = ngsA; pscP = pscA
         #Construct vector mapping
@@ -124,6 +124,7 @@ def __init__(self, a, dofsDescr=None, p=None, solverParameters=None, optionsPref
         self.ksp.setOptionsPrefix(optionsPrefix)
         self.ksp.setFromOptions()
         self.pscX, self.pscB = pscA.createVecs()
+        self.ksp.setUp()
 
     def solve(self, b, x):
         """

ngsPETSc/nullspace.py

@@ -2,11 +2,8 @@
 This module contains all the function and class needed to wrap a PETSc Nullspace in NGSolve
 '''
 from petsc4py import PETSc
-
-
 from ngsPETSc import VectorMapping
 
-
 class NullSpace:
     '''
     This class creates a PETSc Null space from NGSolve vectors