add a few more PC things

Author: JaredCrean2

src/PETSc.jl

@@ -72,6 +72,7 @@ function PetscInitialized()
 end
 
 function PetscFinalize()
+  println("finalizing PETSc")
   gc() # call garbage collection to force all PETSc objects be destroy that are queued up for destruction
   return ccall( (:PetscFinalize,  libpetsclocation),Int32,());
 end

src/petsc_error.jl

@@ -36,7 +36,7 @@ end
 # tell Petsc abut the Error handler
 cfunc = cfunction(error_handler, PetscErrorCode, (comm_type, Int32, Ptr{UInt8}, Ptr{UInt8}, PetscErrorCode, PetscErrorType, Ptr{UInt8}, Ptr{Void}) )
 ctx = C_NULL
-ierr = ccall( (:PetscPushErrorHandler, petsc), PetscErrorCode, (Ptr{Void}, Ptr{Void}), cfunc, ctx)
+#ierr = ccall( (:PetscPushErrorHandler, petsc), PetscErrorCode, (Ptr{Void}, Ptr{Void}), cfunc, ctx)
 
 
 function chkerrq(i::PetscErrorCode)

src/petsc_ksp.jl

@@ -1,5 +1,5 @@
 
-export KSP, KSPSetOperators, KSPSetFromOptions, KSPSolve,  KSPSetUp
+export KSP, KSPSetOperators, KSPSetFromOptions, KSPSolve, KSPSolveTranspose, KSPSetUp, KSP_NULL
 
 export KSPGetConvergedReason, PetscView, KSPSetType, KSPGetType, KSPSetTolerances, KSPGetTolerances, KSPSetInitialGuessNonzero, KSPGetInitialGuessNonzero, KSPGetResidualNorm
 

src/petsc_mat.jl

@@ -192,6 +192,7 @@ end
   end
 
   function MatSetOption(mat::PetscMat,arg2::MatOption,arg3::Union{Integer, Bool})
+      println("mat.pobj = ", mat.pobj)
       ierr = ccall((:MatSetOption, libpetsclocation),PetscErrorCode, (Ptr{Void},MatOption,PetscBool), mat.pobj, arg2, arg3)
 
       if ierr != 0

src/petsc_pc.jl

@@ -1,10 +1,11 @@
-export PC, PCCreate, KSPGetPC, KSPSetPC, PCSetType, PCGetType, PCFactorSetUseInPlace, PCFactorGetUseInPlace, PCSetReusePreconditioner, PCGetReusePreconditioner, PCFactorSetAllowDiagonalFill, PCFactorGetAllowDiagonalFill, PCFactorSetLevels, PCFactorGetLevels, PCSetReusePreconditioner, PCGetReusePreconditioner, PCBJacobiGetSubKSP, PCFactorSetFill, PCJacobiSetType, PCJacobiGetType
+export PC, PCCreate, PCSetFromOptions, KSPGetPC, KSPSetPC, PCSetType, PCGetType, PCFactorSetUseInPlace, PCFactorGetUseInPlace, PCSetReusePreconditioner, PCGetReusePreconditioner, PCFactorSetAllowDiagonalFill, PCFactorGetAllowDiagonalFill, PCFactorSetLevels, PCFactorGetLevels, PCSetReusePreconditioner, PCGetReusePreconditioner, PCBJacobiGetSubKSP, PCFactorSetFill, PCJacobiSetType, PCJacobiGetType
 
 # Shell PC
-export PCShellSetApply, PCShellSetSetUp, PCShellSetContext, PCShellGetContext
+export PCShellSetApply, PCShellSetApplyTranspose, PCShellSetSetUp,
+       PCShellSetContext, PCShellGetContext
 
 # developer PC interface
-export PCApply, PCApplyTranspose, PCApplyTransposeExists
+export PCApply, PCSetUp, PCApplyTranspose, PCApplyTransposeExists
 
 # preconditioner contex
 # the KSP object creates the PC contex, so we don't provide a constructor
@@ -16,11 +17,23 @@ global const PC_NULL = PC(C_NULL)
 
 function PCCreate(comm::MPI.Comm)
   arg2 = Ref{Ptr{Void}}() 
-  ccall((:PCCreate,petsc),PetscErrorCode,(comjm_type,Ptr{PC}),comm,arg2)
+  ccall((:PCCreate,petsc),PetscErrorCode,(comm_type,Ptr{Ptr{Void}}),comm,arg2)
 
   return PC(arg2[])
 end
 
+function PetscDestroy(arg1::PC)
+  if arg1.pobj != C_NULL
+    ccall((:PCDestroy,petsc),PetscErrorCode,(Ptr{Ptr{Void}},),&arg1.pobj)
+  end
+end
+
+
+function PCSetFromOptions(arg1::PC)
+    ccall((:PCSetFromOptions,petsc),PetscErrorCode,(Ptr{Void},),arg1.pobj)
+end
+
+
 function KSPGetPC(ksp::KSP)
     arr = Array(Ptr{Void}, 1)
     ccall((:KSPGetPC,petsc),PetscErrorCode,(Ptr{Void},Ptr{Void}),ksp.pobj, arr)
@@ -47,6 +60,11 @@ function PCShellSetApply(arg1::PC,arg2::Ptr{Void})
     ccall((:PCShellSetApply,petsc),PetscErrorCode,(Ptr{Void},Ptr{Void}),arg1.pobj,arg2)
 end
 
+function PCShellSetApplyTranspose(arg1::PC,arg2::Ptr{Void})
+    ccall((:PCShellSetApplyTranspose,petsc),PetscErrorCode,(Ptr{Void},Ptr{Void}),arg1.pobj,arg2)
+end
+
+
 function PCShellSetSetUp(arg1::PC,arg2::Ptr{Void})
     ccall((:PCShellSetSetUp,petsc),PetscErrorCode,(Ptr{Void},Ptr{Void}),arg1.pobj,arg2)
 end
@@ -151,6 +169,7 @@ function PCJacobiGetType(pc::PC)
 end
 
 function PCApply(arg1::PC,arg2::PetscVec,arg3::PetscVec)
+  println("pc.pobj = ", arg1.pobj)
     ccall((:PCApply,petsc),PetscErrorCode,(Ptr{Void},Ptr{Void},Ptr{Void}),arg1.pobj ,arg2.pobj, arg3.pobj)
 end
 
@@ -165,5 +184,8 @@ function PCApplyTransposeExists(arg1::PC)
     return arg2[] == 1
 end
 
+function PCSetUp(arg1::PC)
+    ccall((:PCSetUp,petsc),PetscErrorCode,(Ptr{Void},),arg1.pobj)
+end