[petsc-users] lock matrices in EPS

[Jose E. Roman]

We cannot reproduce this problem. Can you provide more details? For instance, a sequence of calls to SLEPc in your code. Are you modifying A and B after EPSSolve()?
Jose


> El 13 feb 2017, a las 23:33, Kong, Fande <fande.kong at inl.gov> escribió:
> 
> Hi ALL,
> 
> I am solving a generalized eigenvalue problem Ax=\lambda Bx. I want to retrieve matrices A and B from EPS, but got the following errors:
> 
> [0]PETSC ERROR: No support for this operation for this object type
> [0]PETSC ERROR: Cannot retrieve original matrices (have been modified)
> [0]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
> [0]PETSC ERROR: Petsc Release Version 3.7.5, unknown 
> [0]PETSC ERROR: ./ex14 on a arch-darwin-c-debug named FN604208 by kongf Mon Feb 13 15:29:33 2017
> [0]PETSC ERROR: Configure options --with-clanguage=cxx --with-shared-libraries=1 --with-blaslapack=1 --with-mpi=1 --download-parmetis=1 --download-metis=1 --with-debugging=yes --with-c2html=0 --download-hypre=1 --download-superlu_dist=1 PETSC_ARCH=arch-darwin-c-debug
> [0]PETSC ERROR: #1 STGetOperators() line 308 in /Users/kongf/projects/slepc/src/sys/classes/st/interface/stfunc.c
> [0]PETSC ERROR: #2 EPSGetOperators() line 319 in /Users/kongf/projects/slepc/src/eps/interface/epssetup.c
> 
> 
> My question is how to lock the original matrices I passed to EPS.  "-st_matmode copy" does not help.
> 
> Fande Kong,
>