| 360 | |
| 361 | |
| 362 | void matrixmatrixmultiply(const ap::real_2d_array& a, |
| 363 | int ai1, |
| 364 | int ai2, |
| 365 | int aj1, |
| 366 | int aj2, |
| 367 | bool transa, |
| 368 | const ap::real_2d_array& b, |
| 369 | int bi1, |
| 370 | int bi2, |
| 371 | int bj1, |
| 372 | int bj2, |
| 373 | bool transb, |
| 374 | double alpha, |
| 375 | ap::real_2d_array& c, |
| 376 | int ci1, |
| 377 | int ci2, |
| 378 | int cj1, |
| 379 | int cj2, |
| 380 | double beta, |
| 381 | ap::real_1d_array& work) |
| 382 | { |
| 383 | int arows; |
| 384 | int acols; |
| 385 | int brows; |
| 386 | int bcols; |
| 387 | int crows; |
| 388 | int ccols; |
| 389 | int i; |
| 390 | int j; |
| 391 | int k = 0; |
| 392 | int l; |
| 393 | int r; |
| 394 | double v; |
| 395 | |
| 396 | |
| 397 | // |
| 398 | // Setup |
| 399 | // |
| 400 | if( !transa ) |
| 401 | { |
| 402 | arows = ai2-ai1+1; |
| 403 | acols = aj2-aj1+1; |
| 404 | } |
| 405 | else |
| 406 | { |
| 407 | arows = aj2-aj1+1; |
| 408 | acols = ai2-ai1+1; |
| 409 | } |
| 410 | if( !transb ) |
| 411 | { |
| 412 | brows = bi2-bi1+1; |
| 413 | bcols = bj2-bj1+1; |
| 414 | } |
| 415 | else |
| 416 | { |
| 417 | brows = bj2-bj1+1; |
| 418 | bcols = bi2-bi1+1; |
| 419 | } |
no test coverage detected