/[escript]/trunk/escript/src/LocalOps.h
ViewVC logotype

Diff of /trunk/escript/src/LocalOps.h

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1945 by ksteube, Thu Sep 25 23:11:13 2008 UTC revision 1946 by jfenwick, Wed Oct 29 05:48:53 2008 UTC
# Line 388  void  eigenvalues_and_eigenvectors3(cons Line 388  void  eigenvalues_and_eigenvectors3(cons
388           } else {           } else {
389              const register double S00=A00-(*ev0);              const register double S00=A00-(*ev0);
390              const register double absS00=fabs(S00);              const register double absS00=fabs(S00);
391              if (fabs(S00)>m) {              if (absS00>m) {
392                  vectorInKernel3__nonZeroA00(S00,A01,A02,A01,A11-(*ev0),A12,A02,A12,A22-(*ev0),V00,V10,V20);                  vectorInKernel3__nonZeroA00(S00,A01,A02,A01,A11-(*ev0),A12,A02,A12,A22-(*ev0),V00,V10,V20);
393              } else if (absA02<m) {              } else if (absA02<m) {
394                  vectorInKernel3__nonZeroA00(A01,A11-(*ev0),A12,S00,A01,A02,A02,A12,A22-(*ev0),V00,V10,V20);                  vectorInKernel3__nonZeroA00(A01,A11-(*ev0),A12,S00,A01,A02,A02,A12,A22-(*ev0),V00,V10,V20);
# Line 398  void  eigenvalues_and_eigenvectors3(cons Line 398  void  eigenvalues_and_eigenvectors3(cons
398              normalizeVector3(V00,V10,V20);;              normalizeVector3(V00,V10,V20);;
399              const register double T00=A00-(*ev2);              const register double T00=A00-(*ev2);
400              const register double absT00=fabs(T00);              const register double absT00=fabs(T00);
401              if (fabs(T00)>m) {              if (absT00>m) {
402                   vectorInKernel3__nonZeroA00(T00,A01,A02,A01,A11-(*ev2),A12,A02,A12,A22-(*ev2),V02,V12,V22);                   vectorInKernel3__nonZeroA00(T00,A01,A02,A01,A11-(*ev2),A12,A02,A12,A22-(*ev2),V02,V12,V22);
403              } else if (absA02<m) {              } else if (absA02<m) {
404                   vectorInKernel3__nonZeroA00(A01,A11-(*ev2),A12,T00,A01,A02,A02,A12,A22-(*ev2),V02,V12,V22);                   vectorInKernel3__nonZeroA00(A01,A11-(*ev2),A12,T00,A01,A02,A02,A12,A22-(*ev2),V02,V12,V22);

Legend:
Removed from v.1945  
changed lines
  Added in v.1946

  ViewVC Help
Powered by ViewVC 1.1.26