/[escript]/branches/trilinos_from_5897/dudley/src/Assemble_PDE_System2_2D.cpp
ViewVC logotype

Diff of /branches/trilinos_from_5897/dudley/src/Assemble_PDE_System2_2D.cpp

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

revision 4346 by jfenwick, Tue Apr 2 04:46:45 2013 UTC revision 4521 by jfenwick, Mon Aug 26 11:51:30 2013 UTC
# Line 56  void Dudley_Assemble_PDE_System2_2D(Dudl Line 56  void Dudley_Assemble_PDE_System2_2D(Dudl
56      index_t *row_index;      index_t *row_index;
57      register dim_t q, s, r, k, m;      register dim_t q, s, r, k, m;
58      register double rtmp, rtmp0, rtmp1, rtmp00, rtmp10, rtmp01, rtmp11;      register double rtmp, rtmp0, rtmp1, rtmp00, rtmp10, rtmp01, rtmp11;
59      bool_t add_EM_F, add_EM_S;      bool add_EM_F, add_EM_S;
60    
61      bool_t extendedA = isExpanded(A);      bool extendedA = isExpanded(A);
62      bool_t extendedB = isExpanded(B);      bool extendedB = isExpanded(B);
63      bool_t extendedC = isExpanded(C);      bool extendedC = isExpanded(C);
64      bool_t extendedD = isExpanded(D);      bool extendedD = isExpanded(D);
65      bool_t extendedX = isExpanded(X);      bool extendedX = isExpanded(X);
66      bool_t extendedY = isExpanded(Y);      bool extendedY = isExpanded(Y);
67      double *F_p = (requireWrite(F), getSampleDataRW(F, 0)); /* use comma, to get around the mixed code and declarations thing */      double *F_p = (requireWrite(F), getSampleDataRW(F, 0)); /* use comma, to get around the mixed code and declarations thing */
68      const double *S = p.shapeFns;      const double *S = p.shapeFns;
69      dim_t len_EM_S = p.numShapes * p.numShapes * p.numEqu * p.numComp;      dim_t len_EM_S = p.numShapes * p.numShapes * p.numEqu * p.numComp;

Legend:
Removed from v.4346  
changed lines
  Added in v.4521

  ViewVC Help
Powered by ViewVC 1.1.26