From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001 From: Joseph Hunkeler Date: Wed, 8 Jul 2015 20:46:52 -0400 Subject: Initial commit --- noao/imred/vtel/imratio.x | 29 +++++++++++++++++++++++++++++ 1 file changed, 29 insertions(+) create mode 100644 noao/imred/vtel/imratio.x (limited to 'noao/imred/vtel/imratio.x') diff --git a/noao/imred/vtel/imratio.x b/noao/imred/vtel/imratio.x new file mode 100644 index 00000000..5586d204 --- /dev/null +++ b/noao/imred/vtel/imratio.x @@ -0,0 +1,29 @@ +# IMRATIO -- Divide two images and return the result in a third image. + +procedure imratio (numerator, denominator, ratio, xdim, ydim) + +real numerator[xdim, ydim] # input numerator +real denominator[xdim, ydim] # input denominator +real ratio[xdim, ydim] # output ratio image +int xdim, ydim # dimensions of the image + +int i +real ezero() +extern ezero() + +begin + do i = 1, ydim { + call arltr (denominator[1,i], xdim, 1E-10, 0.0) + call advzr (numerator[1,i], denominator[1,i], ratio[1,i], xdim, + ezero) + } +end + + +real procedure ezero (input) + +real input + +begin + return (0.0) +end -- cgit