Commit 960d9115 authored by Michel Henry's avatar Michel Henry
Browse files

wip

parent ca37deae
......@@ -122,6 +122,8 @@ void fluid_problem_node_force_volume(FluidProblem *problem, const double *soluti
}
double f[D],dfdu,dfddp;
particle_force_f(problem,f,&dfdu,&dfddp,s,ds,sold,c,dc,a,dt,iel,ip);
if(ip == 0) printf("Force : %14.7f\n", f[0]);
if (!local_vector)
continue;
double rho,nu;
......
......@@ -120,7 +120,6 @@ tic = time.time()
while t < tEnd :
time_integration.iterate(fluid, p, dt, min_nsub=5, external_particles_forces=g*p.mass())
t += dt
# Output files writting
if ii %outf == 0 :
ioutput = int(ii/outf) + 1
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment