@@ -156,7 +156,7 @@ double *LagrangeOptimizer::computeLagrangeParameters(const double *reconData) {
156156 for (int i = 0 ; i < myPlaneCount; i++) {
157157#pragma omp parallel for default(none) \
158158 shared (i, myPlaneCount, myNodeCount, myVxCount, myVyCount, reconData, \
159- i_g) private (lindex, rindex)
159+ i_g) private (lindex, rindex)
160160 for (int k = 0 ; k < myNodeCount; k++) {
161161 for (int j = 0 ; j < myVxCount; j++) {
162162 for (int l = 0 ; l < myVyCount; l++) {
@@ -180,7 +180,7 @@ double *LagrangeOptimizer::computeLagrangeParameters(const double *reconData) {
180180 std::vector<double > V4 (myNodeCount * myVxCount * myVyCount, 0 );
181181#pragma omp parallel for default(none) \
182182 shared (myNodeCount, myVxCount, myVyCount, myVolume, myVth, myVp, myMuQoi, \
183- myVth2, myParticleMass, V2, V3, V4) private (i, j, l, m)
183+ myVth2, myParticleMass, V2, V3, V4) private (i, j, l, m)
184184 for (k = 0 ; k < myNodeCount * myVxCount * myVyCount; ++k) {
185185 i = int (k / (myVxCount * myVyCount));
186186 j = int (k % myVyCount);
@@ -210,8 +210,8 @@ double *LagrangeOptimizer::computeLagrangeParameters(const double *reconData) {
210210 std::vector<double > Tperp (myNodeCount, 0 );
211211#pragma omp parallel for default(none) \
212212 shared (myNodeCount, myVxCount, myVyCount, myVolume, myVth, myVp, f0_f, \
213- myMuQoi, myVth2, myParticleMass, mySmallElectronCharge, D, U, \
214- Tperp) private (i, j, l, m)
213+ myMuQoi, myVth2, myParticleMass, mySmallElectronCharge, D, U, \
214+ Tperp) private (i, j, l, m)
215215 for (k = 0 ; k < myNodeCount * myVxCount * myVyCount; ++k) {
216216 i = int (k / (myVxCount * myVyCount));
217217 j = int (k % myVyCount);
@@ -225,9 +225,10 @@ double *LagrangeOptimizer::computeLagrangeParameters(const double *reconData) {
225225 std::vector<double > Tpara (myNodeCount, 0 );
226226 std::vector<double > Rpara (myNodeCount, 0 );
227227 double en;
228- #pragma omp parallel for default(none) shared( \
229- myNodeCount, myVxCount, myVyCount, myVolume, myVth, myVp, f0_f, myVth2, \
230- myParticleMass, mySmallElectronCharge, D, U, Tpara) private (i, j, en)
228+ #pragma omp parallel for default(none) \
229+ shared (myNodeCount, myVxCount, myVyCount, myVolume, myVth, myVp, f0_f, \
230+ myVth2, myParticleMass, mySmallElectronCharge, D, U, \
231+ Tpara) private (i, j, en)
231232 for (k = 0 ; k < myNodeCount * myVxCount * myVyCount; ++k) {
232233 i = int (k / (myVxCount * myVyCount));
233234 j = int (k % myVyCount);
@@ -236,9 +237,9 @@ double *LagrangeOptimizer::computeLagrangeParameters(const double *reconData) {
236237 (f0_f[k] * myVolume[k] * en * myVth2[i] * myParticleMass) /
237238 D[i] / mySmallElectronCharge;
238239 }
239- #pragma omp parallel for default(none) \
240- shared ( myNodeCount, myVxCount, myVyCount, myVolume, myVth, myVth2, \
241- myParticleMass, mySmallElectronCharge, U, Tpara, Rpara) private (i)
240+ #pragma omp parallel for default(none) shared( \
241+ myNodeCount, myVxCount, myVyCount, myVolume, myVth, myVth2, \
242+ myParticleMass, mySmallElectronCharge, U, Tpara, Rpara) private (i)
242243 for (k = 0 ; k < myNodeCount * myVxCount * myVyCount; ++k) {
243244 i = int (k / (myVxCount * myVyCount));
244245 Rpara[i] = mySmallElectronCharge * Tpara[i] +
@@ -282,8 +283,8 @@ double *LagrangeOptimizer::computeLagrangeParameters(const double *reconData) {
282283 int maxIter = 50 ;
283284#pragma omp parallel for default(none) \
284285 shared (reconData, iphi, D, U, V2, V3, V4, f0_f, Tperp, Rpara, DeB, UeB, \
285- TperpEB, TparaEB, PDeB, maxIter, node_unconv, \
286- my_rank) private (count_unLag, breg_recon)
286+ TperpEB, TparaEB, PDeB, maxIter, node_unconv, \
287+ my_rank) private (count_unLag, breg_recon)
287288 for (idx = 0 ; idx < myNodeCount; ++idx) {
288289 int count = 0 ;
289290 double gradients[4 ] = {0.0 , 0.0 , 0.0 , 0.0 };
0 commit comments