« previous next »
Pages: [1]

Author Topic: ParFEAP for Contact Problems  (Read 9189 times)

ad10

  • Jr. Member
  • **
  • Posts: 39
ParFEAP for Contact Problems
« on: February 20, 2024, 08:26:21 PM »
Dear all,

I am trying to simulate a contact problem using ParFEAP, but the solution is not converging. Does FEAP not support in parallel for contact problems?

I used mpirun -n 1 $FEAPHOME8_6/parfeap/feap -ksp_type cg  -pc_type jacobi solver for my problem.

Code: [Select]
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Argument out of range
[0]PETSC ERROR: New nonzero at (30300,53434) 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.16.1, Nov 01, 2021
[0]PETSC ERROR: /home/ajinkya/Desktop/FEAP_Contact/ver86/parfeap/feap on a gnu-opt named ajinkya-HP-Pavilion-Gaming-Laptop-15-dk0xxx by ajinkya Tue Feb 20 22:20:01 2024
[0]PETSC ERROR: Configure options --download-metis --with-debugging=0
[0]PETSC ERROR: #1 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #2 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #3 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #4 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #5 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #6 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #7 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #8 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #9 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #10 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #11 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #12 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #13 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #14 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #15 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #16 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #17 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #18 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #19 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #20 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #21 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #22 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #23 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #24 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #25 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #26 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #27 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #28 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #29 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #30 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #31 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #32 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #33 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #34 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
[0]PETSC ERROR: #35 MatSetValues_SeqAIJ() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/impls/aij/seq/aij.c:520
[0]PETSC ERROR: #36 MatSetValues() at /home/ajinkya/FreeFEM/FreeFEM/FreeFem-sources/3rdparty/ff-petsc/petsc-3.16.1/src/mat/interface/matrix.c:1398
  NO CONVERGENCE REASON:  Indefinite Preconditioner
  NO CONVERGENCE REASON:  NaN in solution         
  NO CONVERGENCE REASON:  NaN in solution         
  NO CONVERGENCE REASON:  NaN in solution         
  NO CONVERGENCE REASON:  NaN in solution         
 Saving PARAVIEW data to Pexplicit_000100001.vtu                                                                                                         
  No history variables output to the PARAVIEW file
  NO CONVERGENCE REASON:  NaN in solution         
Serial code works fine. Do I need to change the solver type??
Logged

Prof. S. Govindjee

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 1164
Re: ParFEAP for Contact Problems
« Reply #1 on: February 21, 2024, 10:55:31 PM »
Contact in parallel should work in version 8.6.  Did you perform a graph partitioning for 1 partition?
Logged
Pages: [1]
« previous next »