LINEAR++ library: AFF to LAPACK

◆ LaLinearSolve()

void LaLinearSolve ( LaTridiagFactDouble AF,
LaGenMatDouble &  X,
LaGenMatDouble &  B 
)
inline

Definition at line 123 of file trfd.h.

References dgttrs(), LaTridiagFactDouble::diag(), F77NAME, LaTridiagFactDouble::pivot(), and LaTridiagFactDouble::size().

125 {
126  char trans = 'N';
127  integer N = AF.size(), nrhs = X.size(1), ldb = B.size(0), info = 0;
128  double *DL = &AF.diag(-1)(0), *D = &AF.diag(0)(0),
129  *DU = &AF.diag(1)(0), *DU2 = &AF.diag(2)(0);
130 
131  X.inject(B);
132  F77NAME(dgttrs)(&trans, &N, &nrhs, DL, D, DU, DU2, &(AF.pivot()(0)),
133  &X(0,0), &ldb, &info);
134 }
void F77NAME() dgttrs(char *trans, integer *N, integer *nrhs, doublereal *DL, doublereal *D, doublereal *DU, doublereal *DU2, integer *ipiv, doublereal *b, integer *ldb, integer *info)
LaVectorDouble diag(int)
Definition: trfd.h:67
LaVectorLongInt & pivot()
Definition: trfd.h:30
long int integer
Definition: f77lapack.h:61
#define F77NAME(x)
Definition: arch.h:17
Here is the call graph for this function: