@@ -40,9 +40,10 @@ void RemoveWeights(G3SkyMap &T, G3SkyMap &Q, G3SkyMap &U, const G3SkyMapWeights
4040 T.ConvertToDense ();
4141 Q.ConvertToDense ();
4242 U.ConvertToDense ();
43+
4344 #ifdef _OPENMP
44- #pragma omp parallel for schedule(static)
45- #endif
45+ #pragma omp parallel for
46+ #endif
4647 for (size_t pix = 0 ; pix < T.size (); pix++) {
4748 StokesVector v (T[pix], Q[pix], U[pix]);
4849 v /= W.at (pix);
@@ -127,9 +128,9 @@ void ApplyWeights(G3SkyMap &T, G3SkyMap &Q, G3SkyMap &U, const G3SkyMapWeights &
127128 g3_assert (!Q.weighted );
128129 g3_assert (!U.weighted );
129130
130- #ifdef _OPENMP
131- #pragma omp parallel for schedule(static)
132- #endif
131+ #ifdef _OPENMP
132+ #pragma omp parallel for
133+ #endif
133134 for (size_t pix = 0 ; pix < T.size (); pix++) {
134135 if (T.at (pix) == 0 && Q.at (pix) == 0 && U.at (pix) == 0 )
135136 continue ;
@@ -163,9 +164,9 @@ py::tuple GetRaDecMap(const G3SkyMap &m)
163164 ra->ConvertToDense ();
164165 dec->ConvertToDense ();
165166
166- #ifdef _OPENMP
167- #pragma omp parallel for schedule(static)
168- #endif
167+ #ifdef _OPENMP
168+ #pragma omp parallel for
169+ #endif
169170 for (size_t i = 0 ; i < m.size (); i++) {
170171 std::vector<double > radec = m.PixelToAngle (i);
171172 (*ra)[i] = radec[0 ];
@@ -200,9 +201,9 @@ G3SkyMapMaskPtr GetRaDecMask(const G3SkyMap &m, double ra_left, double ra_right,
200201 ra_left = wrap_ra (ra_left);
201202 ra_right = wrap_ra (ra_right);
202203
203- #ifdef _OPENMP
204- #pragma omp parallel for schedule(static)
205- #endif
204+ #ifdef _OPENMP
205+ #pragma omp parallel for
206+ #endif
206207 for (size_t i = 0 ; i < m.size (); i++) {
207208 std::vector<double > radec = m.PixelToAngle (i);
208209 double ra = wrap_ra (radec[0 ]);
@@ -228,9 +229,10 @@ G3SkyMapMaskPtr GetGalacticPlaneMask(const G3SkyMap &m, double lat)
228229
229230 if (m.coord_ref == MapCoordReference::Equatorial) {
230231 auto q_rot = get_fk5_j2000_to_gal_quat ();
231- #ifdef _OPENMP
232- #pragma omp parallel for schedule(static)
233- #endif
232+
233+ #ifdef _OPENMP
234+ #pragma omp parallel for
235+ #endif
234236 for (size_t i = 0 ; i < m.size (); i++) {
235237 // compute just latitude part of each coordinate
236238 auto q = q_rot * m.PixelToQuat (i);
@@ -241,8 +243,8 @@ G3SkyMapMaskPtr GetGalacticPlaneMask(const G3SkyMap &m, double lat)
241243 }
242244 } else if (m.coord_ref == MapCoordReference::Galactic) {
243245 #ifdef _OPENMP
244- #pragma omp parallel for schedule(static)
245- #endif
246+ #pragma omp parallel for
247+ #endif
246248 for (size_t i = 0 ; i < m.size (); i++) {
247249 auto q = m.PixelToQuat (i);
248250 if (fabs (q.d ()) <= slat)
0 commit comments