📄 dvbrsm_c.c
字号:
} else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ VBR_MatTriangSlvUU_CaABbC_double(mb, n, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, beta, c, ldc, work, ind_base); break; case 2: /* Left scaling */ VBR_MatTriangSlvUU_CaDABbC_double(mb, n, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, beta, c, ldc, work, ind_base); break; case 3: /* Right scaling */ VBR_MatTriangSlvUU_CaADBbC_double(mb, n, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, beta, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. 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 */ VBR_VecTriangSlvLU_CATBC_double(mb, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CDATBC_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CATDBC_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CATBmC_double(mb, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CDATBmC_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CATDBmC_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CATB_double(mb, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CDATB_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CATDB_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CATBbC_double(mb, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CDATBbC_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CATDBbC_double(mb, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, beta, c, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } } else { /* alpha is general nonzero */ if (beta == 1) { switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CaATBC_double(mb, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CaDATBC_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CaATDBC_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CaATBmC_double(mb, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CaDATBmC_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CaATDBmC_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CaATB_double(mb, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CaDATB_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CaATDB_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, c, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */ } else { /* beta is general nonzero */ switch (unitd) { case 1: /* No scaling */ VBR_VecTriangSlvLU_CaATBbC_double(mb, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, beta, c, work, ind_base); break; case 2: /* Left scaling */ VBR_VecTriangSlvLU_CaDATBbC_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, beta, c, work, ind_base); break; case 3: /* Right scaling */ VBR_VecTriangSlvLU_CaATDBbC_double(mb, dv, alpha, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, beta, c, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. 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 */ VBR_MatTriangSlvLU_CATBC_double(mb, n, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ VBR_MatTriangSlvLU_CDATBC_double(mb, n, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ VBR_MatTriangSlvLU_CATDBC_double(mb, n, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */#if (0) } else if (beta == -1) { switch (unitd) { case 1: /* No scaling */ VBR_MatTriangSlvLU_CATBmC_double(mb, n, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; case 2: /* Left scaling */ VBR_MatTriangSlvLU_CDATBmC_double(mb, n, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ VBR_MatTriangSlvLU_CATDBmC_double(mb, n, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */#endif } else if (beta == 0) { switch (unitd) { case 1: /* No scaling */ VBR_MatTriangSlvLU_CATB_double(mb, n, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, ind_base); break; case 2: /* Left scaling */ VBR_MatTriangSlvLU_CDATB_double(mb, n, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, work, ind_base); break; case 3: /* Right scaling */ VBR_MatTriangSlvLU_CATDB_double(mb, n, dv, val, indx, bindx, rpntr, cpntr, bpntrb, bpntre, b, ldb, c, ldc, ind_base); break; default: printf("Invalid argument unitd in dvbrsm. Use 1-3. \n"); break; } /* end switch on unitd */⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -