A sparse linear equations solver for 3D–regular domains