#pragma rtGlobals=1		// Use modern global access method.

Function allowFCC(H)	// flags whether H=(hkl) is an allowed FCC reflection (not mixed)
	Wave H
	if ( (abs(H[0])+abs(H[1])+abs(H[2]))==0 )	// the (000) is not allowed
		return 0
	endif
	Variable sum = mod(abs(H[0]),2)+mod(abs(H[1]),2)+mod(abs(H[2]),2)
	sum = round(sum)
	return (sum==0 || sum==3)
End


Function allowBCC(H)	// flags whether H=(hkl) is an allowed BCC reflection (sum is even)
	Wave H
	return !mod(sum(H,-inf,inf),2)
End


//  This function suplanted by new built in Igor function,  MatrixDot(a,b)
//
//Function dot(a,b)		// vector dot product
//	Wave a,b
//	Variable i=0, n, sum
//	n = numpnts(a)
//	if (n!=numpnts(b))
//		Abort "in dot product, waves must be same "
//	endif
//
//	do
//		sum += a[i]*b[i]
//		i += 1
//	while (i<n)
//	return sum
//End


Function normV(a)	// returns the magnitude of vector a (a is not changed)
	Wave a
	return sqrt(MatrixDot(a,a))
//	Variable i=0, n, sum=0
//	n = numpnts(a)
//
//	do
//		sum += a[i]*a[i]
//		i += 1
//	while (i<n)
//	return sqrt(sum)
End


Function normalize(a)	// normalize a and return the initial magnitude
	Wave a
	Variable norm = sqrt(MatrixDot(a,a))
	a /= norm
	return norm
//	Variable i=0, n, norm=0
//	n = numpnts(a)
//
//	do
//		norm += a[i]*a[i]
//		i += 1
//	while (i<n)
//	norm = sqrt(norm)
//	if (norm<=0)
//		return 0
//	endif
//
//	i = 0
//	do
//		a[i] /= norm
//		i += 1
//	while (i<n)
//	return norm
End

Function cross(a,b,c)	// vector cross product
	Wave a,b,c
	c[2] = a[0]*b[1] - a[1]*b[0]
	c[0] = a[1]*b[2] - a[2]*b[1]
	c[1] = a[2]*b[0] - a[0]*b[2]
End


//Function printvec(w)		// print a vector to screen
//	Wave w
//	Variable i=0, n
//	n = numpnts(w)
//	printf "%s = {", NameOfWave(w)
//	do
//		if (i<(n-1))
//			printf "%g, ", w[i]
//		else
//			printf "%g}\r",w[i]
//		endif
//		i += 1
//	while (i<n)
//End