40{
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229 char Dir, * one, * zero;
230 Int Afr, conjg, k, kb, kbb, kend, kstart, kstep, ktmp;
231
232
233
235 char * Aptr = NULL;
236
237
238
239
240
241
242
244
247
249 {
251 kstart = 0; kend = ( (
MIN( M, N ) - 1 ) / kb + 1 ) * kb; kstep = kb;
252 }
253 else
254 {
256 kstart = ( (
MIN( M, N ) - 1 ) / kb ) * kb; kend = kstep = -kb;
257 }
258
260 {
262 {
263 if( M >= N )
264 {
265 for( k = kstart; k != kend; k += kstep )
266 {
267 kbb = N - k; kbb =
MIN( kbb, kb ); ktmp = k + kbb;
268
269
270
272 COLUMN, &Aptr, DBUFA, &Afr );
273
274
275
277 DBUFA );
278
279
280
281 if( kbb > 1 )
283 Aptr, k+1, 0, DBUFA );
284
285
286
288 one, C, IC, JC+k, DESCC,
COLUMN );
289
290 if( Afr ) free( Aptr );
291 }
292 }
293 else
294 {
295 for( k = kstart; k != kend; k += kstep )
296 {
297 kbb = M - k; kbb =
MIN( kbb, kb ); ktmp = N - k;
298
299
300
302 DESCA,
ROW, &Aptr, DBUFA, &Afr );
303
304
305
307 DBUFA );
308
309
310
311 if( kbb > 1 )
313 Aptr, 1, 0, DBUFA );
314
315
316
318 one, C, IC+k, JC+k, DESCC,
ROW );
319
320 if( Afr ) free( Aptr );
321 }
322 }
323 }
324 else
325 {
326 if( M >= N )
327 {
328 for( k = kstart; k != kend; k += kstep )
329 {
330 kbb = N - k; kbb =
MIN( kbb, kb ); ktmp = M - k;
331
332
333
335 DESCA,
COLUMN, &Aptr, DBUFA, &Afr );
336
337
338
340 DBUFA );
341
342
343
344 if( kbb > 1 )
346 Aptr, 0, 1, DBUFA );
347
348
349
351 one, C, IC+k, JC+k, DESCC,
COLUMN );
352
353 if( Afr ) free( Aptr );
354 }
355 }
356 else
357 {
358 for( k = kstart; k != kend; k += kstep )
359 {
360 kbb = M - k; kbb =
MIN( kbb, kb ); ktmp = k + kbb;
361
362
363
365 DESCA,
ROW, &Aptr, DBUFA, &Afr );
366
367
368
370 DBUFA );
371
372
373
374 if( kbb > 1 )
376 Aptr, 0, k+1, DBUFA );
377
378
379
381 one, C, IC+k, JC, DESCC,
ROW );
382
383 if( Afr ) free( Aptr );
384 }
385 }
386 }
387 }
388 else
389 {
391
393 {
394 if( M >= N )
395 {
396 for( k = kstart; k != kend; k += kstep )
397 {
398 kbb = N - k; kbb =
MIN( kbb, kb ); ktmp = k + kbb;
399
400
401
403 ROW, &Aptr, DBUFA, &Afr );
404
405
406
407 if( conjg )
409 else
411 0, 0, DBUFA );
412
413
414
415 if( kbb > 1 )
417 Aptr, 0, k+1, DBUFA );
418
419
420
422 one, C, IC, JC+k, DESCC,
COLUMN );
423
424 if( Afr ) free( Aptr );
425 }
426 }
427 else
428 {
429 for( k = kstart; k != kend; k += kstep )
430 {
431 kbb = M - k; kbb =
MIN( kbb, kb ); ktmp = N - k;
432
433
434
436 DESCA,
COLUMN, &Aptr, DBUFA, &Afr );
437
438
439
440 if( conjg )
442 else
444 0, 0, DBUFA );
445
446
447
448 if( kbb > 1 )
450 Aptr, 0, 1, DBUFA );
451
452
453
455 one, C, IC+k, JC+k, DESCC,
ROW );
456
457 if( Afr ) free( Aptr );
458 }
459 }
460 }
461 else
462 {
463 if( M >= N )
464 {
465 for( k = kstart; k != kend; k += kstep )
466 {
467 kbb = N - k; kbb =
MIN( kbb, kb ); ktmp = M - k;
468
469
470
472 DESCA,
ROW, &Aptr, DBUFA, &Afr );
473
474
475
476 if( conjg )
478 else
480 0, 0, DBUFA );
481
482
483
484 if( kbb > 1 )
486 Aptr, 1, 0, DBUFA );
487
488
489
491 one, C, IC+k, JC+k, DESCC,
COLUMN );
492
493 if( Afr ) free( Aptr );
494 }
495 }
496 else
497 {
498 for( k = kstart; k != kend; k += kstep )
499 {
500 kbb = M - k; kbb =
MIN( kbb, kb ); ktmp = k + kbb;
501
502
503
505 DESCA,
COLUMN, &Aptr, DBUFA, &Afr );
506
507
508
509 if( conjg )
511 else
513 0, 0, DBUFA );
514
515
516
517 if( kbb > 1 )
519 Aptr, k+1, 0, DBUFA );
520
521
522
524 one, C, IC+k, JC, DESCC,
ROW );
525
526 if( Afr ) free( Aptr );
527 }
528 }
529 }
530 }
531
532
533
534}