Topic: SLEPC integration

Hello Everyone. I have been trying to use SLEPC with OOFEM and compiled the program accordingly. However, when I use SLEPC as the eigen solver and smtype 7, I seem to get this error. Can anyone help on how to resolve this problem?

Assembling stiffness and mass matrices
Solving ...
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Shift-and-invert requires a target 'which' (see EPSSetWhichEigenpairs), for instance -st_type sinvert -eps_target 0 -eps_target_magnitude
[0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Development GIT revision: v3.14.5-729-g5947d95822  GIT Date: 2021-03-08 14:49:56 -0600
[0]PETSC ERROR: /home/deniz/build/debug/oofem on a arch-linux-c-debug named debian by deniz Wed Mar 24 09:14:46 2021
[0]PETSC ERROR: Configure options 
[0]PETSC ERROR: #1 EPSCheckCompatibleST() line 139 in /home/deniz/Programs/slepc/src/eps/interface/epssetup.c
[0]PETSC ERROR: #2 EPSSetUp() line 355 in /home/deniz/Programs/slepc/src/eps/interface/epssetup.c
[0]PETSC ERROR: #3 EPSSolve() line 136 in /home/deniz/Programs/slepc/src/eps/interface/epssolve.c
[0]PETSC ERROR: #4 solve() line 154 in /home/deniz/Programs/oofem/src/oofemlib/slepcsolver.C
EngngModel info: user time consumed by solution: 0.00s

 


ANALYSIS FINISHED

2

Re: SLEPC integration

Hi,
you should simply use the suggested options when invoking oofem (as they are passed to slepc):

./oofem -f eigen_beam2d.in -st_type sinvert -eps_target 0 -eps_target_magnitude

For the slepc options please follow slepc documentation.
Hope it will help,
Borek