Kenneth, I have worked a bit more on your example and put it in SLEPc https://gitlab.com/slepc/slepc/-/merge_requests/596 This version also has MATOP_DESTROY to avoid memory leaks.
Thanks. Jose > El 12 oct 2023, a las 20:59, Kenneth C Hall <[email protected]> > escribió: > > Jose, > > Thanks very much for this. I will give it a try and let you know how it works. > > Best regards, > Kenneth > > From: Jose E. Roman <[email protected]> > Date: Thursday, October 12, 2023 at 2:12 PM > To: Kenneth C Hall <[email protected]> > Cc: [email protected] <[email protected]> > Subject: Re: [petsc-users] SLEPc/NEP for shell matrice T(lambda) and > T'(lambda) > > I am attaching your example modified with the context stuff. > > With the PETSc branch that I indicated, now it works with NLEIGS, for > instance: > > $ ./test_nep -nep_nleigs_ksp_type gmres -nep_nleigs_pc_type none > -rg_interval_endpoints 0.2,1.1 -nep_target 0.8 -nep_nev 5 -n 400 -nep_monitor > -nep_view -nep_error_relative ::ascii_info_detail > > And also other solvers such as SLP: > > $ ./test_nep -nep_type slp -nep_slp_ksp_type gmres -nep_slp_pc_type none > -nep_target 0.8 -nep_nev 5 -n 400 -nep_monitor -nep_error_relative > ::ascii_info_detail > > I will clean the example code an add it as a SLEPc example. > > Regards, > Jose > > > > El 11 oct 2023, a las 17:27, Kenneth C Hall <[email protected]> > > escribió: > > > > Jose, > > > > Thanks very much for your help with this. Greatly appreciated. I will look > > at the MR. Please let me know if you do get the Fortran example working. > > > > Thanks, and best regards, > > Kenneth > >
