@@ -173,46 +173,46 @@ int main(int argc, char *argv[]){
173
173
#ifndef COMPLEX
174
174
if (uplos & 1 ) {
175
175
for (j = 0 ; j < m ; j ++ ) {
176
- for (i = 0 ; i < j ; i ++ ) a [i + j * m ] = 0. ;
177
- a [j + j * m ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
178
- for (i = j + 1 ; i < m ; i ++ ) a [i + j * m ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
176
+ for (i = 0 ; i < j ; i ++ ) a [( long ) i + ( long ) j * ( long ) m ] = 0. ;
177
+ a [( long ) j + ( long ) j * ( long ) m ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
178
+ for (i = j + 1 ; i < m ; i ++ ) a [( long ) i + ( long ) j * ( long ) m ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
179
179
}
180
180
} else {
181
181
for (j = 0 ; j < m ; j ++ ) {
182
- for (i = 0 ; i < j ; i ++ ) a [i + j * m ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
183
- a [j + j * m ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
184
- for (i = j + 1 ; i < m ; i ++ ) a [i + j * m ] = 0. ;
182
+ for (i = 0 ; i < j ; i ++ ) a [( long ) i + ( long ) j * ( long ) m ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
183
+ a [( long ) j + ( long ) j * ( long ) m ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
184
+ for (i = j + 1 ; i < m ; i ++ ) a [( long ) i + ( long ) j * ( long ) m ] = 0. ;
185
185
}
186
186
}
187
187
#else
188
188
if (uplos & 1 ) {
189
189
for (j = 0 ; j < m ; j ++ ) {
190
190
for (i = 0 ; i < j ; i ++ ) {
191
- a [(i + j * m ) * 2 + 0 ] = 0. ;
192
- a [(i + j * m ) * 2 + 1 ] = 0. ;
191
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 0 ] = 0. ;
192
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 1 ] = 0. ;
193
193
}
194
194
195
- a [(j + j * m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
196
- a [(j + j * m ) * 2 + 1 ] = 0. ;
195
+ a [(( long ) j + ( long ) j * ( long ) m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
196
+ a [(( long ) j + ( long ) j * ( long ) m ) * 2 + 1 ] = 0. ;
197
197
198
198
for (i = j + 1 ; i < m ; i ++ ) {
199
- a [(i + j * m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
200
- a [(i + j * m ) * 2 + 1 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
199
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
200
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 1 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
201
201
}
202
202
}
203
203
} else {
204
204
for (j = 0 ; j < m ; j ++ ) {
205
205
for (i = 0 ; i < j ; i ++ ) {
206
- a [(i + j * m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
207
- a [(i + j * m ) * 2 + 1 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
206
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
207
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 1 ] = ((double ) rand () / (double ) RAND_MAX ) - 0.5 ;
208
208
}
209
209
210
- a [(j + j * m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
211
- a [(j + j * m ) * 2 + 1 ] = 0. ;
210
+ a [(( long ) j + ( long ) j * ( long ) m ) * 2 + 0 ] = ((double ) rand () / (double ) RAND_MAX ) + 8. ;
211
+ a [(( long ) j + ( long ) j * ( long ) m ) * 2 + 1 ] = 0. ;
212
212
213
213
for (i = j + 1 ; i < m ; i ++ ) {
214
- a [(i + j * m ) * 2 + 0 ] = 0. ;
215
- a [(i + j * m ) * 2 + 1 ] = 0. ;
214
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 0 ] = 0. ;
215
+ a [(( long ) i + ( long ) j * ( long ) m ) * 2 + 1 ] = 0. ;
216
216
}
217
217
}
218
218
}
@@ -239,21 +239,27 @@ int main(int argc, char *argv[]){
239
239
for (j = 0 ; j < m ; j ++ ) {
240
240
for (i = 0 ; i <= j ; i ++ ) {
241
241
#ifndef COMPLEX
242
- if (maxerr < fabs (a [i + j * m ] - b [i + j * m ])) maxerr = fabs (a [i + j * m ] - b [i + j * m ]);
242
+ if (maxerr < fabs (a [(long )i + (long )j * (long )m ] - b [(long )i + (long )j * (long )m ]))
243
+ maxerr = fabs (a [(long )i + (long )j * (long )m ] - b [(long )i + (long )j * (long )m ]);
243
244
#else
244
- if (maxerr < fabs (a [(i + j * m ) * 2 + 0 ] - b [(i + j * m ) * 2 + 0 ])) maxerr = fabs (a [(i + j * m ) * 2 + 0 ] - b [(i + j * m ) * 2 + 0 ]);
245
- if (maxerr < fabs (a [(i + j * m ) * 2 + 1 ] - b [(i + j * m ) * 2 + 1 ])) maxerr = fabs (a [(i + j * m ) * 2 + 1 ] - b [(i + j * m ) * 2 + 1 ]);
245
+ if (maxerr < fabs (a [((long )i + (long )j * (long )m ) * 2 + 0 ] - b [((long )i + (long )j * (long )m ) * 2 + 0 ]))
246
+ maxerr = fabs (a [((long )i + (long )j * (long )m ) * 2 + 0 ] - b [((long )i + (long )j * (long )m ) * 2 + 0 ]);
247
+ if (maxerr < fabs (a [((long )i + (long )j * (long )m ) * 2 + 1 ] - b [((long )i + (long )j * (long )m ) * 2 + 1 ]))
248
+ maxerr = fabs (a [((long )i + (long )j * (long )m ) * 2 + 1 ] - b [((long )i + (long )j * (long )m ) * 2 + 1 ]);
246
249
#endif
247
250
}
248
251
}
249
252
} else {
250
253
for (j = 0 ; j < m ; j ++ ) {
251
254
for (i = j ; i < m ; i ++ ) {
252
255
#ifndef COMPLEX
253
- if (maxerr < fabs (a [i + j * m ] - b [i + j * m ])) maxerr = fabs (a [i + j * m ] - b [i + j * m ]);
256
+ if (maxerr < fabs (a [(long )i + (long )j * (long )m ] - b [(long )i + (long )j * (long )m ]))
257
+ maxerr = fabs (a [(long )i + (long )j * (long )m ] - b [(long )i + (long )j * (long )m ]);
254
258
#else
255
- if (maxerr < fabs (a [(i + j * m ) * 2 + 0 ] - b [(i + j * m ) * 2 + 0 ])) maxerr = fabs (a [(i + j * m ) * 2 + 0 ] - b [(i + j * m ) * 2 + 0 ]);
256
- if (maxerr < fabs (a [(i + j * m ) * 2 + 1 ] - b [(i + j * m ) * 2 + 1 ])) maxerr = fabs (a [(i + j * m ) * 2 + 1 ] - b [(i + j * m ) * 2 + 1 ]);
259
+ if (maxerr < fabs (a [((long )i + (long )j * (long )m ) * 2 + 0 ] - b [((long )i + (long )j * (long )m ) * 2 + 0 ]))
260
+ maxerr = fabs (a [((long )i + (long )j * (long )m ) * 2 + 0 ] - b [((long )i + (long )j * (long )m ) * 2 + 0 ]);
261
+ if (maxerr < fabs (a [((long )i + (long )j * (long )m ) * 2 + 1 ] - b [((long )i + (long )j * (long )m ) * 2 + 1 ]))
262
+ maxerr = fabs (a [((long )i + (long )j * (long )m ) * 2 + 1 ] - b [((long )i + (long )j * (long )m ) * 2 + 1 ]);
257
263
#endif
258
264
}
259
265
}
0 commit comments