Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Eliminate Dirichlet boundary conditions differently
Slicing is very expensive operation on sparse matrices. Avoid slicing by calculating matrix product `K_red = C1*K*C1`, where C1 is diagonal matrix. This should be much cheaper.
- Loading branch information