| Line 672... |
Line 672... |
}
|
}
|
ierr = VecMDot(w,nv,V,vals);CHKERRQ(ierr);
|
ierr = VecMDot(w,nv,V,vals);CHKERRQ(ierr);
|
for (j=0;j<nv;j++) {
|
for (j=0;j<nv;j++) {
|
if (lev) *lev += (j==i)? (vals[j]-1.0)*(vals[j]-1.0): vals[j]*vals[j];
|
if (lev) *lev += (j==i)? (vals[j]-1.0)*(vals[j]-1.0): vals[j]*vals[j];
|
else {
|
else {
|
#ifndef PETSC_USE_COMPLEX
|
#if !defined(PETSC_USE_COMPLEX)
|
ierr = PetscPrintf(comm," %12g ",vals[j]);CHKERRQ(ierr);
|
ierr = PetscPrintf(comm," %12g ",vals[j]);CHKERRQ(ierr);
|
#else
|
#else
|
ierr = PetscPrintf(comm," %12g%+12gi ",PetscRealPart(vals[j]),PetscImaginaryPart(vals[j]));CHKERRQ(ierr);
|
ierr = PetscPrintf(comm," %12g%+12gi ",PetscRealPart(vals[j]),PetscImaginaryPart(vals[j]));CHKERRQ(ierr);
|
#endif
|
#endif
|
}
|
}
|