Fast dense solver for hierarchically off-diagonal low-rank matrices: Port from C++ by amiraa127