Commit dba63f25 authored by tilman.metz's avatar tilman.metz

bugfixes

parent 22e19348
......@@ -989,8 +989,8 @@ MPI_Barrier(MPI_COMM_WORLD);
}
} /*hloop*/
if(iteration==1 && HESS && !READ_HESS){
exchange_Fv(fourier_vel_back->Fvx_re,fourier_vel_back->Fvy_re,fourier_vel_back->Fvz_re,nf, velbuff, req_send, req_rec,0);
exchange_Fv(fourier_vel_back->Fvx_im,fourier_vel_back->Fvy_im,fourier_vel_back->Fvz_im,nf, velbuff, req_send, req_rec,0);
exchange_Fv(fourier_vel_back->Fvx_re,fourier_vel_back->Fvy_re,fourier_vel_back->Fvz_re,nf, velbuff, req_send, req_rec,ntr_hess);
exchange_Fv(fourier_vel_back->Fvx_im,fourier_vel_back->Fvy_im,fourier_vel_back->Fvz_im,nf, velbuff, req_send, req_rec,ntr_hess);
MPI_Barrier(MPI_COMM_WORLD);
hess_F(NX,NY,NZ,fourier_vel_fw,fourier_vel_back, hessian,nt,mod, finv, nf,ntr_hess);
......@@ -1046,7 +1046,7 @@ MPI_Barrier(MPI_COMM_WORLD);
modelupdate(NX,NY,NZ,grad,testmod,bfgsmod1, step[steptest],beta,it_group);
matcopy(testmod);
if (FREE_SURF)constant_boundary(xb[0],xb[1],yb[0],yb[1],zb[0],zb[1],testmod);
av_mat(mod,mod_av);
av_mat(testmod,mod_av);
L2=0.0;
L2f=0.0;
......
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