Hi John, Thanks. I'm able to successfully get the MatType in eigenproblems_ex3 after assembling the system and closing the PetscMatrix. I think the problem is either in the linking, or there is some memory corruption.
Thanks, Harshad On Mon, Oct 20, 2014 at 5:18 PM, John Peterson <[email protected]> wrote: > 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
