/[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 586 by gross, Thu Mar 9 23:47:42 2006 UTC revision 587 by gross, Fri Mar 10 02:26:50 2006 UTC
# Line 132  void  vectorInKernel2(const double A00,c Line 132  void  vectorInKernel2(const double A00,c
132                        double* V0, double*V1)                        double* V0, double*V1)
133  {  {
134        register double absA00=fabs(A00);        register double absA00=fabs(A00);
       register double absA01=fabs(A01);  
135        register double absA10=fabs(A10);        register double absA10=fabs(A10);
136          register double absA01=fabs(A01);
137        register double absA11=fabs(A11);        register double absA11=fabs(A11);
138        register double m=absA11>absA01 ? absA11 : absA01;        register double m=absA11>absA10 ? absA11 : absA10;
139        if (absA00>m || absA10>m) {        if (absA00>m || absA01>m) {
140           *V0=-A10;           *V0=-A01;
141           *V1=A00;           *V1=A00;
142        } else {        } else {
143           if (m<=0) {           if (m<=0) {
# Line 145  void  vectorInKernel2(const double A00,c Line 145  void  vectorInKernel2(const double A00,c
145             *V1=0.;             *V1=0.;
146           } else {           } else {
147             *V0=A11;             *V0=A11;
148             *V1=-A01;             *V1=-A10;
149           }           }
150       }       }
151  }  }
# Line 223  void  eigenvalues_and_eigenvectors2(cons Line 223  void  eigenvalues_and_eigenvectors2(cons
223              *V00=-TEMP0*scale;              *V00=-TEMP0*scale;
224              *V10=-TEMP1*scale;              *V10=-TEMP1*scale;
225              if (TEMP1<0.) {              if (TEMP1<0.) {
226                 *V01=  *V10;                 *V01=  *V10;
227                 *V11=-(*V00);                 *V11=-(*V00);
228              } else {              } else {
229                 *V01=-(*V10);                 *V01=-(*V10);
230                 *V11= (*V00);                 *V11= (*V10);
231              }              }
232          } else if (TEMP0>0.) {          } else if (TEMP0>0.) {
233              *V00=TEMP0*scale;              *V00=TEMP0*scale;
234              *V10=TEMP1*scale;              *V10=TEMP1*scale;
235              if (TEMP1<0.) {              if (TEMP1<0.) {
236                 *V01=-(*V10);                 *V01=-(*V10);
237                 *V11= (*V00);                 *V11= (*V00);
238              } else {              } else {
239                 *V01= (*V10);                 *V01= (*V10);
240                 *V11=-(*V00);                 *V11=-(*V00);
241              }              }
242          } else {          } else {
243             *V00=0.;             *V00=0.;
Legend:
Removed from v.586  
changed lines
  Added in v.587
  ViewVC Help
Powered by ViewVC 1.1.26