Skip to content

Remove spam dependency #31

@bryo-han

Description

@bryo-han

Matrix::solve can solve triangular systems/do forward and backward substitution. So, spam may not be needed as a dependency.

To work with Matrix::solve, the triangular matrix must be of a suitable class. pedigreeTools::getLInv returns class dtCMatrix, which is suitable. A matrix imported from MatrixMarket format has to be converted using as(myMat, "dtCMatrix"). Not sure whether a matrix in spam format can be converted w/o going through a dense matrix.

Performance of Matrix::solve compares well with spam both with simulated and Labrador data.

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type
No fields configured for issues without a type.

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions