On Mon, Oct 20, 2014 at 2:44 PM, Harshad Sahasrabudhe <[email protected]> wrote: > Yes, that's true. I'm trying to use an external eigensolver called FEAST, > instead of SLEPc. So instead of > > _system->solve(); > > I call > > assemble_eigensystem(*_equation_systems, _system_name); > _system->matrix_A->close(); > _system->matrix_B->close(); > > where assemble_eigensystem() is the function which assembles the system. > Then I call _system->matrix_A->mat() to get the PETSc matrix.
Assuming nothing strange is happening in assemble_eigensystem(), this seems like it should work. I don't have a SLEPc build handy, but have you tried modifying eigenproblems_ex3 to call mat() on matrix_A after assembly? If that works, I'd have to suspect something is wrong with the crazy library linking stuff you are doing... -- John ------------------------------------------------------------------------------ Comprehensive Server Monitoring with Site24x7. Monitor 10 servers for $9/Month. Get alerted through email, SMS, voice calls or mobile push notifications. Take corrective actions from your mobile device. http://p.sf.net/sfu/Zoho _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
