📄 dbscsm_c.c
字号:
printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSC_MatTriangSlvLU_CaATB_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; case 2: /* Left scaling */ BSC_MatTriangSlvLU_CaDATB_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSC_MatTriangSlvLU_CaATDB_double(mb, n, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ BSC_MatTriangSlvLU_CaATBbC_double(mb, n, alpha, val, bindx, bpntrb, bpntre, lb, b, ldb, beta, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSC_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 */ BSC_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 dbscsm. 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 */ BSC_VecTriangSlvUU_CATBC_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CDATBC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CATDBC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CATBmC_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CDATBmC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CATDBmC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CATB_double(mb, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CDATB_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CATDB_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CATBbC_double(mb, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CDATBbC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CATDBbC_double(mb, dv, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */ } } else { /* alpha is general nonzero */ if (beta == 1) { switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CaATBC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CaDATBC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CaATDBC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CaATBmC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CaDATBmC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CaATDBmC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CaATB_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CaDATB_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CaATDB_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, c, ind_base); break; default: printf("Invalid argument unitd in dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ BSC_VecTriangSlvUU_CaATBbC_double(mb, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ BSC_VecTriangSlvUU_CaDATBbC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ BSC_VecTriangSlvUU_CaATDBbC_double(mb, dv, alpha, val, bindx, bpntrb, bpntre, lb, b, beta, c, work, ind_base); break; default: printf("Invalid argument unitd in dbscsm. 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 */ BSC_MatTriangSlvUU_CATBC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSC_MatTriangSlvUU_CDATBC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSC_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 dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ BSC_MatTriangSlvUU_CATBmC_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ BSC_MatTriangSlvUU_CDATBmC_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSC_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 dbscsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ BSC_MatTriangSlvUU_CATB_double(mb, n, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, ind_base); break; case 2: /* Left scaling */ BSC_MatTriangSlvUU_CDATB_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ BSC_MatTriangSlvUU_CATDB_double(mb, n, dv, val, bindx, bpntrb, bpntre, lb, b, ldb, c, l⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -