| Line 27... |
Line 27... |
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "dvd_schm_basic_preconf"
|
#define __FUNCT__ "dvd_schm_basic_preconf"
|
PetscErrorCode dvd_schm_basic_preconf(dvdDashboard *d, dvdBlackboard *b,
|
PetscErrorCode dvd_schm_basic_preconf(dvdDashboard *d, dvdBlackboard *b,
|
PetscInt max_size_V, PetscInt mpd, PetscInt min_size_V, PetscInt bs,
|
PetscInt max_size_V, PetscInt mpd, PetscInt min_size_V, PetscInt bs,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk, PC pc,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk,
|
HarmType_t harmMode, KSP ksp, InitType_t init, PetscBool allResiduals)
|
HarmType_t harmMode, KSP ksp, InitType_t init, PetscBool allResiduals)
|
{
|
{
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscInt check_sum0, check_sum1;
|
PetscInt check_sum0, check_sum1;
|
|
|
| Line 67... |
Line 67... |
ierr = dvd_calcpairs_qz(d, b, PETSC_NULL); CHKERRQ(ierr);
|
ierr = dvd_calcpairs_qz(d, b, PETSC_NULL); CHKERRQ(ierr);
|
if (harmMode != DVD_HARM_NONE) {
|
if (harmMode != DVD_HARM_NONE) {
|
ierr = dvd_harm_conf(d, b, harmMode, PETSC_FALSE, 0.0); CHKERRQ(ierr);
|
ierr = dvd_harm_conf(d, b, harmMode, PETSC_FALSE, 0.0); CHKERRQ(ierr);
|
}
|
}
|
|
|
/* Setup the preconditioner */
|
|
ierr = dvd_static_precond_PC(d, b, pc); CHKERRQ(ierr);
|
|
|
|
/* Setup the method for improving the eigenvectors */
|
/* Setup the method for improving the eigenvectors */
|
ierr = dvd_improvex_jd(d, b, ksp, bs); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd(d, b, ksp, bs); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_PROJ_KBXX); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_PROJ_KBXX); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_lit_const(d, b, 0, 0.0, 0.0); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_lit_const(d, b, 0, 0.0, 0.0); CHKERRQ(ierr);
|
}
|
}
|
| Line 84... |
Line 81... |
|
|
#undef __FUNCT__
|
#undef __FUNCT__
|
#define __FUNCT__ "dvd_schm_basic_conf"
|
#define __FUNCT__ "dvd_schm_basic_conf"
|
PetscErrorCode dvd_schm_basic_conf(dvdDashboard *d, dvdBlackboard *b,
|
PetscErrorCode dvd_schm_basic_conf(dvdDashboard *d, dvdBlackboard *b,
|
PetscInt max_size_V, PetscInt mpd, PetscInt min_size_V, PetscInt bs,
|
PetscInt max_size_V, PetscInt mpd, PetscInt min_size_V, PetscInt bs,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk, PC pc,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk,
|
IP ip, HarmType_t harmMode, PetscBool fixedTarget, PetscScalar t, KSP ksp,
|
IP ip, HarmType_t harmMode, PetscBool fixedTarget, PetscScalar t, KSP ksp,
|
PetscReal fix, InitType_t init, PetscBool allResiduals)
|
PetscReal fix, InitType_t init, PetscBool allResiduals)
|
{
|
{
|
PetscInt check_sum0, check_sum1, maxits;
|
PetscInt check_sum0, check_sum1, maxits;
|
Vec *fv;
|
Vec *fv;
|
| Line 125... |
Line 122... |
/* Setup Raileigh-Ritz for selecting the best eigenpairs in V */
|
/* Setup Raileigh-Ritz for selecting the best eigenpairs in V */
|
ierr = dvd_calcpairs_qz(d, b, ip); CHKERRQ(ierr);
|
ierr = dvd_calcpairs_qz(d, b, ip); CHKERRQ(ierr);
|
if (harmMode != DVD_HARM_NONE) {
|
if (harmMode != DVD_HARM_NONE) {
|
ierr = dvd_harm_conf(d, b, harmMode, fixedTarget, t); CHKERRQ(ierr);
|
ierr = dvd_harm_conf(d, b, harmMode, fixedTarget, t); CHKERRQ(ierr);
|
}
|
}
|
|
|
/* Setup the preconditioner */
|
|
ierr = dvd_static_precond_PC(d, b, pc); CHKERRQ(ierr);
|
|
|
|
/* Setup the method for improving the eigenvectors */
|
/* Setup the method for improving the eigenvectors */
|
ierr = dvd_improvex_jd(d, b, ksp, bs); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd(d, b, ksp, bs); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_IS(d->sEP, DVD_EP_HERMITIAN)?
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_IS(d->sEP, DVD_EP_HERMITIAN)?
|
DVD_PROJ_KBXZ:DVD_PROJ_KBXY);
|
DVD_PROJ_KBXZ:DVD_PROJ_KBXY);
|