/[escript]/branches/diaplayground/ripley/src/Brick.cpp
ViewVC logotype

Diff of /branches/diaplayground/ripley/src/Brick.cpp

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

revision 3793 by gross, Wed Feb 1 07:39:43 2012 UTC revision 3795 by caltinay, Thu Feb 2 05:54:34 2012 UTC
# Line 9438  void Brick::assemblePDEBoundarySingle(Pa Line 9438  void Brick::assemblePDEBoundarySingle(Pa
9438      {      {
9439          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
9440              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
9441  #pragma omp for nowait  #pragma omp for
9442                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
9443                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
9444                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 9562  void Brick::assemblePDEBoundarySingle(Pa Line 9562  void Brick::assemblePDEBoundarySingle(Pa
9562    
9563          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
9564              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
9565  #pragma omp for nowait  #pragma omp for
9566                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
9567                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
9568                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 9686  void Brick::assemblePDEBoundarySingle(Pa Line 9686  void Brick::assemblePDEBoundarySingle(Pa
9686    
9687          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
9688              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
9689  #pragma omp for nowait  #pragma omp for
9690                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
9691                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
9692                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 9810  void Brick::assemblePDEBoundarySingle(Pa Line 9810  void Brick::assemblePDEBoundarySingle(Pa
9810    
9811          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
9812              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
9813  #pragma omp for nowait  #pragma omp for
9814                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
9815                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
9816                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 9934  void Brick::assemblePDEBoundarySingle(Pa Line 9934  void Brick::assemblePDEBoundarySingle(Pa
9934    
9935          if (m_faceOffset[4] > -1) {          if (m_faceOffset[4] > -1) {
9936              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
9937  #pragma omp for nowait  #pragma omp for
9938                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
9939                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
9940                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10058  void Brick::assemblePDEBoundarySingle(Pa Line 10058  void Brick::assemblePDEBoundarySingle(Pa
10058    
10059          if (m_faceOffset[5] > -1) {          if (m_faceOffset[5] > -1) {
10060              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
10061  #pragma omp for nowait  #pragma omp for
10062                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
10063                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10064                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10202  void Brick::assemblePDEBoundarySingleRed Line 10202  void Brick::assemblePDEBoundarySingleRed
10202      {      {
10203          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
10204              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10205  #pragma omp for nowait  #pragma omp for
10206                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10207                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
10208                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10254  void Brick::assemblePDEBoundarySingleRed Line 10254  void Brick::assemblePDEBoundarySingleRed
10254    
10255          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
10256              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10257  #pragma omp for nowait  #pragma omp for
10258                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10259                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
10260                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10306  void Brick::assemblePDEBoundarySingleRed Line 10306  void Brick::assemblePDEBoundarySingleRed
10306    
10307          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
10308              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10309  #pragma omp for nowait  #pragma omp for
10310                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10311                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10312                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10357  void Brick::assemblePDEBoundarySingleRed Line 10357  void Brick::assemblePDEBoundarySingleRed
10357    
10358          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
10359              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10360  #pragma omp for nowait  #pragma omp for
10361                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10362                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10363                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10409  void Brick::assemblePDEBoundarySingleRed Line 10409  void Brick::assemblePDEBoundarySingleRed
10409    
10410          if (m_faceOffset[4] > -1) {          if (m_faceOffset[4] > -1) {
10411              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
10412  #pragma omp for nowait  #pragma omp for
10413                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
10414                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10415                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10461  void Brick::assemblePDEBoundarySingleRed Line 10461  void Brick::assemblePDEBoundarySingleRed
10461    
10462          if (m_faceOffset[5] > -1) {          if (m_faceOffset[5] > -1) {
10463              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
10464  #pragma omp for nowait  #pragma omp for
10465                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
10466                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10467                          vector<double> EM_S(8*8, 0);                          vector<double> EM_S(8*8, 0);
# Line 10570  void Brick::assemblePDEBoundarySystem(Pa Line 10570  void Brick::assemblePDEBoundarySystem(Pa
10570      {      {
10571          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
10572              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10573  #pragma omp for nowait  #pragma omp for
10574                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10575                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
10576                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 10708  void Brick::assemblePDEBoundarySystem(Pa Line 10708  void Brick::assemblePDEBoundarySystem(Pa
10708    
10709          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
10710              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10711  #pragma omp for nowait  #pragma omp for
10712                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10713                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
10714                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 10846  void Brick::assemblePDEBoundarySystem(Pa Line 10846  void Brick::assemblePDEBoundarySystem(Pa
10846    
10847          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
10848              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10849  #pragma omp for nowait  #pragma omp for
10850                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10851                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10852                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 10984  void Brick::assemblePDEBoundarySystem(Pa Line 10984  void Brick::assemblePDEBoundarySystem(Pa
10984    
10985          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
10986              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
10987  #pragma omp for nowait  #pragma omp for
10988                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
10989                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
10990                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11122  void Brick::assemblePDEBoundarySystem(Pa Line 11122  void Brick::assemblePDEBoundarySystem(Pa
11122    
11123          if (m_faceOffset[4] > -1) {          if (m_faceOffset[4] > -1) {
11124              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
11125  #pragma omp for nowait  #pragma omp for
11126                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
11127                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
11128                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11260  void Brick::assemblePDEBoundarySystem(Pa Line 11260  void Brick::assemblePDEBoundarySystem(Pa
11260    
11261          if (m_faceOffset[5] > -1) {          if (m_faceOffset[5] > -1) {
11262              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
11263  #pragma omp for nowait  #pragma omp for
11264                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
11265                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
11266                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11425  void Brick::assemblePDEBoundarySystemRed Line 11425  void Brick::assemblePDEBoundarySystemRed
11425      {      {
11426          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
11427              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
11428  #pragma omp for nowait  #pragma omp for
11429                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
11430                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
11431                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11483  void Brick::assemblePDEBoundarySystemRed Line 11483  void Brick::assemblePDEBoundarySystemRed
11483    
11484          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
11485              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
11486  #pragma omp for nowait  #pragma omp for
11487                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
11488                      for (index_t k1=0; k1<m_NE1; ++k1) {                      for (index_t k1=0; k1<m_NE1; ++k1) {
11489                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11541  void Brick::assemblePDEBoundarySystemRed Line 11541  void Brick::assemblePDEBoundarySystemRed
11541    
11542          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
11543              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
11544  #pragma omp for nowait  #pragma omp for
11545                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
11546                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
11547                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11599  void Brick::assemblePDEBoundarySystemRed Line 11599  void Brick::assemblePDEBoundarySystemRed
11599    
11600          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
11601              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring              for (index_t k2_0=0; k2_0<2; k2_0++) { // colouring
11602  #pragma omp for nowait  #pragma omp for
11603                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {                  for (index_t k2=k2_0; k2<m_NE2; k2+=2) {
11604                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
11605                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11657  void Brick::assemblePDEBoundarySystemRed Line 11657  void Brick::assemblePDEBoundarySystemRed
11657    
11658          if (m_faceOffset[4] > -1) {          if (m_faceOffset[4] > -1) {
11659              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
11660  #pragma omp for nowait  #pragma omp for
11661                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
11662                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
11663                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);
# Line 11715  void Brick::assemblePDEBoundarySystemRed Line 11715  void Brick::assemblePDEBoundarySystemRed
11715    
11716          if (m_faceOffset[5] > -1) {          if (m_faceOffset[5] > -1) {
11717              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
11718  #pragma omp for nowait  #pragma omp for
11719                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
11720                      for (index_t k0=0; k0<m_NE0; ++k0) {                      for (index_t k0=0; k0<m_NE0; ++k0) {
11721                          vector<double> EM_S(8*8*numEq*numComp, 0);                          vector<double> EM_S(8*8*numEq*numComp, 0);

Legend:
Removed from v.3793  
changed lines
  Added in v.3795

  ViewVC Help
Powered by ViewVC 1.1.26