%%	options

copyright owner	=	Dirk Krause
copyright year	=	2012-2014
license		=	bsd


%%	header

#ifdef __cplusplus
extern "C" {
#endif

/**	Coordinates transformation: X value.
	@param	ct	Transformation data.
	@param	x	X value.
	@param	ec	Pointer to error code variable, may be NULL.
	@return	Transformation result.
*/
double
dk3ct_2d_x(dk3_ct_2d_t const *ct, double x, int *ec);

/**	Coordinates transformation: Y value.
	@param	ct	Transformation data.
	@param	y	X value.
	@param	ec	Pointer to error code variable, may be NULL.
	@return	Transformation result.
*/
double
dk3ct_2d_y(dk3_ct_2d_t const *ct, double y, int *ec);

/**	Coordinates transformation: Radius.
	@param	ct	Transformation data.
	@param	r	X value.
	@param	ec	Pointer to error code variable, may be NULL.
	@return	Transformation result.
*/
double
dk3ct_2d_r(dk3_ct_2d_t const *ct, double r, int *ec);

#ifdef __cplusplus
}
#endif

%%	module

#include "dk3all.h"



$!trace-include


double
dk3ct_2d_x(dk3_ct_2d_t const *ct, double x, int *ec)
{
  double back = 0.0;
  $? "+ dk3ct_2d_x %lg", x
  if(ct) {
    back = dk3ma_d_add_ok(
      dk3ma_d_mul_ok(ct->mx, x, ec),
      ct->nx,
      ec
    );
  } $? "- dk3ct_2d_x %lg", back
  return back;
}



double
dk3ct_2d_y(dk3_ct_2d_t const *ct, double y, int *ec)
{
  double back = 0.0;
  $? "+ dk3ct_2d_y %lg", y
  if(ct) {
    back = dk3ma_d_add_ok(
      dk3ma_d_mul_ok(ct->my, y, ec),
      ct->ny,
      ec
    );
  } $? "- dk3ct_2d_y %lg", back
  return back;
}



double
dk3ct_2d_r(dk3_ct_2d_t const *ct, double r, int *ec)
{
  double back = 0.0;
  $? "+ dk3ct_2d_r %lg", r
  if(ct) {
    back = dk3ma_d_mul_ok(fabs(ct->mx), r, ec);
  } $? "- dk3ct_2d_r %lg", back
  return back;
}