/trunk/src/ps/psgnhep.c
154,7 → 154,7
| #endif |
| if (info) SETERRQ1(((PetscObject)ps)->comm,PETSC_ERR_LIB,"Error in Lapack xGGES %i",info); |
| for (i=0;i<n;i++) { |
| if (beta[i]==0.0) wr[i] = (wr[i]>0.0)? PETSC_MAX_REAL: PETSC_MIN_REAL; |
| if (beta[i]==0.0) wr[i] = (PetscRealPart(wr[i])>0.0)? PETSC_MAX_REAL: PETSC_MIN_REAL; |
| else wr[i] /= beta[i]; |
| #if !defined(PETSC_USE_COMPLEX) |
| if (beta[i]==0.0) wi[i] = (wi[i]>0.0)? PETSC_MAX_REAL: PETSC_MIN_REAL; |
240,7 → 240,7
| } else |
| #endif |
| { |
| if (T[j*ld+j] == 0.0) wi[j] = (S[j*ld+j]>0.0)? PETSC_MAX_REAL: PETSC_MIN_REAL; |
| if (T[j*ld+j] == 0.0) wr[j] = (PetscRealPart(S[j*ld+j])>0.0)? PETSC_MAX_REAL: PETSC_MIN_REAL; |
| else wr[j] = S[j*ld+j] / T[j*ld+j]; |
| wi[j] = 0.0; |
| } |
/trunk/src/eps/impls/krylov/krylov.c
97,7 → 97,7
| ierr = PSGetLeadingDimension(eps->ps,&ld);CHKERRQ(ierr); |
| ierr = PetscObjectTypeCompare((PetscObject)eps->OP,STSHIFT,&isshift);CHKERRQ(ierr); |
| marker = -1; |
| getall = getall || eps->trackall; |
| if (eps->trackall) getall = PETSC_TRUE; |
| for (k=kini;k<kini+nits;k++) { |
| /* eigenvalue */ |
| re = eps->eigr[k]; |