/[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 587 by gross, Fri Mar 10 02:26:50 2006 UTC revision 588 by gross, Fri Mar 10 04:45:04 2006 UTC
# Line 177  void  vectorInKernel3__nonZeroA00(const Line 177  void  vectorInKernel3__nonZeroA00(const
177      register const double I00=1./A00;      register const double I00=1./A00;
178      register const double IA10=I00*A10;      register const double IA10=I00*A10;
179      register const double IA20=I00*A20;      register const double IA20=I00*A20;
180      vectorInKernel2(A11-IA10*A01,A21-IA20*A01,A12-IA10*A02,A22-IA20*A02,&TEMP0,&TEMP1);      vectorInKernel2(A11-IA10*A01,A12-IA10*A02,
181                        A21-IA20*A01,A22-IA20*A02,&TEMP0,&TEMP1);
182      *V0=-(A10*TEMP0+A20*TEMP1);      *V0=-(A10*TEMP0+A20*TEMP1);
183      *V1=A00*TEMP0;      *V1=A00*TEMP0;
184      *V2=A00*TEMP1;      *V2=A00*TEMP1;
# Line 331  void  eigenvalues_and_eigenvectors3(cons Line 332  void  eigenvalues_and_eigenvectors3(cons
332              *ev1=TEMP_EV0;              *ev1=TEMP_EV0;
333              *ev2=TEMP_EV1;              *ev2=TEMP_EV1;
334          } else if (A00>TEMP_EV1) {          } else if (A00>TEMP_EV1) {
335              *V00=TEMP_V00;              *V02=1.;
             *V10=TEMP_V10;  
             *V20=0.;  
             *V01=TEMP_V01;  
             *V11=TEMP_V11;  
             *V21=0.;  
             *V02=0.;  
336              *V12=0.;              *V12=0.;
337              *V22=1.;              *V22=0.;
338                *V00=0.;
339                *V10=TEMP_V00;
340                *V20=TEMP_V10;
341                *V01=0.;
342                *V11=TEMP_V01;
343                *V21=TEMP_V11;
344              *ev0=TEMP_EV0;              *ev0=TEMP_EV0;
345              *ev1=TEMP_EV1;              *ev1=TEMP_EV1;
346              *ev0=A00;              *ev2=A00;
347          } else {          } else {
348              *V00=TEMP_V00;              *V01=1.;
349              *V10=0;              *V11=0.;
             *V20=TEMP_V10;  
             *V01=0.;  
             *V11=1.;  
350              *V21=0.;              *V21=0.;
351              *V02=TEMP_V01;              *V00=0.;
352              *V12=0.;              *V10=TEMP_V00;
353                *V20=TEMP_V10;
354                *V02=0.;
355                *V12=TEMP_V01;
356              *V22=TEMP_V11;              *V22=TEMP_V11;
357              *ev0=TEMP_EV0;              *ev0=TEMP_EV0;
358              *ev1=A00;              *ev1=A00;

Legend:
Removed from v.587  
changed lines
  Added in v.588

  ViewVC Help
Powered by ViewVC 1.1.26