The first derivative matrix is calculated using the ray paths computed after the source
locations in the 3D model. Each element of the matrix. is equal to the time
deviation along the i-th ray due to a unit velocity perturbation in the j-th node/block.
Inversion of the entire sparse A matrix is performed using an iterative LSQR code