[petsc-users] SLEPc/NEP for shell matrice T(lambda) and T'(lambda)

[Jose E. Roman]

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 <kenneth.c.hall at duke.edu> 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 <jroman at dsic.upv.es>
> Date: Thursday, October 12, 2023 at 2:12 PM
> To: Kenneth C Hall <kenneth.c.hall at duke.edu>
> Cc: petsc-users at mcs.anl.gov <petsc-users at mcs.anl.gov>
> 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 <kenneth.c.hall at duke.edu> 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
> >