Hello everyone, I hope this email finds you well.
My Name is Sawsan Shatanawi, and I am currently working on developing a Fortran code for simulating groundwater flow in a 3D system. The code involves solving a nonlinear system, and I have created the matrix to be solved using the PCG solver and Picard iteration. However, when I tried to assign it as a PETSc matrix I started getting a lot of error messages. I am kindly asking if someone can help me, I would be happy to share my code with him/her. Please find the attached file contains a list of errors I have gotten Thank you in advance for your time and assistance. Best regards, Sawsan
make: `GW.exe' is up to date. 5.2000000E-03 50.00000 10.00000 0.0000000E+00 [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: Argument out of range [0]PETSC ERROR: New nonzero at (1,1) caused a malloc Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting. [0]PETSC ERROR: Petsc Release Version 3.20.0, Sep 28, 2023 [0]PETSC ERROR: ./GW.exe on a named cn32 by sawsan.shatanawi Tue Dec 19 11:30:08 2023 [0]PETSC ERROR: Configure options --prefix=/opt/apps/petsc/3.20.0/intel/20.2/mpich/3.3.2 --with-blaslapack-dir=/opt/apps/intel/20.2/compilers_and_libraries_2020.2.254/linux/mkl --with-cc=mpicc --with-cxx=mpicxx --with-fc=mpif90 [0]PETSC ERROR: #1 MatSetValues_SeqAIJ() at /weka/data/kamiak/rohit.dhariwal/software_installs/petsc/3.20.0/petsc-3.20.0/src/mat/impls/aij/seq/aij.c:458 [0]PETSC ERROR: #2 MatSetValues() at /weka/data/kamiak/rohit.dhariwal/software_installs/petsc/3.20.0/petsc-3.20.0/src/mat/interface/matrix.c:1503 [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: It appears a new error in the code was triggered after a previous error, possibly because: [0]PETSC ERROR: - The first error was not properly handled via (for example) the use of [0]PETSC ERROR: PetscCall(TheFunctionThatErrors()); or [0]PETSC ERROR: - The second error was triggered while handling the first error. [0]PETSC ERROR: Above is the traceback for the previous unhandled error, below the traceback for the next error [0]PETSC ERROR: ALL ERRORS in the PETSc libraries are fatal, you should add the appropriate error checking to the code [0]PETSC ERROR: Argument out of range [0]PETSC ERROR: New nonzero at (1,2) caused a malloc Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check [0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting. [0]PETSC ERROR: Petsc Release Version 3.20.0, Sep 28, 2023 [0]PETSC ERROR: ./GW.exe on a named cn32 by sawsan.shatanawi Tue Dec 19 11:30:08 2023 [0]PETSC ERROR: Configure options --prefix=/opt/apps/petsc/3.20.0/intel/20.2/mpich/3.3.2 --with-blaslapack-dir=/opt/apps/intel/20.2/compilers_and_libraries_2020.2.254/linux/mkl --with-cc=mpicc --with-cxx=mpicxx --with-fc=mpif90 [0]PETSC ERROR: #1 MatSetValues_SeqAIJ() at /weka/data/kamiak/rohit.dhariwal/software_installs/petsc/3.20.0/petsc-3.20.0/src/mat/impls/aij/seq/aij.c:458 [0]PETSC ERROR: #2 MatSetValues() at /weka/data/kamiak/rohit.dhariwal/software_installs/petsc/3.20.0/petsc-3.20.0/src/mat/interface/matrix.c:1503 [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: It appears a new error in the code was triggered after a previous error, possibly because: [0]PETSC ERROR: - The first error was not properly handled via (for example) the use of [0]PETSC ERROR: PetscCall(TheFunctionThatErrors()); or [0]PETSC ERROR: - The second error was triggered while handling the first error. [0]PETSC ERROR: Above is the traceback for the previous unhandled error, below the traceback for the next error [0]PETSC ERROR: ALL ERRORS in the PETSc libraries are fatal, you should add the appropriate error checking to the code [0]PETSC ERROR: Argument out of range [0]PETSC ERROR: New nonzero at (1,16) caused a malloc