@@ -188,8 +188,8 @@ scrapeParticles (PC& pc, const amrex::Vector<const amrex::MultiFab*>& distance_t
188
188
189
189
amrex::ParticleReal xp, yp, zp;
190
190
getPosition (ip, xp, yp, zp);
191
-
192
-
191
+
192
+
193
193
194
194
int i, j, k;
195
195
amrex::Real W[AMREX_SPACEDIM][2 ];
@@ -200,12 +200,13 @@ scrapeParticles (PC& pc, const amrex::Vector<const amrex::MultiFab*>& distance_t
200
200
if (phi_value < 0.0 )
201
201
{
202
202
amrex::Real x_temp, y_temp, z_temp;
203
- amrex::Real dt_fraction = bisect ( 0 , 1 ,
203
+ amrex::Real dt_fraction = amrex:: bisect ( 0.0 , 1.0 ,
204
204
[=] AMREX_GPU_DEVICE (amrex::Real dt_frac) mutable {
205
205
x_temp = xp; y_temp = yp; z_temp = zp;
206
206
UpdatePosition (x_temp, y_temp, z_temp, ux[ip], uy[ip], uz[ip], -dt_frac*dt);
207
207
ablastr::particles::compute_weights_nodal (x_temp, y_temp, z_temp, plo, dxi, i, j, k, W);
208
208
phi_value = ablastr::particles::interp_field_nodal (i, j, k, W, phi);
209
+ return phi_value;
209
210
} );
210
211
211
212
amrex::RealVect normal = DistanceToEB::interp_normal (i, j, k, W, phi, dxi);
0 commit comments