/[escript]/trunk/finley/src/Util.cpp
ViewVC logotype

Diff of /trunk/finley/src/Util.cpp

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

revision 4788 by sshaw, Thu Feb 13 01:10:29 2014 UTC revision 4789 by sshaw, Fri Mar 21 00:30:32 2014 UTC
# Line 196  void normalVector(int len, int dim, int Line 196  void normalVector(int len, int dim, int
196                  A11=A[INDEX3(0,0,q,2,dim1)];                  A11=A[INDEX3(0,0,q,2,dim1)];
197                  A21=A[INDEX3(1,0,q,2,dim1)];                  A21=A[INDEX3(1,0,q,2,dim1)];
198                  length = sqrt(A11*A11+A21*A21);                  length = sqrt(A11*A11+A21*A21);
199                  if (! length>0) {                  if (length <= 0) {
200                      setError(ZERO_DIVISION_ERROR, __FILE__ ": area equals zero.");                      setError(ZERO_DIVISION_ERROR, __FILE__ ": area equals zero.");
201                      return;                      return;
202                  } else {                  } else {
# Line 218  void normalVector(int len, int dim, int Line 218  void normalVector(int len, int dim, int
218                  CO_A23=A31*A12-A11*A32;                  CO_A23=A31*A12-A11*A32;
219                  CO_A33=A11*A22-A21*A12;                  CO_A33=A11*A22-A21*A12;
220                  length=sqrt(CO_A13*CO_A13+CO_A23*CO_A23+CO_A33*CO_A33);                  length=sqrt(CO_A13*CO_A13+CO_A23*CO_A23+CO_A33*CO_A33);
221                  if (! length>0) {                  if (length <= 0) {
222                      setError(ZERO_DIVISION_ERROR, __FILE__ ": area equals zero.");                      setError(ZERO_DIVISION_ERROR, __FILE__ ": area equals zero.");
223                      return;                      return;
224                  } else {                  } else {

Legend:
Removed from v.4788  
changed lines
  Added in v.4789

  ViewVC Help
Powered by ViewVC 1.1.26