Level 2 variant of CMatrixRightTRSM *************************************************************************/
| 1674 | Level 2 variant of CMatrixRightTRSM |
| 1675 | *************************************************************************/ |
| 1676 | static void cmatrixrighttrsm2(int m, |
| 1677 | int n, |
| 1678 | const ap::complex_2d_array& a, |
| 1679 | int i1, |
| 1680 | int j1, |
| 1681 | bool isupper, |
| 1682 | bool isunit, |
| 1683 | int optype, |
| 1684 | ap::complex_2d_array& x, |
| 1685 | int i2, |
| 1686 | int j2) |
| 1687 | { |
| 1688 | int i; |
| 1689 | int j; |
| 1690 | ap::complex vc; |
| 1691 | ap::complex vd; |
| 1692 | |
| 1693 | |
| 1694 | // |
| 1695 | // Special case |
| 1696 | // |
| 1697 | if( n*m==0 ) |
| 1698 | { |
| 1699 | return; |
| 1700 | } |
| 1701 | |
| 1702 | // |
| 1703 | // Try to call fast TRSM |
| 1704 | // |
| 1705 | if( cmatrixrighttrsmf(m, n, a, i1, j1, isupper, isunit, optype, x, i2, j2) ) |
| 1706 | { |
| 1707 | return; |
| 1708 | } |
| 1709 | |
| 1710 | // |
| 1711 | // General case |
| 1712 | // |
| 1713 | if( isupper ) |
| 1714 | { |
| 1715 | |
| 1716 | // |
| 1717 | // Upper triangular matrix |
| 1718 | // |
| 1719 | if( optype==0 ) |
| 1720 | { |
| 1721 | |
| 1722 | // |
| 1723 | // X*A^(-1) |
| 1724 | // |
| 1725 | for(i = 0; i <= m-1; i++) |
| 1726 | { |
| 1727 | for(j = 0; j <= n-1; j++) |
| 1728 | { |
| 1729 | if( isunit ) |
| 1730 | { |
| 1731 | vd = 1; |
| 1732 | } |
| 1733 | else |
no test coverage detected