Subversion Repositories slepc-dev

Rev

Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 982 Rev 1009
Line 277... Line 277... 
    }
 
    }
 
    else {
 
    else {
 
      if (W) w = W[i];
 
      if (W) w = W[i];
 
      else w = V[i];
 
      else w = V[i];
 
    }
 
    }
 
    ierr = VecMDot(nv,w,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
 
#ifndef PETSC_USE_COMPLEX
 
        ierr = PetscPrintf(comm," %12g  ",vals[j]);CHKERRQ(ierr);
 
        ierr = PetscPrintf(comm," %12g  ",vals[j]);CHKERRQ(ierr);