cal/cal/normcomp.cal

{
        Compute overlap between two normal distribution functions,
        defined by mean1, mean2, stdev1 and stdev2.

        2/10/95         Greg Ward

        relies on gauss.cal and norm.cal
}
min(a,b) : if(a-b, b, a);
max(a,b) : if(a-b, a, b);

Z1(u) = Z((u - mean1)/stdev1) / stdev1;
Z2(u) = Z((u - mean2)/stdev2) / stdev2;

f(u) = min(Z1(u),Z2(u));

xmin = min(mean1-3*stdev1, mean2-3*stdev2);
xmax = max(mean1+3*stdev1, mean2+3*stdev2);
Revision:	1.1
Committed:	Sat Feb 22 02:07:21 2003 UTC (21 years, 2 months ago) by greg
Branch:	MAIN
CVS Tags:	rad5R2, rad4R2P2, rad5R0, rad5R1, rad3R7P2, rad3R7P1, rad4R2, rad4R1, rad4R0, rad3R5, rad3R6, rad3R6P1, rad3R8, rad3R9, rad4R2P1
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	Content
1	{
2	Compute overlap between two normal distribution functions,
3	defined by mean1, mean2, stdev1 and stdev2.
4
5	2/10/95 Greg Ward
6
7	relies on gauss.cal and norm.cal
8	}
9	min(a,b) : if(a-b, b, a);
10	max(a,b) : if(a-b, a, b);
11
12	Z1(u) = Z((u - mean1)/stdev1) / stdev1;
13	Z2(u) = Z((u - mean2)/stdev2) / stdev2;
14
15	f(u) = min(Z1(u),Z2(u));
16
17	xmin = min(mean1-3stdev1, mean2-3stdev2);
18	xmax = max(mean1+3stdev1, mean2+3stdev2);