This commit is contained in:
张壹 2024-12-19 12:18:23 +08:00
parent 0988ff4a7f
commit 974f5cdd53

View File

@ -46,7 +46,7 @@ void gctl::gobser_wrt_mass_depth(double &depth_der, double &mass_der, const poin
int i; int i;
int o_size = obsp.size(); int o_size = obsp.size();
#pragma omp parallel for private (i, tmp_val, tmp_grad, tmp_obs, depth_mass, ders) schedule(guided) //#pragma omp parallel for private (i, tmp_val, tmp_grad, tmp_obs, depth_mass, ders) schedule(guided)
for (i = 0; i < o_size; i++) for (i = 0; i < o_size; i++)
{ {
depth_mass.resize(2); depth_mass.resize(2);
@ -59,12 +59,12 @@ void gctl::gobser_wrt_mass_depth(double &depth_der, double &mass_der, const poin
tmp_val = autodiff::val(tmp_obs); tmp_val = autodiff::val(tmp_obs);
tmp_grad = lf.gradient(tmp_val, i); tmp_grad = lf.gradient(tmp_val, i);
#pragma omp critical //#pragma omp critical
{ //{
lf.evaluate(tmp_val, i); lf.evaluate(tmp_val, i);
depth_der += tmp_grad*ders[0]; depth_der += tmp_grad*ders[0];
mass_der += tmp_grad*ders[1]; mass_der += tmp_grad*ders[1];
} //}
} }
return; return;