dbsrsm_c.c
来自「用于求解大型稀疏线性方程组Ax=b的数值计算库.」· C语言 代码 · 共 1,474 行 · 第 1/5 页
C
1,474 行
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_MatTriangSlvLU_CaATBbC_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvLU_CaDATBbC_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvLU_CaATDBbC_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; case 2: /* Upper triangular, Unit diagonal */ switch (n) { case 1: /* Vec Mult */ if (alpha == 1) { if (beta == 1) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CATBC_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CDATBC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_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_VecTriangSlvUU_CATBmC_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CDATBmC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CATDBmC_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 */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CATB_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CDATB_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CATDB_double(mb, dv, 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_CATBbC_double(mb, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CDATBbC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CATDBbC_double(mb, dv, 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 */ } } else { /* alpha is general nonzero */ if (beta == 1) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CaATBC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDATBC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaATDBC_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 */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSR_VecTriangSlvUU_CaATBmC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDATBmC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaATDBmC_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_CaATB_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDATB_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaATDB_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_CaATBbC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ BSR_VecTriangSlvUU_CaDATBbC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ BSR_VecTriangSlvUU_CaATDBbC_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_CATBC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDATBC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CATDBC_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_CATBmC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDATBmC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CATDBmC_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_CATB_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDATB_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CATDB_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_CATBbC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSR_MatTriangSlvUU_CDATBbC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSR_MatTriangSlvUU_CATDBbC_double(mb, n, dv, val, bindx, ⌨️ 快捷键说明
复制代码Ctrl + C
搜索代码Ctrl + F
全屏模式F11
增大字号Ctrl + =
减小字号Ctrl + -
显示快捷键?