/[escript]/trunk/ripley/src/WaveAssembler2D.cpp
ViewVC logotype

Diff of /trunk/ripley/src/WaveAssembler2D.cpp

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

revision 4644 by sshaw, Fri Jan 24 03:29:25 2014 UTC revision 4645 by sshaw, Mon Feb 3 00:06:24 2014 UTC
# Line 516  const double X_10_0 = X_01_0, X_10_1 = X Line 516  const double X_10_0 = X_01_0, X_10_1 = X
516                              EM_F[INDEX2(1,2,numEq)]+=Btmp10 + Btmp11 + Btmp8 + Btmp9;                              EM_F[INDEX2(1,2,numEq)]+=Btmp10 + Btmp11 + Btmp8 + Btmp9;
517                              EM_F[INDEX2(1,3,numEq)]+=Btmp12 + Btmp13 + Btmp14 + Btmp15;                              EM_F[INDEX2(1,3,numEq)]+=Btmp12 + Btmp13 + Btmp14 + Btmp15;
518                          } else { // constant data                          } else { // constant data
519                              for (index_t k=0; k<numEq; k++) {  
520                                  throw RipleyException("WaveAssembler for constant data not yet implemented");                              const double wX_00 = -(du_p[INDEX2(0,0,numEq)] * c11_p[0]
521  //                                const double wX0 = X_p[INDEX2(k, 0, numEq)]*w18*6;                                                  + du_p[INDEX2(1,1,numEq)] * c13_p[0])*w18;
522  //                                const double wX1 = X_p[INDEX2(k, 1, numEq)]*w19*6;                              const double wX_01 = -(c44_p[0] *
523  //                                EM_F[INDEX2(k,0,numEq)]+= wX0 + wX1;                                                  (du_p[INDEX2(1,0,numEq)] + du_p[INDEX2(0,1,numEq)]))*w19;
524  //                                EM_F[INDEX2(k,1,numEq)]+=-wX0 + wX1;                              const double wX_10 = -(c44_p[0] *
525  //                                EM_F[INDEX2(k,2,numEq)]+= wX0 - wX1;                                                  (du_p[INDEX2(1,0,numEq)] + du_p[INDEX2(0,1,numEq)]))*w18;
526  //                                EM_F[INDEX2(k,3,numEq)]+=-wX0 - wX1;                              const double wX_11 = -(du_p[INDEX2(0,0,numEq)] * c13_p[0]
527                              }                                                  + du_p[INDEX2(1,1,numEq)] * c33_p[0])*w19;
528                                EM_F[INDEX2(0,0,numEq)]+= wX_00 + wX_01;
529                                EM_F[INDEX2(0,1,numEq)]+=-wX_00 + wX_01;
530                                EM_F[INDEX2(0,2,numEq)]+= wX_00 - wX_01;
531                                EM_F[INDEX2(0,3,numEq)]+=-wX_00 - wX_01;
532                                EM_F[INDEX2(1,0,numEq)]+= wX_10 + wX_11;
533                                EM_F[INDEX2(1,1,numEq)]+=-wX_10 + wX_11;
534                                EM_F[INDEX2(1,2,numEq)]+= wX_10 - wX_11;
535                                EM_F[INDEX2(1,3,numEq)]+=-wX_10 - wX_11;
536                          }                          }
537                      }                      }
538                      ///////////////                      ///////////////

Legend:
Removed from v.4644  
changed lines
  Added in v.4645

  ViewVC Help
Powered by ViewVC 1.1.26