Function gmtapi_grid2matrix

src/gmt_api.c:6873–6912

}

GMT_LOCAL void *gmtapi_grid2matrix (struct GMTAPI_CTRL *API, struct GMT_GRID *In, struct GMT_MATRIX *Out) {
	bool alloc = (Out == NULL);
	openmp_int row, col;
	uint64_t ij, ij_M;
	double d;
	GMT_putfunction api_put_val = NULL;
	p_func_uint64_t GMT_2D_to_index = NULL;

	if (alloc) Out = gmtlib_create_matrix (API->GMT, 1U, 0);

	gmtapi_grdheader_to_matrixinfo (API->GMT, In->header, Out);
	if (alloc) {	/* Allocate the matrix itself */
		int error;
		Out->type = API->GMT->current.setting.export_type;
		Out->registration = In->header->registration;
		Out->shape = GMT_IS_ROW_FORMAT;	/* For now */
		Out->dim = (Out->shape == GMT_IS_ROW_FORMAT) ? Out->n_columns : Out->n_rows;	/* Matrix layout order */

		if ((error = gmtlib_alloc_univector (API->GMT, &(Out->data), Out->type, Out->n_rows * Out->n_columns)) != GMT_NOERROR) {
			gmt_M_free (API->GMT, Out);
			return_null (API, error);
		}
	}
	if ((GMT_2D_to_index = gmtapi_get_2d_to_index (API, Out->shape, GMT_GRID_IS_REAL)) == NULL) {
		if (alloc) gmt_M_free (API->GMT, Out);
		return_null (API, GMT_WRONG_MATRIX_SHAPE);
	}
	if ((api_put_val = gmtapi_select_put_function (API, Out->type)) == NULL) {
		if (alloc) gmt_M_free (API->GMT, Out);
		return_null (API, GMT_NOT_A_VALID_TYPE);
	}

	gmt_M_grd_loop (API->GMT, In, row, col, ij) {
		d = In->data[ij];
		ij_M = GMT_2D_to_index (row, col, Out->dim);
		api_put_val (&(Out->data), ij_M, d);
	}

	return Out;
}

GMT_LOCAL void *gmtapi_matrix2grid (struct GMTAPI_CTRL *API, struct GMT_MATRIX *In, struct GMT_GRID *Out) {
	bool alloc = (Out == NULL);