Changeset 68 for branches/branche-mb/src/dta_lfexp.h

Timestamp:

08/30/11 16:13:04 (13 years ago)

Author:

berrada

Message:

commit final avec filtre

File:

• branches/branche-mb/src/dta_lfexp.h (modified) (4 diffs)

branches/branche-mb/src/dta_lfexp.h

@@ -37 +37 @@
         tim1= 0.;
 #ifdef K_NEU
-        tim1= - bim1 * x2;
+        tim1= bim1 * x2;
 #endif
       }
       else{
-        tim1= -bim1 * x3;
+        tim1= bim1 * x3;
       }
       if (Yi==NX-2){
         tip1= 0.;
 #ifdef K_NEU
-        tip1= - bi * x2;
+        tip1= bi * x2;
 #endif
       }
       else{
-        tip1= -bi * x4;
+        tip1= bi * x4;
       }
       //cdt au bord j
@@ -56 +56 @@
         tjm1= 0.;
 #ifdef K_NEU
-        tjm1= - bjm1 * x2;
+        tjm1= bjm1 * x2;
 #endif
       }
       else{
-        tjm1= -bjm1 * x5;
+        tjm1= bjm1 * x5;
       }
       if (Yj==NY-2){
         tjp1= 0.;
 #ifdef K_NEU
-        tjp1= - bj * x2;
+        tjp1= bj * x2;
 #endif
       }
       else{
-        tjp1= -bj * x6;
+        tjp1= bj * x6;
       }
 
-      tij= (1./zbtr+bi+bim1+bj+bjm1) * x2;
+      tij= (1./zbtr-bi-bim1-bj-bjm1) * x2;
       YS1 = zbtr*(tip1 + tim1 + tjp1 + tjm1 + tij);
     }
@@ -105 +105 @@
         dtim1_x2= 0.;
 #ifdef K_NEU
-        dtim1_x2= - bim1;
+        dtim1_x2= bim1;
 #endif
       }
       else{
-        dtim1_x3= -bim1;
+        dtim1_x3= bim1;
       }
       if (Yi==NX-2){
         dtip1_x2= 0.;
 #ifdef K_NEU
-        dtip1_x2= - bi;
+        dtip1_x2= bi;
 #endif
       }
       else{
-        dtip1_x4= -bi;
+        dtip1_x4= bi;
       }
       //cdt au bord j
@@ -124 +124 @@
         dtjm1_x2= 0.;
 #ifdef K_NEU
-        dtjm1_x2= - bjm1;
+        dtjm1_x2= bjm1;
 #endif
       }
       else{
-        dtjm1_x5= -bjm1;
+        dtjm1_x5= bjm1;
       }
       if (Yj==NY-2){
         dtjp1_x2= 0.;
 #ifdef K_NEU
-        dtjp1_x2= - bj;
+        dtjp1_x2= bj;
 #endif
       }
       else{
-        dtjp1_x6= -bj;
+        dtjp1_x6= bj;
       }
-      dtij_x2= (1./zbtr+bi+bim1+bj+bjm1);
+      dtij_x2= (1./zbtr-bi-bim1-bj-bjm1);
 
       YJ1I1 = 0.;