Assertion fail

[robertopsouto]

I'm using Strumpack compiled in Debug mode, and trying to solve a system with a unsymmetric matrix using Fortran code, similar to the example ./examples/sparse/fexample.f90.
However, it returns the error Program received signal SIGABRT: Process abort signal. when executing the STRUMPACK_factor(S) function. I'm attaching the code where this error occurs:

strumpack_unsym.f90.md

By using Strumpack compiled in Release mode, the error is Program received signal SIGSEGV: Segmentation fault - invalid memory reference., at the same local in the code.

Please, what could I be forgetting to do to make the code work correctly?
Thank you in advance.

# Initializing STRUMPACK
# using 1 OpenMP thread(s)
 n+1, ia(n+1):         4001      397501
# matching job: none
# matrix equilibration, r_cond = 0.144869 , c_cond = 1 , type = N
# Matrix padded with zeros to get symmetric pattern.
# Number of nonzeros increased from 397,500 to 496,355.
# initial matrix:
#   - number of unknowns = 4,000
#   - number of nonzeros = 496,355
# nested dissection reordering:
#   - Natural reordering
#   - strategy parameter = 8
#   - number of separators = 3
#   - number of levels = 2
#   - nd time = 0.0563012
#   - symmetrization time = 0.167539
# symbolic factorization:
#   - nr of dense Frontal matrices = 2
#   - nr of BLR Frontal matrices = 1
#   - symb-factor time = 6.5885e-05
#   - sep-reorder time = 0.25354
# multifrontal factorization:
#   - estimated memory usage (exact solver) = 127.874 MB
#   - minimum pivot, sqrt(eps)*|A|_1 = 1.05367e-08
#   - replacing of small pivots is not enabled
ltsn: /home/rpsouto/programming/strumpack/github/STRUMPACK/src/dense/DenseMatrix.hpp:295: scalar_t& strumpack::DenseMatrix<scalar_t>::operator()(std::size_t, std::size_t) [with scalar_t = double; std::size_t = long unsigned int]: Assertion `i<=rows() && j<=cols()' failed.

Program received signal SIGABRT: Process abort signal.