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

Diff of /trunk/ripley/src/Rectangle.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 3376  void Rectangle::assemblePDEBoundarySingl Line 3376  void Rectangle::assemblePDEBoundarySingl
3376      {      {
3377          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
3378              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
3379  #pragma omp for nowait  #pragma omp for
3380                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
3381                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3382                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3438  void Rectangle::assemblePDEBoundarySingl Line 3438  void Rectangle::assemblePDEBoundarySingl
3438    
3439          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
3440              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring                          for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring            
3441  #pragma omp for nowait  #pragma omp for
3442                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
3443                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3444                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3500  void Rectangle::assemblePDEBoundarySingl Line 3500  void Rectangle::assemblePDEBoundarySingl
3500    
3501          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
3502              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
3503  #pragma omp for nowait  #pragma omp for
3504                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
3505                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3506                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3562  void Rectangle::assemblePDEBoundarySingl Line 3562  void Rectangle::assemblePDEBoundarySingl
3562    
3563          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
3564              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
3565  #pragma omp for nowait  #pragma omp for
3566                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
3567                      const index_t e = m_faceOffset[3]+k0;                      const index_t e = m_faceOffset[3]+k0;
3568                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
# Line 3641  void Rectangle::assemblePDEBoundarySingl Line 3641  void Rectangle::assemblePDEBoundarySingl
3641      {      {
3642          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
3643              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
3644  #pragma omp for nowait  #pragma omp for
3645                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
3646                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3647                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3675  void Rectangle::assemblePDEBoundarySingl Line 3675  void Rectangle::assemblePDEBoundarySingl
3675    
3676          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
3677              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring                          for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring            
3678  #pragma omp for nowait  #pragma omp for
3679                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
3680                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3681                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3709  void Rectangle::assemblePDEBoundarySingl Line 3709  void Rectangle::assemblePDEBoundarySingl
3709    
3710          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
3711              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
3712  #pragma omp for nowait  #pragma omp for
3713                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
3714                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3715                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3742  void Rectangle::assemblePDEBoundarySingl Line 3742  void Rectangle::assemblePDEBoundarySingl
3742    
3743          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
3744              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
3745  #pragma omp for nowait  #pragma omp for
3746                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
3747                      vector<double> EM_S(4*4, 0);                      vector<double> EM_S(4*4, 0);
3748                      vector<double> EM_F(4, 0);                      vector<double> EM_F(4, 0);
# Line 3811  void Rectangle::assemblePDEBoundarySyste Line 3811  void Rectangle::assemblePDEBoundarySyste
3811      {      {
3812          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
3813              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
3814  #pragma omp for nowait  #pragma omp for
3815                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
3816                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
3817                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 3886  void Rectangle::assemblePDEBoundarySyste Line 3886  void Rectangle::assemblePDEBoundarySyste
3886    
3887          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
3888              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring                          for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring            
3889  #pragma omp for nowait  #pragma omp for
3890                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
3891                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
3892                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 3961  void Rectangle::assemblePDEBoundarySyste Line 3961  void Rectangle::assemblePDEBoundarySyste
3961    
3962          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
3963              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
3964  #pragma omp for nowait  #pragma omp for
3965                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
3966                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
3967                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 4036  void Rectangle::assemblePDEBoundarySyste Line 4036  void Rectangle::assemblePDEBoundarySyste
4036    
4037          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
4038              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
4039  #pragma omp for nowait  #pragma omp for
4040                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
4041                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
4042                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 4135  void Rectangle::assemblePDEBoundarySyste Line 4135  void Rectangle::assemblePDEBoundarySyste
4135      {      {
4136          if (m_faceOffset[0] > -1) {          if (m_faceOffset[0] > -1) {
4137              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring
4138  #pragma omp for nowait  #pragma omp for
4139                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
4140                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
4141                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 4176  void Rectangle::assemblePDEBoundarySyste Line 4176  void Rectangle::assemblePDEBoundarySyste
4176    
4177          if (m_faceOffset[1] > -1) {          if (m_faceOffset[1] > -1) {
4178              for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring                          for (index_t k1_0=0; k1_0<2; k1_0++) { // colouring            
4179  #pragma omp for nowait  #pragma omp for
4180                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {                  for (index_t k1=k1_0; k1<m_NE1; k1+=2) {
4181                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
4182                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 4217  void Rectangle::assemblePDEBoundarySyste Line 4217  void Rectangle::assemblePDEBoundarySyste
4217    
4218          if (m_faceOffset[2] > -1) {          if (m_faceOffset[2] > -1) {
4219              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
4220  #pragma omp for nowait  #pragma omp for
4221                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
4222                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
4223                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);
# Line 4258  void Rectangle::assemblePDEBoundarySyste Line 4258  void Rectangle::assemblePDEBoundarySyste
4258    
4259          if (m_faceOffset[3] > -1) {          if (m_faceOffset[3] > -1) {
4260              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring              for (index_t k0_0=0; k0_0<2; k0_0++) { // colouring
4261  #pragma omp for nowait  #pragma omp for
4262                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {                  for (index_t k0 = k0_0; k0 < m_NE0; k0+=2) {
4263                      vector<double> EM_S(4*4*numEq*numComp, 0);                      vector<double> EM_S(4*4*numEq*numComp, 0);
4264                      vector<double> EM_F(4*numEq, 0);                      vector<double> EM_F(4*numEq, 0);

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

  ViewVC Help
Powered by ViewVC 1.1.26