#include "slepcbv.h" PetscErrorCode BVTraceQuadrature(BV Y,BV V,PetscInt L,PetscInt L_max,PetscScalar *w,VecScatter scat,PetscSubcomm subcomm,PetscInt npoints,PetscBool useconj,PetscReal *est_eig)Collective
Y | - first basis vectors | |
V | - second basis vectors | |
L | - block size | |
L_max | - maximum block size | |
w | - quadrature weights | |
scat | - (optional) VecScatter object to communicate between subcommunicators | |
subcomm | - subcommunicator layout | |
npoints | - number of points to process by the subcommunicator | |
useconj | - whether conjugate points can be used or not |
est_eig | - estimated eigenvalue count |
When using subcommunicators, Y is stored in the subcommunicators for a subset of integration points. In that case, the computation is done in the subcomm and then scattered to the whole communicator in S using the VecScatter scat. The value npoints is the number of points to be processed in this subcomm and the flag useconj indicates whether symmetric points can be reused.