From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- sys/mwcs/gen/mwc1trand.x | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) create mode 100644 sys/mwcs/gen/mwc1trand.x (limited to 'sys/mwcs/gen/mwc1trand.x') diff --git a/sys/mwcs/gen/mwc1trand.x b/sys/mwcs/gen/mwc1trand.x new file mode 100644 index 00000000..af46e02d --- /dev/null +++ b/sys/mwcs/gen/mwc1trand.x @@ -0,0 +1,24 @@ +include "../mwcs.h" + +# MW_C1TRAN -- Optimized 1D coordinate transformation. + +double procedure mw_c1trand (a_ct, x) + +pointer a_ct #I pointer to CTRAN descriptor +double x #I coordinates in input system + +double y +pointer ct + +begin + # Get real or double version of descriptor. + ct = CT_D(a_ct) + + # Perform the transformation; LNR is a simple linear transformation. + if (CT_TYPE(ct) == LNR) { + return (Memd[CT_LTM(ct)] * x + Memd[CT_LTV(ct)]) + } else { + call mw_ctrand (a_ct, x, y, 1) + return (y) + } +end -- cgit