new KSPSystemSize current value

Interface/ProData.cpp

@@ -220,6 +220,7 @@ void Init_ProblemStructure()
   Current.KSPIterations = 0.0;
   Current.KSPIteration = 0.0;
   Current.KSPResidual = 0.0;
+  Current.KSPSystemSize = 0.0;
 }
 
 // FIXME: TODO to remove parser memory leaks!

Interface/ProData.h

@@ -1647,7 +1647,7 @@ struct CurrentData {
   double  RelaxFac, NbrTestedFac, SolveJacAdaptFailed; //+++
 
   // Iterative linear system solvers
-  double  KSPIterations, KSPIteration, KSPResidual;
+  double  KSPIterations, KSPIteration, KSPResidual, KSPSystemSize;
 };
 
 /* ------------------------------------------------------------------------ */

Interface/ProDefines.h

@@ -332,6 +332,7 @@ struct StringXPointer  Current_Value[] = {
   {"KSPIterations", &Current.KSPIterations}, {"KSPIts", &Current.KSPIterations},
   {"KSPIteration", &Current.KSPIteration},
   {"KSPResidual", &Current.KSPResidual},
+  {"KSPSystemSize", &Current.KSPSystemSize},
   {NULL       , NULL}
 } ;
 

Kernel/LinAlg_PETSC.cpp

@@ -1182,6 +1182,8 @@ static void _zitsol(gMatrix *A, gVector *B, gVector *X)
   PetscInt n;
   _try(VecGetLocalSize(B->V, &n));
 
+  Current.KSPSystemSize = n;
+
   int *row = (int*)Malloc(nnz * sizeof(int));
   int *col = (int*)Malloc(nnz * sizeof(int));
   double *valr = (double*)Malloc(nnz * sizeof(double));
@@ -1291,6 +1293,7 @@ static void _solve(gMatrix *A, gVector *B, gSolver *Solver, gVector *X,
 
   PetscInt i, j;
   _try(MatGetSize(A->M, &i, &j));
+  Current.KSPSystemSize = i;
   if(!i){
     Message::Warning("Zero-size system: skipping solve!");
     return;

doc/texinfo/getdp.texi

@@ -1569,6 +1569,8 @@ Current iteration in a Krylov subspace solver.
 Current residual in a Krylov subspace solver.
 @item $KSPIterations
 Total number of iterations of Krylov subspace solver.
+@item $KSPSystemSize
+System size of Krylov subspace solver.
 @end table
 
 @noindent Note: