Skip to content

GitLab

Commit 228e0d5a authored by Kirill Terekhov's avatar Kirill Terekhov
Browse files

Fix warnings from intel compiler

parent dedca4f9
Hide whitespace changes
Inline Side-by-side
Source/Solver/solver_inner/solver_mlmptiluc/solver_mlmtiluc2.cpp
View file @ 228e0d5a
......@@ -1848,9 +1848,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
//T = Timer();
#if defined(USE_OMP_FACT)
#pragma omp parallel for private(l,u,i,j)
#pragma omp parallel for private(k,l,u,i,j)
#endif
for (int k = cbeg; k < cend; ++k)
for (k = cbeg; k < cend; ++k)
{
if( V[k] == std::numeric_limits<INMOST_DATA_REAL_TYPE>::max() ) l = 1;
else l = exp(V[k]);
......@@ -2318,9 +2318,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
#if defined(USE_OMP_FACT)
#pragma omp parallel for
#pragma omp parallel for private(k)
#endif
for(int k = wbeg; k < wend; ++k)
for(k = wbeg; k < wend; ++k)
{
for (INMOST_DATA_ENUM_TYPE jt = A_Address[k].first; jt < A_Address[k].last; ++jt)
{
......@@ -3037,9 +3037,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
B_Entries.push_back(Sparse::Row::make_entry(-1,0.0));
#if defined(USE_OMP_FACT)
#pragma omp parallel for
#pragma omp parallel for private(k)
#endif
for (int k = cbeg; k < cend; ++k)
for (k = cbeg; k < cend; ++k)
std::sort(B_Entries.begin() + B_Address[k].first, B_Entries.begin() + B_Address[k].last);
double sparsity = nzA/(double)(cend-cbeg)/(double)(cend-cbeg);
......@@ -3131,9 +3131,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = cbeg; k < cend; k++)
for (k = cbeg; k < cend; k++)
{
for (INMOST_DATA_ENUM_TYPE r = B_Address[k].first; r < B_Address[k].last; ++r)
B_Entries[r].second *= (DL[k] * DR[B_Entries[r].first]);
......@@ -3338,9 +3338,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
}
#if defined(USE_OMP_FACT)
#pragma omp for private(u)
#pragma omp for private(k,u)
#endif
for(int k = cbeg; k < cend; ++k)
for(k = cbeg; k < cend; ++k)
{
for (INMOST_DATA_ENUM_TYPE r = B_Address[k].first; r < B_Address[k].last; ++r)
{
......@@ -3355,18 +3355,18 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
//std::fill(temp.begin() + cbeg - mobeg, temp.begin() + cend - mobeg, 0.0);
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for(int k = cbeg; k < cend; ++k) temp[k] = 0.0;
for(k = cbeg; k < cend; ++k) temp[k] = 0.0;
for (INMOST_DATA_ENUM_TYPE iter = 0; iter < sciters; iter++)
{
//std::fill(DL.Begin() + cbeg - mobeg, DL.Begin() + cend - mobeg, std::numeric_limits<INMOST_DATA_REAL_TYPE>::max());
//std::fill(DR.Begin() + cbeg - mobeg, DR.Begin() + cend - mobeg, std::numeric_limits<INMOST_DATA_REAL_TYPE>::max());
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for(int k = cbeg; k < cend; ++k) DL[k] = DR[k] = std::numeric_limits<INMOST_DATA_REAL_TYPE>::max();
for(k = cbeg; k < cend; ++k) DL[k] = DR[k] = std::numeric_limits<INMOST_DATA_REAL_TYPE>::max();
#if defined(USE_OMP_FACT)
//#pragma omp for private(k,i,u)
......@@ -3388,9 +3388,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
}
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = cbeg; k < cend; k++)
for (k = cbeg; k < cend; k++)
{
if( DR[k] != std::numeric_limits<INMOST_DATA_REAL_TYPE>::max() &&
DL[k] != std::numeric_limits<INMOST_DATA_REAL_TYPE>::max() )
......@@ -3400,9 +3400,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
}
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = cbeg; k < cend; k++)
for (k = cbeg; k < cend; k++)
{
DL[k] = exp(-temp[k]);
DR[k] = exp(temp[k]);
......@@ -3412,9 +3412,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
//if( isnan(DR[k]) || DR[k] != DR[k] || fabs(DR[k]) < 1.0e-12 ) std::cout << __FILE__ << ":" << __LINE__ << " DR[" << k << "] is " << DR[k] << std::endl;
}
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = cbeg; k < cend; k++)
for (k = cbeg; k < cend; k++)
{
for (INMOST_DATA_ENUM_TYPE r = B_Address[k].first; r < B_Address[k].last; ++r)
B_Entries[r].second *= (DL[k] * DR[B_Entries[r].first]);
......@@ -3442,9 +3442,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
}
*/
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = cbeg; k < cend; k++)
for (k = cbeg; k < cend; k++)
{
DL[k] *= U[k];
DR[k] *= V[k];
......@@ -3456,9 +3456,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
//stack scaling
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = cbeg; k < cend; k++)
for (k = cbeg; k < cend; k++)
{
DL0[k] *= DL[k];
DR0[k] *= DR[k];
......@@ -3489,9 +3489,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
//std::cout << "rescaling [" << cbeg << "," << cend << "] ";
//std::cout << "selected [" << kbeg << ":" << kend << "]" << std::endl;
#if defined(USE_OMP_FACT)
#pragma omp for
#pragma omp for private(k)
#endif
for (int k = kbeg; k < kend; ++k)
for (k = kbeg; k < kend; ++k)
{
if( F_Address[level]->at(k).Size() )
{
......@@ -4558,9 +4558,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
printf("Assemble LF\n");
#if defined(USE_OMP_FACT)
#pragma omp parallel for private(Li, u, j, curr, next, LFmin, LFmax, LFnorm, LFnum, LFtau, LFdrop) reduction(+:ndrops_lf)
#pragma omp parallel for private(k, Li, u, j, curr, next, LFmin, LFmax, LFnorm, LFnum, LFtau, LFdrop) reduction(+:ndrops_lf)
#endif
for (int k = cend; k < wend; ++k)
for (k = cend; k < wend; ++k)
{
int nthr = Thread();
interval<INMOST_DATA_ENUM_TYPE,INMOST_DATA_REAL_TYPE> & LineValues = LineValuesS[nthr];
......@@ -4842,9 +4842,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
if( verbosity > 1 )
printf("Assemble EU\n");
#if defined(USE_OMP_FACT)
#pragma omp parallel for private(Li, u, i, j, curr, next, EUmin, EUmax, EUnorm, EUnum, EUtau, EUdrop) reduction(+:ndrops_eu)
#pragma omp parallel for private(k, Li, u, i, j, curr, next, EUmin, EUmax, EUnorm, EUnum, EUtau, EUdrop) reduction(+:ndrops_eu)
#endif
for(int k = cend; k < wend; ++k)
for(k = cend; k < wend; ++k)
{
int nthr = Thread();
interval<INMOST_DATA_ENUM_TYPE,INMOST_DATA_REAL_TYPE> & LineValues = LineValuesS[nthr];
......@@ -5063,7 +5063,7 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
#if defined(USE_OMP_FACT)
#pragma omp for
#endif
for(int k = cend; k < wend; ++k)
for(k = cend; k < wend; ++k)
{
INMOST_DATA_REAL_TYPE v, coef;
INMOST_DATA_ENUM_TYPE Sbeg = EOL, i, j;
......@@ -5102,7 +5102,7 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
#if defined(USE_OMP_FACT)
#pragma omp for
#endif
for(int k = cend; k < wend; ++k)
for(k = cend; k < wend; ++k)
{
for(unsigned j = 0; j < Scolmax_local.size(); ++j)
Scolmax[k] = std::max(Scolmax[k],Scolmax_local[j][k]);
......@@ -5124,9 +5124,9 @@ const INMOST_DATA_ENUM_TYPE UNDEF = ENUMUNDEF, EOL = ENUMUNDEF - 1;
nnzA = 0;
for(k = cend; k < wend; ++k) nnzA += A_Address[k].Size();
#if defined(USE_OMP_FACT)
#pragma omp parallel for private(Sbeg, Li, Ui, i,j, u,l, Smin,Smax,Snorm,Snum,Stau,Sdrop) reduction(+:ndrops_s)
#pragma omp parallel for private(k,Sbeg, Li, Ui, i,j, u,l, Smin,Smax,Snorm,Snum,Stau,Sdrop) reduction(+:ndrops_s)
#endif
for(int k = cend; k < wend; ++k)
for(k = cend; k < wend; ++k)
{
int nthr = Thread();
interval<INMOST_DATA_ENUM_TYPE,INMOST_DATA_REAL_TYPE> & LineValues = LineValuesS[nthr];
......
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