1) Calculate decomposition of 
with the factorization step of the direct solver;
2) Calculate 
with the solving step applied to multiple right-hand sides;
3) Calculate 
as sparse-dense matrix-matrix multiplication;
4) Calculate 
as a difference.