/* 2D transforms */
CAMLprim value
ml_gsl_wavelet2d_transform_matrix (value w, value ordering,
value dir, value a, value ws)
{
_DECLARE_MATRIX(a);
_CONVERT_MATRIX(a);
if (Int_val (ordering) == 0)
gsl_wavelet2d_transform_matrix (Wavelet_val (w), &m_a,
gsl_direction_val (dir), WS_val (ws));
else
gsl_wavelet2d_nstransform_matrix (Wavelet_val (w), &m_a,
gsl_direction_val (dir), WS_val (ws));
return Val_unit;
}
/**
* C++ version of gsl_wavelet2d_transform_matrix().
* @param a A matrix
* @param dir Direction: gsl_wavelet-forward() or gsl_wavelet_backward()
* @param work Supply a workspace of appropriate (undocumented) size
* @return Error code on failure
*/
int transform_matrix( matrix& a, gsl_wavelet_direction dir,
wavelet::workspace& work ) const {
return gsl_wavelet2d_transform_matrix( get(), a.get(), dir, work.get() ); }
