Please go to api/linalg.finc for the API documentation. Since GSL follows the C convention for ordering of elements, all matrices must be set up and read out transposed.