37       static PetscErrorCode ierr;  variable
202         ierr = MatCreateBAIJ(PETSC_COMM_WORLD, petsc_bs, petsc_m, petsc_m, petsc_M, petsc_M,  in SolGMRp()
205         ierr = MatSetOption(lhsP, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);  in SolGMRp()
209         ierr = MatSetOption(lhsP, MAT_SUBSET_OFF_PROC_ENTRIES, PETSC_TRUE);  in SolGMRp()
211         ierr = MatSetUp(lhsP);  in SolGMRp()
214       ierr = MatGetOwnershipRange(lhsP, &myMatStart, &myMatEnd);  in SolGMRp()
218       if(genpar.lhs ==1) ierr = MatZeroEntries(lhsP);  in SolGMRp()
287         ierr = ISCreateGeneral(PETSC_COMM_SELF, petsc_n, indexsetary,  in SolGMRp()
293       ierr = MatAssemblyBegin(lhsP, MAT_FINAL_ASSEMBLY);  in SolGMRp()
294       ierr = MatAssemblyEnd(lhsP, MAT_FINAL_ASSEMBLY);  in SolGMRp()
302         ierr = MatGetLocalSize(lhsP, &LocalRow, &LocalCol);  in SolGMRp()
303         ierr = VecCreateMPI(PETSC_COMM_WORLD, LocalRow, petsc_M, &resP);  in SolGMRp()
304         ierr = VecCreateMPI(PETSC_COMM_WORLD, LocalRow, petsc_M, &DyP);  in SolGMRp()
306       ierr = VecZeroEntries(resP);  in SolGMRp()
308         ierr = VecCreateSeq(PETSC_COMM_SELF, petsc_n, &DyPLocal);  in SolGMRp()
328           ierr =  VecSetValue(resP, PetscRow, valtoinsert, ADD_VALUES);  in SolGMRp()
331       ierr = VecAssemblyBegin(resP);  in SolGMRp()
332       ierr = VecAssemblyEnd(resP);  in SolGMRp()
333       ierr = VecNorm(resP,NORM_2,&resNrm);  in SolGMRp()
342         ierr = KSPCreate(PETSC_COMM_WORLD, &ksp);  in SolGMRp()
343         ierr = KSPSetOperators(ksp, lhsP, lhsP);  in SolGMRp()
344         ierr = KSPGetPC(ksp, &pc);  in SolGMRp()
345         ierr = PCSetType(pc, PCPBJACOBI);  in SolGMRp()
348         ierr = KSPSetTolerances(ksp, timdat.etol, PETSC_DEFAULT, PETSC_DEFAULT, maxits);  in SolGMRp()
349         ierr = KSPSetFromOptions(ksp);   in SolGMRp()
352         ierr = KSPSetResidualHistory(ksp,resHist,maxitsHist, PETSC_TRUE);   in SolGMRp()
354       ierr = KSPSolve(ksp, resP, DyP);  in SolGMRp()
355       ierr = KSPGetResidualHistory(ksp,&resHist,&maxitsHist);   in SolGMRp()
363       ierr = VecScatterCreate(DyP, LocalIndexSet, DyPLocal, NULL, &scatter7);  in SolGMRp()
365       ierr = VecScatterBegin(scatter7, DyP, DyPLocal, INSERT_VALUES, SCATTER_FORWARD);  in SolGMRp()
366       ierr = VecScatterEnd(scatter7, DyP, DyPLocal, INSERT_VALUES, SCATTER_FORWARD);  in SolGMRp()
367       ierr = VecGetArray(DyPLocal, &parray);  in SolGMRp()
376       ierr = VecRestoreArray(DyPLocal, &parray);  in SolGMRp()
384       ierr = KSPGetIterationNumber(ksp, &its);  in SolGMRp()
529         ierr = MatCreateAIJ(PETSC_COMM_WORLD, petsc_m, petsc_m, petsc_M, petsc_M,  in SolGMRpSclr()
532         ierr = MatSetOption(lhsPs, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);  in SolGMRpSclr()
536         ierr = MatSetOption(lhsPs, MAT_SUBSET_OFF_PROC_ENTRIES, PETSC_TRUE);  in SolGMRpSclr()
538         ierr = MatSetUp(lhsPs);  in SolGMRpSclr()
541       ierr = MatGetOwnershipRange(lhsPs, &myMatStart, &myMatEnd);  in SolGMRpSclr()
546       if(genpar.lhs == 1) ierr = MatZeroEntries(lhsPs);  in SolGMRpSclr()
602         ierr = ISCreateGeneral(PETSC_COMM_SELF, petsc_n, indexsetarys,  in SolGMRpSclr()
607       ierr = MatAssemblyBegin(lhsPs, MAT_FINAL_ASSEMBLY);  in SolGMRpSclr()
608       ierr = MatAssemblyEnd(lhsPs, MAT_FINAL_ASSEMBLY);  in SolGMRpSclr()
611         ierr = MatGetLocalSize(lhsPs, &LocalRow, &LocalCol);  in SolGMRpSclr()
612         ierr = VecCreateMPI(PETSC_COMM_WORLD, LocalRow, petsc_M, &resPs);  in SolGMRpSclr()
613         ierr = VecCreateMPI(PETSC_COMM_WORLD, LocalRow, petsc_M, &DyPs);  in SolGMRpSclr()
615       ierr = VecZeroEntries(resPs);  in SolGMRpSclr()
617         ierr = VecCreateSeq(PETSC_COMM_SELF, petsc_n, &DyPLocals);  in SolGMRpSclr()
635           ierr =  VecSetValue(resPs, PetscRow, valtoinsert, ADD_VALUES);  in SolGMRpSclr()
637       ierr = VecAssemblyBegin(resPs);  in SolGMRpSclr()
638       ierr = VecAssemblyEnd(resPs);  in SolGMRpSclr()
639       ierr = VecNorm(resPs,NORM_2,&resNrm);  in SolGMRpSclr()
642         ierr = KSPCreate(PETSC_COMM_WORLD, &ksps);  in SolGMRpSclr()
643         ierr = KSPSetOperators(ksps, lhsPs, lhsPs);  in SolGMRpSclr()
644         ierr = KSPGetPC(ksps, &pcs);  in SolGMRpSclr()
645         ierr = PCSetType(pcs, PCPBJACOBI);  in SolGMRpSclr()
648         ierr = KSPSetTolerances(ksps, timdat.etol, PETSC_DEFAULT, PETSC_DEFAULT, maxits);  in SolGMRpSclr()
649         ierr = KSPSetFromOptions(ksps);   in SolGMRpSclr()
651       ierr = KSPSolve(ksps, resPs, DyPs);  in SolGMRpSclr()
653       ierr = VecScatterCreate(DyPs, LocalIndexSets, DyPLocals, NULL, &scatter7s);  in SolGMRpSclr()
655       ierr = VecScatterBegin(scatter7s, DyPs, DyPLocals, INSERT_VALUES, SCATTER_FORWARD);  in SolGMRpSclr()
656       ierr = VecScatterEnd(scatter7s, DyPs, DyPLocals, INSERT_VALUES, SCATTER_FORWARD);  in SolGMRpSclr()
657       ierr = VecGetArray(DyPLocals, &parray);  in SolGMRpSclr()
664       ierr = VecRestoreArray(DyPLocals, &parray);  in SolGMRpSclr()
672       ierr = KSPGetIterationNumber(ksps, &its);  in SolGMRpSclr()