remove on_axis_J useless function

Author: beck-llr

src/ElectroMagn/ElectroMagnAM.cpp

@@ -1292,121 +1292,3 @@ void ElectroMagnAM::initAntennas( Patch *patch )
 
 }
 
-//! Evaluating EM fields modes correctly on axis
-void ElectroMagnAM::on_axis_J( bool diag_flag )
-{
-
-    if( isYmin ) {
-    
-        cField2D *Jl ;
-        cField2D *Jr ;
-        cField2D *Jt ;
-        
-        for( unsigned int imode=0 ; imode<nmodes ; imode++ ) {
-        
-            //static cast of the total currents and densities
-            Jl    = Jl_[imode];
-            Jr    = Jr_[imode];
-            Jt    = Jt_[imode];
-            
-            // Set Jr below axis to zero for all modes
-            for( unsigned int i=0; i<nl_p; i++ ) {
-                for( unsigned int j=0; j<oversize[1]+1; j++ ) {
-                    ( *Jr )( i, j ) = 0. ;
-                }
-            }
-            // Set Jt below and on axis to zero for all modes except mode 1
-            if( imode != 1 ) {
-                for( unsigned int i=0; i<nl_p; i++ ) {
-                    for( unsigned int j=0; j<oversize[1]+1; j++ ) {
-                        ( *Jt )( i, j ) = 0. ;
-                    }
-                }
-            } else {
-                // Set Jt below axis to zero for all modes
-                for( unsigned int i=0; i<nl_p; i++ ) {
-                    for( unsigned int j=0; j<oversize[1]; j++ ) {
-                        ( *Jt )( i, j ) = 0. ;
-                    }
-                    // Set Jt on axis for mode 1 for continuity equation
-                    ( *Jt )( i, oversize[1] ) = - 1./3.* ( 4.* Icpx * ( *Jr )( i, oversize[1]+1 ) + ( *Jt )( i, oversize[1]+1 ) );
-                }
-            }
-            if( imode != 0 ) {
-                // Set Jx below and on axis to zero for all modes except mode 0
-                for( unsigned int i=0; i<nl_d; i++ ) {
-                    for( unsigned int j=0; j<oversize[1]+1; j++ ) {
-                        ( *Jl )( i, j ) = 0. ;
-                    }
-                }
-            } else {
-                // Set Jx below axis to zero for all modes
-                for( unsigned int i=0; i<nl_d; i++ ) {
-                    for( unsigned int j=0; j<oversize[1]; j++ ) {
-                        ( *Jl )( i, j ) = 0. ;
-                    }
-                }
-                
-            }
-            
-        }
-        //if(diag_flag){
-        //    for ( unsigned int imode=0 ; imode<nmodes ; imode++ ) {
-        //        cField2D* rho   = rho_AM_[imode];
-        //        if (imode == 0){
-        //            for (unsigned int ism=0; ism < n_species; ism++){
-        //                Jt    = Jt_s[ism];
-        //                Jr    = Jr_s[ism];
-        //                if ( ( Jt != NULL )  )
-        //                    for (unsigned int i=0; i<nl_p; i++)
-        //                        (*Jt)(i,oversize[1]) = 0. ;
-        //                if ( ( Jr != NULL )  )
-        //                    for (unsigned int i=0; i<nl_p; i++)
-        //                        (*Jr)(i,oversize[1]) =  -(*Jr)(i,oversize[1]) ;
-        
-        //            }
-        //        }
-        //        else if (imode == 1){
-        //            for (unsigned int i=0; i<nl_p; i++)
-        //                (*rho)(i,oversize[1])= 0.;
-        //            //Loop on all modes and species for J_s
-        //            for (unsigned int ism=n_species; ism <  2*n_species; ism++){
-        //                Jl    = Jl_s[ism];
-        //                Jt    = Jt_s[ism];
-        //                Jr    = Jr_s[ism];
-        //                //if ( ( Jt != NULL ) && (Jr != NULL ) ) {
-        //                //    for (unsigned int i=0; i<nl_p; i++)
-        //                //        (*Jt)(i,oversize[1]) = - 1./3.* (4.* Icpx * (*Jr)(i,oversize[1]+1) + (*Jt)(i,oversize[1]+1));
-        //                //    for (unsigned int i=0; i<nl_p; i++)
-        //                //        (*Jr)(i,oversize[1])= 2.*Icpx* (*Jt)(i,oversize[1])-(*Jr)(i,oversize[1]+1) ;
-        //                //}
-        //                if ( Jl != NULL )
-        //                    for (unsigned int i=0; i<nl_d; i++)
-        //                        (*Jl)(i,oversize[1])= 0. ;
-        //            }
-        //        }
-        //        else {  // imode > 1
-        //            //Loop on all modes and species for J_s
-        //            for (unsigned int ism=2*n_species; ism <  n_species*nmodes ; ism++){
-        //                Jt    = Jt_s[ism];
-        //                if ( Jt != NULL )
-        //                    for (unsigned int i=0; i<nl_p; i++)
-        //                        (*Jt)(i,oversize[1]) = 0.;
-        
-        //                Jr    = Jr_s[ism];
-        //                if ( Jr != NULL )
-        //                    for (unsigned int i=0; i<nl_p; i++)
-        //                        (*Jr)(i,oversize[1]) = -(*Jr)(i,oversize[1]+1) ;
-        
-        //                Jl    = Jl_s[ism];
-        //                if ( Jl != NULL )
-        //                    for (unsigned int i=0; i<nl_d; i++)
-        //                        (*Jl)(i,oversize[1])= 0. ;
-        //            }
-        //        }
-        //    }
-        //}
-    }
-    return;
-}
-

src/ElectroMagn/ElectroMagnAM.h

@@ -200,7 +200,6 @@ class ElectroMagnAM : public ElectroMagn
         return norm2;
     }
 
-    void on_axis_J( bool diag_flag );
     //! from smpi is ymax
     const bool isYmin;