adding normal only

Author: EyaDammak

Source/Particles/ParticleBoundaryBuffer.cpp

@@ -167,9 +167,6 @@ struct CopyAndTimestamp {
     int m_delta_index;
     int m_normal_index;
     int m_step;
-    const amrex::Real m_dt;
-    amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> m_plo;
-    amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> m_phi;
     int m_idim;
     int m_iside;
 
@@ -189,59 +186,16 @@ struct CopyAndTimestamp {
         for (int j = 0; j < src.m_num_runtime_int; ++j) {
             dst.m_runtime_idata[j][dst_i] = src.m_runtime_idata[j][src_i];
         }
-        
-        amrex::ParticleReal const ux = dst.m_rdata[PIdx::ux][dst_i];
-        amrex::ParticleReal const uy = dst.m_rdata[PIdx::uy][dst_i];
-        amrex::ParticleReal const uz = dst.m_rdata[PIdx::uz][dst_i];
 
         dst.m_runtime_idata[m_step_index][dst_i] = m_step;
-        amrex::Real delta_t;
-        amrex::Real exit_velocity = 0;
-        if (m_idim == 0) {
-            exit_velocity = ux;
-        } else if (m_idim == 1) {
-            exit_velocity = uy;
-        } else {
-            exit_velocity = uz;
-        }
+        dst.m_runtime_rdata[m_delta_index][dst_i] = 0._rt; 
 
-        if (m_iside==1){ //iside=1, we are in the right side
-            delta_t = std::abs(dst.m_rdata[PIdx::z][dst_i]-m_phi[m_idim])/exit_velocity; 
-            dst.m_runtime_rdata[m_delta_index][dst_i] = delta_t;
-        } else{
-            delta_t = std::abs(dst.m_rdata[PIdx::z][dst_i]-m_plo[m_idim])/exit_velocity; 
-            dst.m_runtime_rdata[m_delta_index][dst_i] = delta_t; 
-        }
-    }
         //calculation of the normal to the boundary
         std::array<double, 3> n = {0.0, 0.0, 0.0}; 
         n[m_idim]=1-2*m_iside;
         dst.m_runtime_rdata[m_normal_index][dst_i]= n[0];
         dst.m_runtime_rdata[m_normal_index+1][dst_i]= n[1];
         dst.m_runtime_rdata[m_normal_index+2][dst_i]= n[2];
-
-#if (defined WARPX_DIM_3D)
-        [dst_i]-=(m_dt-delta_t)*ux;
-        dst.m_rdata[PIdx::y][dst_i]-=(m_dt-dedst.m_rdata[PIdx::x]lta_t)*uy;
-        dst.m_rdata[PIdx::z][dst_i]-=(m_dt-delta_t)*uz;
-#elif (defined WARPX_DIM_XZ)
-        dst.m_rdata[PIdx::x][dst_i]-=(m_dt-delta_t)*ux;
-        dst.m_rdata[PIdx::z][dst_i]-=(m_dt-delta_t)*uz;
-#elif (defined WARPX_DIM_RZ)
-        amrex::ParticleReal theta_temp=dst.m_rdata[PIdx::theta][dst_i];
-        amrex::ParticleReal x_temp=dst.m_rdata[PIdx::x][dst_i]*std::cos(theta_temp);
-        amrex::ParticleReal y_temp=dst.m_rdata[PIdx::y][dst_i]*std::sin(theta);
-        amrex::ParticleReal z_temp=dst.m_rdata[PIdx::z][dst_i];
-        x_temp -=(m_dt-delta_t)*ux;
-        y_temp -=(m_dt-delta_t)*uy;
-        z_temp -=(m_dt-delta_t)*uz;
-        theta_temp= std::atan2(y_temp, x_temp);
-        dst.m_rdata[PIdx::theta][dst_i] = theta_temp;
-        dst.m_rdata[PIdx::x][dst_i] =x_temp/std::cos(theta_temp);
-        dst.m_rdata[PIdx::z][dst_i] =z_temp/std::sin(theta_temp);
-#else (defined WARPX_DIM_1D_Z)
-        dst.m_rdata[PIdx::z][dst_i] -=(m_dt-delta_t)*uz;
-#endif
     }
 };
 
@@ -467,17 +421,15 @@ void ParticleBoundaryBuffer::gatherParticles (MultiParticleContainer& mypc,
                         {
                           WARPX_PROFILE("ParticleBoundaryBuffer::gatherParticles::filterAndTransform");
                           auto& warpx = WarpX::GetInstance();
-                          const auto dt = warpx.getdt(pti.GetLevel());
                           auto string_to_index_intcomp = buffer[i].getParticleRuntimeiComps();
                           const int step_scraped_index = string_to_index_intcomp.at("stepScraped");
                           auto string_to_index_realcomp = buffer[i].getParticleRuntimeComps();
                           const int delta_index = string_to_index_realcomp.at("deltaTimeScraped");
                           const int normal_index = string_to_index_realcomp.at("nx");
                           const int step = warpx_instance.getistep(0);
-
                           amrex::filterAndTransformParticles(ptile_buffer, ptile,
                                                              predicate,
-                                                             CopyAndTimestamp{step_scraped_index, delta_index, normal_index, step, dt, plo, phi, idim, iside},
+                                                             CopyAndTimestamp{step_scraped_index, delta_index, normal_index, step, idim, iside},
                                                              0, dst_index);
                         }
                     }
@@ -612,4 +564,4 @@ ParticleBoundaryBuffer::getParticleBufferPointer(const std::string species_name,
     auto index = WarpX::GetInstance().GetPartContainer().getSpeciesID(species_name);
 
     return &buffer[index];
-}
+}