Initial commit

1 files changed, 522 insertions, 0 deletions

diff --git a/noao/onedspec/odcombine/srcwt/generic/icaverage.x b/noao/onedspec/odcombine/srcwt/generic/icaverage.x
new file mode 100644
index 00000000..4b464e91
--- /dev/null
+++ b/noao/onedspec/odcombine/srcwt/generic/icaverage.x
@@ -0,0 +1,522 @@
+# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.
+
+include	<imhdr.h>
+include	"../icombine.h"
+
+
+# IC_AVERAGE -- Compute the average (or summed) image line.
+# Options include a weighted average/sum.
+
+procedure ic_averages (d, m, n, wts, w, npts, doblank, doaverage, average)
+
+pointer	d[ARB]			# Data pointers
+pointer	m[ARB]			# Image ID pointers
+int	n[npts]			# Number of points
+real	wts[ARB]		# Weights
+pointer	w[ARB]			# Weight data pointers
+int	npts			# Number of output points per line
+int	doblank			# Set blank values?
+int	doaverage		# Do average?
+real	average[npts]		# Average (returned)
+
+int	i, j, k
+real	sumwt, wt
+real	sum
+
+include	"../icombine.com"
+
+begin
+	# If no data has been excluded do the average/sum without checking
+	# the number of points and using the fact that the weights are
+	# normalized.  If all the data has been excluded set the average/sum
+	# to the blank value if requested.
+
+	if (dflag == D_ALL && w[1] == NULL) {
+	    if (dowts) {
+		do i = 1, npts {
+		    k = i - 1
+		    wt = wts[Memi[m[1]+k]]
+		    sum = Mems[d[1]+k] * wt
+		    do j = 2, n[i] {
+			wt = wts[Memi[m[j]+k]]
+			sum = sum + Mems[d[j]+k] * wt
+		    }
+		    average[i] = sum
+		}
+	    } else {
+		do i = 1, npts {
+		    k = i - 1
+		    sum = Mems[d[1]+k]
+		    do j = 2, n[i]
+			sum = sum + Mems[d[j]+k]
+		    if (doaverage == YES)
+			average[i] = sum / n[i]
+		    else
+			average[i] = sum
+		}
+	    }
+	} else if (dflag == D_NONE) {
+	    if (doblank == YES) {
+		do i = 1, npts
+		    average[i] = blank
+	    }
+	} else {
+	    if (dowts) {
+		if (w[1] == NULL) {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = wts[Memi[m[1]+k]]
+			    sum = Mems[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = wts[Memi[m[j]+k]]
+				sum = sum + Mems[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Mems[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Mems[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		} else {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = Memr[w[Memi[m[1]+k]]+k]
+			    sum = Mems[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = Memr[w[Memi[m[j]+k]]+k]
+				sum = sum + Mems[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Mems[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Mems[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		}
+	    } else {
+		do i = 1, npts {
+		    if (n[i] > 0) {
+			k = i - 1
+			sum = Mems[d[1]+k]
+			do j = 2, n[i]
+			    sum = sum + Mems[d[j]+k]
+			if (doaverage == YES)
+			    average[i] = sum / n[i]
+			else
+			    average[i] = sum
+		    } else if (doblank == YES)
+			average[i] = blank
+		}
+	    }
+	}
+end
+
+# IC_AVERAGE -- Compute the average (or summed) image line.
+# Options include a weighted average/sum.
+
+procedure ic_averagei (d, m, n, wts, w, npts, doblank, doaverage, average)
+
+pointer	d[ARB]			# Data pointers
+pointer	m[ARB]			# Image ID pointers
+int	n[npts]			# Number of points
+real	wts[ARB]		# Weights
+pointer	w[ARB]			# Weight data pointers
+int	npts			# Number of output points per line
+int	doblank			# Set blank values?
+int	doaverage		# Do average?
+real	average[npts]		# Average (returned)
+
+int	i, j, k
+real	sumwt, wt
+real	sum
+
+include	"../icombine.com"
+
+begin
+	# If no data has been excluded do the average/sum without checking
+	# the number of points and using the fact that the weights are
+	# normalized.  If all the data has been excluded set the average/sum
+	# to the blank value if requested.
+
+	if (dflag == D_ALL && w[1] == NULL) {
+	    if (dowts) {
+		do i = 1, npts {
+		    k = i - 1
+		    wt = wts[Memi[m[1]+k]]
+		    sum = Memi[d[1]+k] * wt
+		    do j = 2, n[i] {
+			wt = wts[Memi[m[j]+k]]
+			sum = sum + Memi[d[j]+k] * wt
+		    }
+		    average[i] = sum
+		}
+	    } else {
+		do i = 1, npts {
+		    k = i - 1
+		    sum = Memi[d[1]+k]
+		    do j = 2, n[i]
+			sum = sum + Memi[d[j]+k]
+		    if (doaverage == YES)
+			average[i] = sum / n[i]
+		    else
+			average[i] = sum
+		}
+	    }
+	} else if (dflag == D_NONE) {
+	    if (doblank == YES) {
+		do i = 1, npts
+		    average[i] = blank
+	    }
+	} else {
+	    if (dowts) {
+		if (w[1] == NULL) {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = wts[Memi[m[1]+k]]
+			    sum = Memi[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = wts[Memi[m[j]+k]]
+				sum = sum + Memi[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Memi[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Memi[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		} else {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = Memr[w[Memi[m[1]+k]]+k]
+			    sum = Memi[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = Memr[w[Memi[m[j]+k]]+k]
+				sum = sum + Memi[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Memi[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Memi[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		}
+	    } else {
+		do i = 1, npts {
+		    if (n[i] > 0) {
+			k = i - 1
+			sum = Memi[d[1]+k]
+			do j = 2, n[i]
+			    sum = sum + Memi[d[j]+k]
+			if (doaverage == YES)
+			    average[i] = sum / n[i]
+			else
+			    average[i] = sum
+		    } else if (doblank == YES)
+			average[i] = blank
+		}
+	    }
+	}
+end
+
+# IC_AVERAGE -- Compute the average (or summed) image line.
+# Options include a weighted average/sum.
+
+procedure ic_averager (d, m, n, wts, w, npts, doblank, doaverage, average)
+
+pointer	d[ARB]			# Data pointers
+pointer	m[ARB]			# Image ID pointers
+int	n[npts]			# Number of points
+real	wts[ARB]		# Weights
+pointer	w[ARB]			# Weight data pointers
+int	npts			# Number of output points per line
+int	doblank			# Set blank values?
+int	doaverage		# Do average?
+real	average[npts]		# Average (returned)
+
+int	i, j, k
+real	sumwt, wt
+real	sum
+
+include	"../icombine.com"
+
+begin
+	# If no data has been excluded do the average/sum without checking
+	# the number of points and using the fact that the weights are
+	# normalized.  If all the data has been excluded set the average/sum
+	# to the blank value if requested.
+
+	if (dflag == D_ALL && w[1] == NULL) {
+	    if (dowts) {
+		do i = 1, npts {
+		    k = i - 1
+		    wt = wts[Memi[m[1]+k]]
+		    sum = Memr[d[1]+k] * wt
+		    do j = 2, n[i] {
+			wt = wts[Memi[m[j]+k]]
+			sum = sum + Memr[d[j]+k] * wt
+		    }
+		    average[i] = sum
+		}
+	    } else {
+		do i = 1, npts {
+		    k = i - 1
+		    sum = Memr[d[1]+k]
+		    do j = 2, n[i]
+			sum = sum + Memr[d[j]+k]
+		    if (doaverage == YES)
+			average[i] = sum / n[i]
+		    else
+			average[i] = sum
+		}
+	    }
+	} else if (dflag == D_NONE) {
+	    if (doblank == YES) {
+		do i = 1, npts
+		    average[i] = blank
+	    }
+	} else {
+	    if (dowts) {
+		if (w[1] == NULL) {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = wts[Memi[m[1]+k]]
+			    sum = Memr[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = wts[Memi[m[j]+k]]
+				sum = sum + Memr[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Memr[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Memr[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		} else {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = Memr[w[Memi[m[1]+k]]+k]
+			    sum = Memr[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = Memr[w[Memi[m[j]+k]]+k]
+				sum = sum + Memr[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Memr[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Memr[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		}
+	    } else {
+		do i = 1, npts {
+		    if (n[i] > 0) {
+			k = i - 1
+			sum = Memr[d[1]+k]
+			do j = 2, n[i]
+			    sum = sum + Memr[d[j]+k]
+			if (doaverage == YES)
+			    average[i] = sum / n[i]
+			else
+			    average[i] = sum
+		    } else if (doblank == YES)
+			average[i] = blank
+		}
+	    }
+	}
+end
+
+# IC_AVERAGE -- Compute the average (or summed) image line.
+# Options include a weighted average/sum.
+
+procedure ic_averaged (d, m, n, wts, w, npts, doblank, doaverage, average)
+
+pointer	d[ARB]			# Data pointers
+pointer	m[ARB]			# Image ID pointers
+int	n[npts]			# Number of points
+real	wts[ARB]		# Weights
+pointer	w[ARB]			# Weight data pointers
+int	npts			# Number of output points per line
+int	doblank			# Set blank values?
+int	doaverage		# Do average?
+double	average[npts]		# Average (returned)
+
+int	i, j, k
+real	sumwt, wt
+double	sum
+
+include	"../icombine.com"
+
+begin
+	# If no data has been excluded do the average/sum without checking
+	# the number of points and using the fact that the weights are
+	# normalized.  If all the data has been excluded set the average/sum
+	# to the blank value if requested.
+
+	if (dflag == D_ALL && w[1] == NULL) {
+	    if (dowts) {
+		do i = 1, npts {
+		    k = i - 1
+		    wt = wts[Memi[m[1]+k]]
+		    sum = Memd[d[1]+k] * wt
+		    do j = 2, n[i] {
+			wt = wts[Memi[m[j]+k]]
+			sum = sum + Memd[d[j]+k] * wt
+		    }
+		    average[i] = sum
+		}
+	    } else {
+		do i = 1, npts {
+		    k = i - 1
+		    sum = Memd[d[1]+k]
+		    do j = 2, n[i]
+			sum = sum + Memd[d[j]+k]
+		    if (doaverage == YES)
+			average[i] = sum / n[i]
+		    else
+			average[i] = sum
+		}
+	    }
+	} else if (dflag == D_NONE) {
+	    if (doblank == YES) {
+		do i = 1, npts
+		    average[i] = blank
+	    }
+	} else {
+	    if (dowts) {
+		if (w[1] == NULL) {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = wts[Memi[m[1]+k]]
+			    sum = Memd[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = wts[Memi[m[j]+k]]
+				sum = sum + Memd[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Memd[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Memd[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		} else {
+		    do i = 1, npts {
+			if (n[i] > 0) {
+			    k = i - 1
+			    wt = Memr[w[Memi[m[1]+k]]+k]
+			    sum = Memd[d[1]+k] * wt
+			    sumwt = wt
+			    do j = 2, n[i] {
+				wt = Memr[w[Memi[m[j]+k]]+k]
+				sum = sum + Memd[d[j]+k] * wt
+				sumwt = sumwt + wt
+			    }
+			    if (doaverage == YES) {
+				if (sumwt > 0)
+				    average[i] = sum / sumwt
+				else {
+				    sum = Memd[d[1]+k]
+				    do j = 2, n[i]
+					sum = sum + Memd[d[j]+k]
+				    average[i] = sum / n[i]
+				}
+			    } else
+				average[i] = sum
+			} else if (doblank == YES)
+			    average[i] = blank
+		    }
+		}
+	    } else {
+		do i = 1, npts {
+		    if (n[i] > 0) {
+			k = i - 1
+			sum = Memd[d[1]+k]
+			do j = 2, n[i]
+			    sum = sum + Memd[d[j]+k]
+			if (doaverage == YES)
+			    average[i] = sum / n[i]
+			else
+			    average[i] = sum
+		    } else if (doblank == YES)
+			average[i] = blank
+		}
+	    }
+	}
+end
+