printf(" total j: %f %f %f\n",aj->f.x,aj->f.y,aj->f.z);
printf(" energy: %f = %f %f %f %f\n",0.5*f_c*(b*f_a+f_r),
f_c,b,f_a,f_r);
printf(" %f %f %f\n",exchange->zeta_ij,.0,.0);
}
printf(" total j: %f %f %f\n",aj->f.x,aj->f.y,aj->f.z);
printf(" energy: %f = %f %f %f %f\n",0.5*f_c*(b*f_a+f_r),
f_c,b,f_a,f_r);
printf(" %f %f %f\n",exchange->zeta_ij,.0,.0);
}