<div dir="ltr">Hello all,<div><br></div><div>I'm trying to load my laplacian matrix into a fortran module, and i have implemented it and it works for the first iteration of laplacian solver, but when starts the second step the laplacian matrix object becomes corrupts and looks like it loses one of it's dimensions. </div><div><br></div><div>Can you help me understand whats happening?<br></div><div><br></div><div>The modules are attached, the error i get is the following, i bolded the lines where i detected corruption:</div><div><font face="monospace, monospace"><br></font></div><div><div><font face="monospace, monospace"> ucmsSeamount Entering MAIN loop.</font></div><div><font face="monospace, monospace"> RHS loaded, size:       213120 /      213120</font></div><div><font face="monospace, monospace"> <b>CSRMAt loaded, sizes:      213120 x      213120</b></font></div><div><font face="monospace, monospace">   8.39198399     s</font></div><div><font face="monospace, monospace"> solveP pass:           1      !Iteration number</font></div><div><font face="monospace, monospace"> RHS loaded, size:       213120 /      213120</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Invalid argument</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Wrong type of object: Parameter # 1</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Petsc Release Version 3.7.3, unknown </font></div><div><font face="monospace, monospace">[0]PETSC ERROR: ./ucmsSeamount                                                                                                                                                                                                                                                  �J� on a arch-linux2-c-debug named valera-HP-xw4600-Workstation by valera Fri Sep 23 10:27:21 2016</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-fblaslapack=1 --download-mpich=1 --download-ml=1</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #1 MatGetSize() line 6295 in /home/valera/v5PETSc/petsc/petsc/src/mat/interface/matrix.c</font></div><div><font face="monospace, monospace"><b> CSRMAt loaded, sizes:      213120 x           0</b></font></div><div><font face="monospace, monospace">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Invalid argument</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Wrong type of object: Parameter # 2</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Petsc Release Version 3.7.3, unknown </font></div><div><font face="monospace, monospace">[0]PETSC ERROR: ./ucmsSeamount                                                                                                                                                                                                                                                  �J� on a arch-linux2-c-debug named valera-HP-xw4600-Workstation by valera Fri Sep 23 10:27:21 2016</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-fblaslapack=1 --download-mpich=1 --download-ml=1</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #2 KSPSetOperators() line 531 in /home/valera/v5PETSc/petsc/petsc/src/ksp/ksp/interface/itcreate.c</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Nonconforming object sizes</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Preconditioner number of local rows -1 does not equal resulting vector number of rows 213120</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Petsc Release Version 3.7.3, unknown </font></div><div><font face="monospace, monospace">[0]PETSC ERROR: ./ucmsSeamount                                                                                                                                                                                                                                                  �J� on a arch-linux2-c-debug named valera-HP-xw4600-Workstation by valera Fri Sep 23 10:27:21 2016</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-fblaslapack=1 --download-mpich=1 --download-ml=1</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #3 PCApply() line 474 in /home/valera/v5PETSc/petsc/petsc/src/ksp/pc/interface/precon.c</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Object is in wrong state</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Mat object's type is not set: Argument # 1</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Petsc Release Version 3.7.3, unknown </font></div><div><font face="monospace, monospace">[0]PETSC ERROR: ./ucmsSeamount                                                                                                                                                                                                                                                  �J� on a arch-linux2-c-debug named valera-HP-xw4600-Workstation by valera Fri Sep 23 10:27:21 2016</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-fblaslapack=1 --download-mpich=1 --download-ml=1</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #4 MatGetFactorAvailable() line 4286 in /home/valera/v5PETSc/petsc/petsc/src/mat/interface/matrix.c</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #5 PCGetDefaultType_Private() line 28 in /home/valera/v5PETSc/petsc/petsc/src/ksp/pc/interface/precon.c</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #6 PCSetFromOptions() line 159 in /home/valera/v5PETSc/petsc/petsc/src/ksp/pc/interface/pcset.c</font></div><div><font face="monospace, monospace">[0]PETSC ERROR: #7 KSPSetFromOptions() line 400 in /home/valera/v5PETSc/petsc/petsc/src/ksp/ksp/interface/itcl.c</font></div><div><font face="monospace, monospace">application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0</font></div><div><font face="monospace, monospace">[unset]: aborting job:</font></div><div><font face="monospace, monospace">application called MPI_Abort(MPI_COMM_WORLD, 73) - process 0</font></div></div><div><br></div></div>