R/fmesher.R

In fmesher: Triangle Meshes and Related Geometry Tools

# # Sparse partial inverse
# #
# # @param A A symmetric matrix
# #
# # @export
# # @keywords internal
#
# fm_qinv <- function(A) {
#   A_C <- fm_as_dgCMatrix(A)
#   stopifnot(nrow(A_C) == ncol(A_C))
#   if (!identical(A_C, Matrix::t(A_C))) {
#     warning("Asymmetric matrix A detected, but only lower left triangle will be used.")
#   }
#   C_qinv(A_C)
# }