Merge branch 'develop' of llrgit.in2p3.fr:smilei/smilei into develop

SmileiPIC · Dec 15, 2023 · 7543fe8 · 7543fe8
2 parents 965c7f8 + 2dd6bb6
commit 7543fe8
Show file tree

Hide file tree

Showing 3 changed files with 2 additions and 18 deletions.
diff --git a/src/ElectroMagnBC/ElectroMagnBCAM_BM.cpp b/src/ElectroMagnBC/ElectroMagnBCAM_BM.cpp
@@ -188,23 +188,13 @@ void ElectroMagnBCAM_BM::apply( ElectroMagn *EMfields, double, Patch *patch )
 
             unsigned int j= n_d[1]-2;
 
-
-            //   MESSAGE("JGLOB "<< patch->getCellStartingGlobalIndex(1)+j);
-            //std::cout<<"come heree "<<patch->getCellStartingGlobalIndex(1)<<"  "<<j<<" \n " ;
-            //std::cout<<"come here "<<n_p[1] <<" nr*d[1] "<<n_p[1]*d[1]<<" \n " ;
             // for Bl^(p,d)
             for( unsigned int i=0 ; i<n_p[0]-1; i++ ) {
                 ( *Bl )( i, j+1 ) = ( *Bl_old )( i, j )
                                     -      Alpha_Bl_Rmax * ( ( *Bl )( i, j ) - ( *Bl_old )( i, j+1 ) )
                                     +      Gamma_Bl_Rmax * ( ( *Br )( i+1, j ) + ( *Br_old )( i+1, j ) - ( *Br )( i, j ) - ( *Br_old )( i, j ) )
                                     -      Beta_Bl_Rmax * Icpx * ( double )imode * ( ( *Er )( i, j+1 ) + ( *Er )( i, j ) )
                                     - 2. * Beta_Bl_Rmax * ( *Et )( i, j );
-                //if (std::abs((*Bl)(i,j+1))>1.){
-                //MESSAGE("BlBM");
-                //MESSAGE(i);
-                //MESSAGE(j+1);
-                //MESSAGE((*Bl)(i,j+1));
-                //}
             }//i  ---end Bl
 
             // for Bt^(d,d)
@@ -215,12 +205,6 @@ void ElectroMagnBCAM_BM::apply( ElectroMagn *EMfields, double, Patch *patch )
                                         + Gamma_Bt_Rmax * ( *Bt_old )( i, j )
                                         - Icpx * ( double )imode * CB_BM * Epsilon_Bt_Rmax  * ( ( *Br )( i, j ) + ( *Br_old )( i, j ) )
                                         - CE_BM * Delta_Bt_Rmax * ( ( *Er )( i, j+1 )+( *Er )( i, j )-( *Er )( i-1, j+1 ) -( *Er )( i-1, j ) ) ;
-                //if (std::abs((*Bt)(i,j+1))>1.){
-                //    MESSAGE("BtMF");
-                //    MESSAGE(i);
-                //    MESSAGE(j+1);
-                //    MESSAGE((*Bt)(i,j+1));
-                //}
             }//i  ---end Bt
         }
     }

diff --git a/src/ElectroMagnBC/ElectroMagnBCAM_SM.cpp b/src/ElectroMagnBC/ElectroMagnBCAM_SM.cpp
@@ -203,8 +203,8 @@ void ElectroMagnBCAM_SM::apply( ElectroMagn *EMfields, double time_dual, Patch *
                 unsigned int i=0;
                 ( *Br )( i, j ) = Alpha_Xmin   * ( *Et )( i, j )
                                   +              Beta_Xmin    * ( *Br )( i+1, j )
-                                  +              Gamma_Xmin   * byW;
-                +              Delta_Xmin   *( ( *Bl )( i, j+1 )- ( *Bl )( i, j ) );
+                                  +              Gamma_Xmin   * byW
+                                  +              Delta_Xmin   *( ( *Bl )( i, j+1 )- ( *Bl )( i, j ) );
             }//j  ---end compute Br
 
 

diff --git a/validation/references/tstAM_04b_laser_propagation_numpy.py.txt b/validation/references/tstAM_04b_laser_propagation_numpy.py.txt