📄 dbsrsm_c.c
字号:
break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CaABmC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDABmC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaADBmC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CaAB_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDAB_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaADB_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CaABbC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDABbC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaADBbC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } } break; default: /* Mat Mult */ if (alpha == 1) { if (beta == 1) { switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CABC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDABC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CADBC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CABmC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDABmC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CADBmC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CAB_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDAB_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CADB_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CABbC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDABbC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CADBbC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } } else { /* alpha is general nonzero */ if (beta == 1) { switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CaABC_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CaDABC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CaADBC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CaABmC_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CaDABmC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CaADBmC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CaAB_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CaDAB_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CaADB_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ BSR_MatTriangSlvUU_CaABbC_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CaDABbC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CaADBbC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } } break; } /* end switch on n */ break; default: /* Invalid descra[1] */ printf("Check value of descra[1] (descra(2) in Fortran.\n"); printf("Valid values of descra[1] (descra(2)) are 1 and 2.\n"); break; } /* end switch on descra[1] */ break; case 1: switch ( descra[1] ) { case 1: /* Lower triangular, Unit diagonal */ switch (n) { case 1: /* Vec Mult */ if (alpha == 1) { if (beta == 1) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvLU_CATBC_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvLU_CDATBC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvLU_CATDBC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dbsrsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvLU_CATBmC_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvLU_CDATBmC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, ⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -