| 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, Vec *initV, PetscInt size_initV, PetscInt plusk, PC pc,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk, PC pc,
|
HarmType_t harmMode, KSP ksp, InitType_t init, PetscTruth allResiduals)
|
HarmType_t harmMode, KSP ksp, InitType_t init, PetscTruth allResiduals)
|
{
|
{
|
PetscErrorCode ierr;
|
PetscErrorCode ierr;
|
PetscInt check_sum0, check_sum1;
|
PetscInt check_sum0, check_sum1;
|
|
|
| Line 49... |
Line 49... |
harmMode==DVD_HARM_NONE?PETSC_FALSE:PETSC_TRUE,
|
harmMode==DVD_HARM_NONE?PETSC_FALSE:PETSC_TRUE,
|
allResiduals);
|
allResiduals);
|
CHKERRQ(ierr);
|
CHKERRQ(ierr);
|
|
|
/* Setup the initial subspace for V */
|
/* Setup the initial subspace for V */
|
if (initV) {
|
if (size_initV) {
|
ierr = dvd_initV_user(d, b, initV, size_initV, ini_size_V); CHKERRQ(ierr);
|
ierr = dvd_initV_user(d, b, size_initV, ini_size_V); CHKERRQ(ierr);
|
} else switch(init) {
|
} else switch(init) {
|
case DVD_INITV_CLASSIC:
|
case DVD_INITV_CLASSIC:
|
ierr = dvd_initV_classic(d, b, ini_size_V); CHKERRQ(ierr); break;
|
ierr = dvd_initV_classic(d, b, ini_size_V); CHKERRQ(ierr); break;
|
case DVD_INITV_KRYLOV:
|
case DVD_INITV_KRYLOV:
|
ierr = dvd_initV_krylov(d, b, ini_size_V); CHKERRQ(ierr); break;
|
ierr = dvd_initV_krylov(d, b, ini_size_V); CHKERRQ(ierr); break;
|
}
|
}
|
|
|
/* Setup the convergence in order to use the SLEPc convergence test */
|
/* Setup the convergence in order to use the SLEPc convergence test */
|
ierr = dvd_testconv_slepc(d, b);CHKERRQ(ierr);
|
ierr = dvd_testconv_slepc(d, b);CHKERRQ(ierr);
|
|
|
/* Setup Raileigh-Ritz for selecting the best eigenpairs in V */
|
/* Setup Raileigh-Ritz for selecting the best eigenpairs in V */
|
| Line 87... |
Line 87... |
|
|
#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, Vec *initV, PetscInt size_initV, PetscInt plusk, PC pc,
|
PetscInt ini_size_V, PetscInt size_initV, PetscInt plusk, PC pc,
|
IP ip, HarmType_t harmMode, PetscTruth fixedTarget, PetscScalar t, KSP ksp,
|
IP ip, HarmType_t harmMode, PetscTruth fixedTarget, PetscScalar t, KSP ksp,
|
PetscReal fix, InitType_t init, PetscTruth allResiduals)
|
PetscReal fix, InitType_t init, PetscTruth allResiduals)
|
{
|
{
|
PetscInt check_sum0, check_sum1, maxits;
|
PetscInt check_sum0, check_sum1, maxits;
|
Vec *fv;
|
Vec *fv;
|
| Line 111... |
Line 111... |
harmMode==DVD_HARM_NONE?PETSC_FALSE:PETSC_TRUE,
|
harmMode==DVD_HARM_NONE?PETSC_FALSE:PETSC_TRUE,
|
allResiduals);
|
allResiduals);
|
CHKERRQ(ierr);
|
CHKERRQ(ierr);
|
|
|
/* Setup the initial subspace for V */
|
/* Setup the initial subspace for V */
|
if (initV) {
|
if (size_initV) {
|
ierr = dvd_initV_user(d, b, initV, size_initV, ini_size_V); CHKERRQ(ierr);
|
ierr = dvd_initV_user(d, b, size_initV, ini_size_V); CHKERRQ(ierr);
|
} else switch(init) {
|
} else switch(init) {
|
case DVD_INITV_CLASSIC:
|
case DVD_INITV_CLASSIC:
|
ierr = dvd_initV_classic(d, b, ini_size_V); CHKERRQ(ierr); break;
|
ierr = dvd_initV_classic(d, b, ini_size_V); CHKERRQ(ierr); break;
|
case DVD_INITV_KRYLOV:
|
case DVD_INITV_KRYLOV:
|
ierr = dvd_initV_krylov(d, b, ini_size_V); CHKERRQ(ierr); break;
|
ierr = dvd_initV_krylov(d, b, ini_size_V); CHKERRQ(ierr); break;
|