2 #include "../gmx_blas.h"
3 #include "../gmx_lapack.h"
4 #include "lapack_limits.h"
6 #include "types/simple.h"
9 F77_FUNC(dtrtri,DTRTRI)(const char *uplo,
16 int a_dim1, a_offset, i__1, i__3, i__4, i__5;
25 a_offset = 1 + a_dim1;
29 upper = (*uplo=='U' || *uplo=='u');
30 nounit = (*diag=='N' || *diag=='n');
43 for (*info = 1; *info <= i__1; ++(*info)) {
44 if (fabs(a[*info + *info * a_dim1])<GMX_DOUBLE_MIN) {
51 nb = DTRTRI_BLOCKSIZE;
52 if (nb <= 1 || nb >= *n) {
54 F77_FUNC(dtrti2,DTRTI2)(uplo, diag, n, &a[a_offset], lda, info);
61 for (j = 1; i__3 < 0 ? j >= i__1 : j <= i__1; j += i__3) {
62 i__4 = nb, i__5 = *n - j + 1;
63 jb = (i__4<i__5) ? i__4 : i__5;
66 F77_FUNC(dtrmm,DTRMM)("Left", "Upper", "No transpose", diag, &i__4, &jb, &
67 c_b18, &a[a_offset], lda, &a[j * a_dim1 + 1], lda);
69 F77_FUNC(dtrsm,DTRSM)("Right", "Upper", "No transpose", diag, &i__4, &jb, &
70 c_b22, &a[j + j * a_dim1], lda, &a[j * a_dim1 + 1],
73 F77_FUNC(dtrti2,DTRTI2)("Upper", diag, &jb, &a[j + j * a_dim1], lda, info);
77 nn = (*n - 1) / nb * nb + 1;
79 for (j = nn; i__3 < 0 ? j >= 1 : j <= 1; j += i__3) {
80 i__1 = nb, i__4 = *n - j + 1;
81 jb = (i__1<i__4) ? i__1 : i__4;
84 i__1 = *n - j - jb + 1;
85 F77_FUNC(dtrmm,DTRMM)("Left", "Lower", "No transpose", diag, &i__1, &jb,
86 &c_b18, &a[j + jb + (j + jb) * a_dim1], lda, &a[j
87 + jb + j * a_dim1], lda);
88 i__1 = *n - j - jb + 1;
89 F77_FUNC(dtrsm,DTRSM)("Right", "Lower", "No transpose", diag, &i__1, &jb,
90 &c_b22, &a[j + j * a_dim1], lda, &a[j + jb + j *
94 F77_FUNC(dtrti2,DTRTI2)("Lower", diag, &jb, &a[j + j * a_dim1], lda, info);