| Line 64... |
Line 64... |
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 method for improving the eigenvectors */
|
/* Setup the method for improving the eigenvectors */
|
ierr = dvd_improvex_jd(d, b, ksp, bs, cX_proj); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd(d, b, ksp, bs, cX_proj, PETSC_FALSE); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_PROJ_KZX); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_PROJ_KZX); 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);
|
}
|
}
|
|
|
PetscFunctionReturn(0);
|
PetscFunctionReturn(0);
|
| Line 80... |
Line 80... |
PetscErrorCode dvd_schm_basic_conf(dvdDashboard *d, dvdBlackboard *b,
|
PetscErrorCode dvd_schm_basic_conf(dvdDashboard *d, dvdBlackboard *b,
|
PetscInt mpd, PetscInt min_size_V, PetscInt bs,
|
PetscInt mpd, PetscInt min_size_V, PetscInt bs,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk,
|
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, orthoV_type_t orth,
|
PetscReal fix, InitType_t init, PetscBool allResiduals, orthoV_type_t orth,
|
PetscInt cX_proj, PetscInt cX_impr)
|
PetscInt cX_proj, PetscInt cX_impr, PetscBool dynamic)
|
{
|
{
|
PetscInt check_sum0, check_sum1, maxits;
|
PetscInt check_sum0, check_sum1, maxits;
|
Vec *fv;
|
Vec *fv;
|
PetscScalar *fs;
|
PetscScalar *fs;
|
PetscReal tol;
|
PetscReal tol;
|
| Line 116... |
Line 116... |
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 method for improving the eigenvectors */
|
/* Setup the method for improving the eigenvectors */
|
ierr = dvd_improvex_jd(d, b, ksp, bs, cX_proj); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd(d, b, ksp, bs, cX_proj, dynamic); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_PROJ_KZX);
|
ierr = dvd_improvex_jd_proj_uv(d, b, DVD_PROJ_KZX);
|
CHKERRQ(ierr);
|
CHKERRQ(ierr);
|
ierr = KSPGetTolerances(ksp, &tol, PETSC_NULL, PETSC_NULL, &maxits);
|
ierr = KSPGetTolerances(ksp, &tol, PETSC_NULL, PETSC_NULL, &maxits);
|
CHKERRQ(ierr);
|
CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_lit_const(d, b, maxits, tol, fix); CHKERRQ(ierr);
|
ierr = dvd_improvex_jd_lit_const(d, b, maxits, tol, fix); CHKERRQ(ierr);
|