Context Navigation

source: branches/HeuristicLab.Problems.GaussianProcessTuning/ILNumerics.2.14.4735.573/Functions/builtin/or.cs @ 9102

Visit:

Last change on this file since 9102 was 9102, checked in by gkronber, 12 years ago
#1967: ILNumerics source for experimentation
File size: 94.7 KB

Rev	Line
[9102]	1	///
	2	/// This file is part of ILNumerics Community Edition.
	3	///
	4	/// ILNumerics Community Edition - high performance computing for applications.
	5	/// Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net
	6	///
	7	/// ILNumerics Community Edition is free software: you can redistribute it and/or modify
	8	/// it under the terms of the GNU General Public License version 3 as published by
	9	/// the Free Software Foundation.
	10	///
	11	/// ILNumerics Community Edition is distributed in the hope that it will be useful,
	12	/// but WITHOUT ANY WARRANTY; without even the implied warranty of
	13	/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
	14	/// GNU General Public License for more details.
	15	///
	16	/// You should have received a copy of the GNU General Public License
	17	/// along with ILNumerics Community Edition. See the file License.txt in the root
	18	/// of your distribution package. If not, see <http://www.gnu.org/licenses/>.
	19	///
	20	/// In addition this software uses the following components and/or licenses:
	21	///
	22	/// =================================================================================
	23	/// The Open Toolkit Library License
	24	///
	25	/// Copyright (c) 2006 - 2009 the Open Toolkit library.
	26	///
	27	/// Permission is hereby granted, free of charge, to any person obtaining a copy
	28	/// of this software and associated documentation files (the "Software"), to deal
	29	/// in the Software without restriction, including without limitation the rights to
	30	/// use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
	31	/// the Software, and to permit persons to whom the Software is furnished to do
	32	/// so, subject to the following conditions:
	33	///
	34	/// The above copyright notice and this permission notice shall be included in all
	35	/// copies or substantial portions of the Software.
	36	///
	37	/// =================================================================================
	38	///
	39
	40	using System;
	41	using System.Collections.Generic;
	42	using System.Text;
	43	using System.Runtime.InteropServices;
	44	using ILNumerics.Storage;
	45	using ILNumerics.Misc;
	46	using ILNumerics.Native;
	47	using ILNumerics.Exceptions;
	48
	49
	50
	51	namespace ILNumerics {
	52
	53	public partial class ILMath {
	54
	55
	56
	57
	58	#region HYCALPER AUTO GENERATED CODE
	59
	60	/// <summary>Elementwise logical 'or' operator</summary>
	61	/// <param name="A">Input array A</param>
	62	/// <param name="B">Input array B</param>
	63	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	64	/// <remarks><para>On empty input an empty array will be returned.</para>
	65	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	66	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	67	public unsafe static ILRetLogical or(ILInArray<Int64> A, ILInArray<Int64> B) {
	68	using (ILScope.Enter(A, B)) {
	69	int outLen;
	70	BinOpItMode mode;
	71	byte[] retArr;
	72
	73	Int64[] arrA = A.GetArrayForRead();
	74
	75	Int64[] arrB = B.GetArrayForRead();
	76	ILSize outDims;
	77	if (A.IsScalar) {
	78	if (B.IsScalar) {
	79	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	80	} else if (B.IsEmpty) {
	81	return new ILRetLogical(B.Size);
	82	} else {
	83	outLen = B.S.NumberOfElements;
	84	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	85	mode = BinOpItMode.SAN;
	86	}
	87	outDims = B.Size;
	88	} else {
	89	outDims = A.Size;
	90	if (B.IsScalar) {
	91	if (A.IsEmpty) {
	92	return new ILRetLogical(A.Size);
	93	}
	94	outLen = A.S.NumberOfElements;
	95	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	96	mode = BinOpItMode.ASN;
	97	} else {
	98	// array + array
	99	if (!A.Size.IsSameSize(B.Size)) {
	100	return orEx(A,B);
	101	}
	102	outLen = A.S.NumberOfElements;
	103	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	104	mode = BinOpItMode.AAN;
	105	}
	106	}
	107	int workerCount = 1;
	108	Action<object> worker = data => {
	109	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	110	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	111	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	112	Int64 scalar;
	113	cLast = cp + range.Item2;
	114	#region loops
	115	switch (mode) {
	116	case BinOpItMode.AAN:
	117	Int64* ap = ((Int64*)range.Item3 + range.Item1);
	118	Int64* bp = ((Int64*)range.Item4 + range.Item1);
	119	while (cp < cLast) {
	120	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	121	cp++;
	122	ap++;
	123	bp++;
	124	}
	125	break;
	126	case BinOpItMode.ASN:
	127	ap = ((Int64*)range.Item3 + range.Item1);
	128	scalar = ((Int64)range.Item4);
	129	while (cp < cLast) {
	130	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	131	cp++;
	132	ap++;
	133
	134	}
	135	break;
	136	case BinOpItMode.SAN:
	137	scalar = ((Int64)range.Item3);
	138	bp = ((Int64*)range.Item4 + range.Item1);
	139	while (cp < cLast) {
	140	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	141	cp++;
	142	bp++;
	143	}
	144	break;
	145	default:
	146	break;
	147	}
	148	#endregion
	149	System.Threading.Interlocked.Decrement(ref workerCount);
	150	};
	151
	152	#region do the work
	153	fixed (Int64* arrAP = arrA)
	154	fixed (Int64* arrBP = arrB)
	155	fixed (byte* retArrP = retArr) {
	156	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	157	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	158	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	159	workItemLength = outLen / workItemCount;
	160	} else {
	161	workItemLength = outLen / 2;
	162	workItemCount = 2;
	163	}
	164	} else {
	165	workItemLength = outLen;
	166	workItemCount = 1;
	167	}
	168
	169	for (; i < workItemCount - 1; i++) {
	170	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	171	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	172	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	173	System.Threading.Interlocked.Increment(ref workerCount);
	174	ILThreadPool.QueueUserWorkItem(i, worker, range);
	175	}
	176	// the last (or may the only) chunk is done right here
	177	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	178	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	179	ILThreadPool.Wait4Workers(ref workerCount);
	180	}
	181
	182	#endregion
	183	return new ILRetLogical(retArr, outDims);
	184	}
	185	}
	186
	187	private static unsafe ILRetLogical orEx(ILInArray<Int64> A, ILInArray<Int64> B) {
	188	using (ILScope.Enter(A, B)) {
	189	#region parameter checking
	190	if (isnull(A) \|\| isnull(B))
	191	return new ILRetLogical(ILSize.Empty00);
	192	if (A.IsEmpty) {
	193	return new ILRetLogical(B.S);
	194	} else if (B.IsEmpty) {
	195	return new ILRetLogical(A.S);
	196	}
	197	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	198	// return add(A,B);
	199	int dim = -1;
	200	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	201	if (A.S[l] != B.S[l]) {
	202	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	203	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	204	}
	205	dim = l;
	206	}
	207	}
	208	if (dim > 1)
	209	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	210	#endregion
	211
	212	#region parameter preparation
	213	byte[] retArr;
	214
	215	Int64[] arrA = A.GetArrayForRead();
	216
	217	Int64[] arrB = B.GetArrayForRead();
	218	ILSize outDims;
	219	BinOptItExMode mode;
	220	int arrInc = 0;
	221	int arrStepInc = 0;
	222	int dimLen = 0;
	223	if (A.IsVector) {
	224	outDims = B.S;
	225	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	226	mode = BinOptItExMode.VAN;
	227	dimLen = A.Length;
	228	} else if (B.IsVector) {
	229	outDims = A.S;
	230	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	231	mode = BinOptItExMode.AVN;
	232	dimLen = B.Length;
	233	} else {
	234	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	235	}
	236	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	237	arrStepInc = outDims.SequentialIndexDistance(dim);
	238	#endregion
	239
	240	#region worker loops definition
	241	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	242	int workerCount = 1;
	243	Action<object> worker = data => {
	244	// expects: iStart, iLen, ap, bp, cp
	245	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	246	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	247
	248	Int64* ap;
	249
	250	Int64* bp;
	251	byte* cp;
	252	switch (mode) {
	253	case BinOptItExMode.VAN:
	254	for (int s = 0; s < range.Item2; s++) {
	255	ap = (Int64*)range.Item3;
	256	bp = (Int64)range.Item4 + range.Item1 + s arrStepInc; ;
	257	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	258	for (int l = 0; l < dimLen; l++) {
	259	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	260	ap++;
	261	bp += arrInc;
	262	cp += arrInc;
	263	}
	264	}
	265	break;
	266	case BinOptItExMode.AVN:
	267	for (int s = 0; s < range.Item2; s++) {
	268	ap = (Int64)range.Item3 + range.Item1 + s arrStepInc;
	269	bp = (Int64*)range.Item4;
	270	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	271	for (int l = 0; l < dimLen; l++) {
	272	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	273	ap += arrInc;
	274	bp++;
	275	cp += arrInc;
	276	}
	277	}
	278	break;
	279	}
	280	System.Threading.Interlocked.Decrement(ref workerCount);
	281	};
	282	#endregion
	283
	284	#region work distribution
	285	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	286	int outLen = outDims.NumberOfElements;
	287	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	288	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	289	workItemLength = outLen / dimLen / workItemCount;
	290	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	291	} else {
	292	workItemLength = outLen / dimLen / 2;
	293	workItemCount = 2;
	294	}
	295	} else {
	296	workItemLength = outLen / dimLen;
	297	workItemCount = 1;
	298	}
	299
	300	fixed ( Int64* arrAP = arrA)
	301	fixed ( Int64* arrBP = arrB)
	302	fixed (byte* retArrP = retArr) {
	303
	304	for (; i < workItemCount - 1; i++) {
	305	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	306	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	307	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	308	System.Threading.Interlocked.Increment(ref workerCount);
	309	ILThreadPool.QueueUserWorkItem(i, worker, range);
	310	}
	311	// the last (or may the only) chunk is done right here
	312	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	313	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	314	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	315
	316	ILThreadPool.Wait4Workers(ref workerCount);
	317	}
	318	#endregion
	319
	320	return new ILRetLogical(retStorage);
	321	}
	322	}
	323
	324	/// <summary>Elementwise logical 'or' operator</summary>
	325	/// <param name="A">Input array A</param>
	326	/// <param name="B">Input array B</param>
	327	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	328	/// <remarks><para>On empty input an empty array will be returned.</para>
	329	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	330	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	331	public unsafe static ILRetLogical or(ILInArray<Int32> A, ILInArray<Int32> B) {
	332	using (ILScope.Enter(A, B)) {
	333	int outLen;
	334	BinOpItMode mode;
	335	byte[] retArr;
	336
	337	Int32[] arrA = A.GetArrayForRead();
	338
	339	Int32[] arrB = B.GetArrayForRead();
	340	ILSize outDims;
	341	if (A.IsScalar) {
	342	if (B.IsScalar) {
	343	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	344	} else if (B.IsEmpty) {
	345	return new ILRetLogical(B.Size);
	346	} else {
	347	outLen = B.S.NumberOfElements;
	348	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	349	mode = BinOpItMode.SAN;
	350	}
	351	outDims = B.Size;
	352	} else {
	353	outDims = A.Size;
	354	if (B.IsScalar) {
	355	if (A.IsEmpty) {
	356	return new ILRetLogical(A.Size);
	357	}
	358	outLen = A.S.NumberOfElements;
	359	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	360	mode = BinOpItMode.ASN;
	361	} else {
	362	// array + array
	363	if (!A.Size.IsSameSize(B.Size)) {
	364	return orEx(A,B);
	365	}
	366	outLen = A.S.NumberOfElements;
	367	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	368	mode = BinOpItMode.AAN;
	369	}
	370	}
	371	int workerCount = 1;
	372	Action<object> worker = data => {
	373	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	374	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	375	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	376	Int32 scalar;
	377	cLast = cp + range.Item2;
	378	#region loops
	379	switch (mode) {
	380	case BinOpItMode.AAN:
	381	Int32* ap = ((Int32*)range.Item3 + range.Item1);
	382	Int32* bp = ((Int32*)range.Item4 + range.Item1);
	383	while (cp < cLast) {
	384	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	385	cp++;
	386	ap++;
	387	bp++;
	388	}
	389	break;
	390	case BinOpItMode.ASN:
	391	ap = ((Int32*)range.Item3 + range.Item1);
	392	scalar = ((Int32)range.Item4);
	393	while (cp < cLast) {
	394	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	395	cp++;
	396	ap++;
	397
	398	}
	399	break;
	400	case BinOpItMode.SAN:
	401	scalar = ((Int32)range.Item3);
	402	bp = ((Int32*)range.Item4 + range.Item1);
	403	while (cp < cLast) {
	404	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	405	cp++;
	406	bp++;
	407	}
	408	break;
	409	default:
	410	break;
	411	}
	412	#endregion
	413	System.Threading.Interlocked.Decrement(ref workerCount);
	414	};
	415
	416	#region do the work
	417	fixed (Int32* arrAP = arrA)
	418	fixed (Int32* arrBP = arrB)
	419	fixed (byte* retArrP = retArr) {
	420	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	421	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	422	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	423	workItemLength = outLen / workItemCount;
	424	} else {
	425	workItemLength = outLen / 2;
	426	workItemCount = 2;
	427	}
	428	} else {
	429	workItemLength = outLen;
	430	workItemCount = 1;
	431	}
	432
	433	for (; i < workItemCount - 1; i++) {
	434	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	435	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	436	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	437	System.Threading.Interlocked.Increment(ref workerCount);
	438	ILThreadPool.QueueUserWorkItem(i, worker, range);
	439	}
	440	// the last (or may the only) chunk is done right here
	441	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	442	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	443	ILThreadPool.Wait4Workers(ref workerCount);
	444	}
	445
	446	#endregion
	447	return new ILRetLogical(retArr, outDims);
	448	}
	449	}
	450
	451	private static unsafe ILRetLogical orEx(ILInArray<Int32> A, ILInArray<Int32> B) {
	452	using (ILScope.Enter(A, B)) {
	453	#region parameter checking
	454	if (isnull(A) \|\| isnull(B))
	455	return new ILRetLogical(ILSize.Empty00);
	456	if (A.IsEmpty) {
	457	return new ILRetLogical(B.S);
	458	} else if (B.IsEmpty) {
	459	return new ILRetLogical(A.S);
	460	}
	461	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	462	// return add(A,B);
	463	int dim = -1;
	464	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	465	if (A.S[l] != B.S[l]) {
	466	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	467	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	468	}
	469	dim = l;
	470	}
	471	}
	472	if (dim > 1)
	473	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	474	#endregion
	475
	476	#region parameter preparation
	477	byte[] retArr;
	478
	479	Int32[] arrA = A.GetArrayForRead();
	480
	481	Int32[] arrB = B.GetArrayForRead();
	482	ILSize outDims;
	483	BinOptItExMode mode;
	484	int arrInc = 0;
	485	int arrStepInc = 0;
	486	int dimLen = 0;
	487	if (A.IsVector) {
	488	outDims = B.S;
	489	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	490	mode = BinOptItExMode.VAN;
	491	dimLen = A.Length;
	492	} else if (B.IsVector) {
	493	outDims = A.S;
	494	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	495	mode = BinOptItExMode.AVN;
	496	dimLen = B.Length;
	497	} else {
	498	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	499	}
	500	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	501	arrStepInc = outDims.SequentialIndexDistance(dim);
	502	#endregion
	503
	504	#region worker loops definition
	505	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	506	int workerCount = 1;
	507	Action<object> worker = data => {
	508	// expects: iStart, iLen, ap, bp, cp
	509	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	510	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	511
	512	Int32* ap;
	513
	514	Int32* bp;
	515	byte* cp;
	516	switch (mode) {
	517	case BinOptItExMode.VAN:
	518	for (int s = 0; s < range.Item2; s++) {
	519	ap = (Int32*)range.Item3;
	520	bp = (Int32)range.Item4 + range.Item1 + s arrStepInc; ;
	521	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	522	for (int l = 0; l < dimLen; l++) {
	523	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	524	ap++;
	525	bp += arrInc;
	526	cp += arrInc;
	527	}
	528	}
	529	break;
	530	case BinOptItExMode.AVN:
	531	for (int s = 0; s < range.Item2; s++) {
	532	ap = (Int32)range.Item3 + range.Item1 + s arrStepInc;
	533	bp = (Int32*)range.Item4;
	534	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	535	for (int l = 0; l < dimLen; l++) {
	536	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	537	ap += arrInc;
	538	bp++;
	539	cp += arrInc;
	540	}
	541	}
	542	break;
	543	}
	544	System.Threading.Interlocked.Decrement(ref workerCount);
	545	};
	546	#endregion
	547
	548	#region work distribution
	549	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	550	int outLen = outDims.NumberOfElements;
	551	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	552	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	553	workItemLength = outLen / dimLen / workItemCount;
	554	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	555	} else {
	556	workItemLength = outLen / dimLen / 2;
	557	workItemCount = 2;
	558	}
	559	} else {
	560	workItemLength = outLen / dimLen;
	561	workItemCount = 1;
	562	}
	563
	564	fixed ( Int32* arrAP = arrA)
	565	fixed ( Int32* arrBP = arrB)
	566	fixed (byte* retArrP = retArr) {
	567
	568	for (; i < workItemCount - 1; i++) {
	569	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	570	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	571	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	572	System.Threading.Interlocked.Increment(ref workerCount);
	573	ILThreadPool.QueueUserWorkItem(i, worker, range);
	574	}
	575	// the last (or may the only) chunk is done right here
	576	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	577	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	578	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	579
	580	ILThreadPool.Wait4Workers(ref workerCount);
	581	}
	582	#endregion
	583
	584	return new ILRetLogical(retStorage);
	585	}
	586	}
	587
	588	/// <summary>Elementwise logical 'or' operator</summary>
	589	/// <param name="A">Input array A</param>
	590	/// <param name="B">Input array B</param>
	591	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	592	/// <remarks><para>On empty input an empty array will be returned.</para>
	593	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	594	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	595	public unsafe static ILRetLogical or(ILInArray<float> A, ILInArray<float> B) {
	596	using (ILScope.Enter(A, B)) {
	597	int outLen;
	598	BinOpItMode mode;
	599	byte[] retArr;
	600
	601	float[] arrA = A.GetArrayForRead();
	602
	603	float[] arrB = B.GetArrayForRead();
	604	ILSize outDims;
	605	if (A.IsScalar) {
	606	if (B.IsScalar) {
	607	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	608	} else if (B.IsEmpty) {
	609	return new ILRetLogical(B.Size);
	610	} else {
	611	outLen = B.S.NumberOfElements;
	612	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	613	mode = BinOpItMode.SAN;
	614	}
	615	outDims = B.Size;
	616	} else {
	617	outDims = A.Size;
	618	if (B.IsScalar) {
	619	if (A.IsEmpty) {
	620	return new ILRetLogical(A.Size);
	621	}
	622	outLen = A.S.NumberOfElements;
	623	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	624	mode = BinOpItMode.ASN;
	625	} else {
	626	// array + array
	627	if (!A.Size.IsSameSize(B.Size)) {
	628	return orEx(A,B);
	629	}
	630	outLen = A.S.NumberOfElements;
	631	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	632	mode = BinOpItMode.AAN;
	633	}
	634	}
	635	int workerCount = 1;
	636	Action<object> worker = data => {
	637	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	638	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	639	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	640	float scalar;
	641	cLast = cp + range.Item2;
	642	#region loops
	643	switch (mode) {
	644	case BinOpItMode.AAN:
	645	float* ap = ((float*)range.Item3 + range.Item1);
	646	float* bp = ((float*)range.Item4 + range.Item1);
	647	while (cp < cLast) {
	648	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	649	cp++;
	650	ap++;
	651	bp++;
	652	}
	653	break;
	654	case BinOpItMode.ASN:
	655	ap = ((float*)range.Item3 + range.Item1);
	656	scalar = ((float)range.Item4);
	657	while (cp < cLast) {
	658	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	659	cp++;
	660	ap++;
	661
	662	}
	663	break;
	664	case BinOpItMode.SAN:
	665	scalar = ((float)range.Item3);
	666	bp = ((float*)range.Item4 + range.Item1);
	667	while (cp < cLast) {
	668	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	669	cp++;
	670	bp++;
	671	}
	672	break;
	673	default:
	674	break;
	675	}
	676	#endregion
	677	System.Threading.Interlocked.Decrement(ref workerCount);
	678	};
	679
	680	#region do the work
	681	fixed (float* arrAP = arrA)
	682	fixed (float* arrBP = arrB)
	683	fixed (byte* retArrP = retArr) {
	684	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	685	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	686	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	687	workItemLength = outLen / workItemCount;
	688	} else {
	689	workItemLength = outLen / 2;
	690	workItemCount = 2;
	691	}
	692	} else {
	693	workItemLength = outLen;
	694	workItemCount = 1;
	695	}
	696
	697	for (; i < workItemCount - 1; i++) {
	698	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	699	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	700	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	701	System.Threading.Interlocked.Increment(ref workerCount);
	702	ILThreadPool.QueueUserWorkItem(i, worker, range);
	703	}
	704	// the last (or may the only) chunk is done right here
	705	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	706	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	707	ILThreadPool.Wait4Workers(ref workerCount);
	708	}
	709
	710	#endregion
	711	return new ILRetLogical(retArr, outDims);
	712	}
	713	}
	714
	715	private static unsafe ILRetLogical orEx(ILInArray<float> A, ILInArray<float> B) {
	716	using (ILScope.Enter(A, B)) {
	717	#region parameter checking
	718	if (isnull(A) \|\| isnull(B))
	719	return new ILRetLogical(ILSize.Empty00);
	720	if (A.IsEmpty) {
	721	return new ILRetLogical(B.S);
	722	} else if (B.IsEmpty) {
	723	return new ILRetLogical(A.S);
	724	}
	725	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	726	// return add(A,B);
	727	int dim = -1;
	728	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	729	if (A.S[l] != B.S[l]) {
	730	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	731	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	732	}
	733	dim = l;
	734	}
	735	}
	736	if (dim > 1)
	737	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	738	#endregion
	739
	740	#region parameter preparation
	741	byte[] retArr;
	742
	743	float[] arrA = A.GetArrayForRead();
	744
	745	float[] arrB = B.GetArrayForRead();
	746	ILSize outDims;
	747	BinOptItExMode mode;
	748	int arrInc = 0;
	749	int arrStepInc = 0;
	750	int dimLen = 0;
	751	if (A.IsVector) {
	752	outDims = B.S;
	753	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	754	mode = BinOptItExMode.VAN;
	755	dimLen = A.Length;
	756	} else if (B.IsVector) {
	757	outDims = A.S;
	758	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	759	mode = BinOptItExMode.AVN;
	760	dimLen = B.Length;
	761	} else {
	762	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	763	}
	764	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	765	arrStepInc = outDims.SequentialIndexDistance(dim);
	766	#endregion
	767
	768	#region worker loops definition
	769	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	770	int workerCount = 1;
	771	Action<object> worker = data => {
	772	// expects: iStart, iLen, ap, bp, cp
	773	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	774	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	775
	776	float* ap;
	777
	778	float* bp;
	779	byte* cp;
	780	switch (mode) {
	781	case BinOptItExMode.VAN:
	782	for (int s = 0; s < range.Item2; s++) {
	783	ap = (float*)range.Item3;
	784	bp = (float)range.Item4 + range.Item1 + s arrStepInc; ;
	785	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	786	for (int l = 0; l < dimLen; l++) {
	787	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	788	ap++;
	789	bp += arrInc;
	790	cp += arrInc;
	791	}
	792	}
	793	break;
	794	case BinOptItExMode.AVN:
	795	for (int s = 0; s < range.Item2; s++) {
	796	ap = (float)range.Item3 + range.Item1 + s arrStepInc;
	797	bp = (float*)range.Item4;
	798	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	799	for (int l = 0; l < dimLen; l++) {
	800	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	801	ap += arrInc;
	802	bp++;
	803	cp += arrInc;
	804	}
	805	}
	806	break;
	807	}
	808	System.Threading.Interlocked.Decrement(ref workerCount);
	809	};
	810	#endregion
	811
	812	#region work distribution
	813	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	814	int outLen = outDims.NumberOfElements;
	815	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	816	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	817	workItemLength = outLen / dimLen / workItemCount;
	818	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	819	} else {
	820	workItemLength = outLen / dimLen / 2;
	821	workItemCount = 2;
	822	}
	823	} else {
	824	workItemLength = outLen / dimLen;
	825	workItemCount = 1;
	826	}
	827
	828	fixed ( float* arrAP = arrA)
	829	fixed ( float* arrBP = arrB)
	830	fixed (byte* retArrP = retArr) {
	831
	832	for (; i < workItemCount - 1; i++) {
	833	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	834	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	835	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	836	System.Threading.Interlocked.Increment(ref workerCount);
	837	ILThreadPool.QueueUserWorkItem(i, worker, range);
	838	}
	839	// the last (or may the only) chunk is done right here
	840	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	841	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	842	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	843
	844	ILThreadPool.Wait4Workers(ref workerCount);
	845	}
	846	#endregion
	847
	848	return new ILRetLogical(retStorage);
	849	}
	850	}
	851
	852	/// <summary>Elementwise logical 'or' operator</summary>
	853	/// <param name="A">Input array A</param>
	854	/// <param name="B">Input array B</param>
	855	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	856	/// <remarks><para>On empty input an empty array will be returned.</para>
	857	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	858	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	859	public unsafe static ILRetLogical or(ILInArray<fcomplex> A, ILInArray<fcomplex> B) {
	860	using (ILScope.Enter(A, B)) {
	861	int outLen;
	862	BinOpItMode mode;
	863	byte[] retArr;
	864
	865	fcomplex[] arrA = A.GetArrayForRead();
	866
	867	fcomplex[] arrB = B.GetArrayForRead();
	868	ILSize outDims;
	869	if (A.IsScalar) {
	870	if (B.IsScalar) {
	871	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	872	} else if (B.IsEmpty) {
	873	return new ILRetLogical(B.Size);
	874	} else {
	875	outLen = B.S.NumberOfElements;
	876	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	877	mode = BinOpItMode.SAN;
	878	}
	879	outDims = B.Size;
	880	} else {
	881	outDims = A.Size;
	882	if (B.IsScalar) {
	883	if (A.IsEmpty) {
	884	return new ILRetLogical(A.Size);
	885	}
	886	outLen = A.S.NumberOfElements;
	887	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	888	mode = BinOpItMode.ASN;
	889	} else {
	890	// array + array
	891	if (!A.Size.IsSameSize(B.Size)) {
	892	return orEx(A,B);
	893	}
	894	outLen = A.S.NumberOfElements;
	895	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	896	mode = BinOpItMode.AAN;
	897	}
	898	}
	899	int workerCount = 1;
	900	Action<object> worker = data => {
	901	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	902	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	903	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	904	fcomplex scalar;
	905	cLast = cp + range.Item2;
	906	#region loops
	907	switch (mode) {
	908	case BinOpItMode.AAN:
	909	fcomplex* ap = ((fcomplex*)range.Item3 + range.Item1);
	910	fcomplex* bp = ((fcomplex*)range.Item4 + range.Item1);
	911	while (cp < cLast) {
	912	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	913	cp++;
	914	ap++;
	915	bp++;
	916	}
	917	break;
	918	case BinOpItMode.ASN:
	919	ap = ((fcomplex*)range.Item3 + range.Item1);
	920	scalar = ((fcomplex)range.Item4);
	921	while (cp < cLast) {
	922	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	923	cp++;
	924	ap++;
	925
	926	}
	927	break;
	928	case BinOpItMode.SAN:
	929	scalar = ((fcomplex)range.Item3);
	930	bp = ((fcomplex*)range.Item4 + range.Item1);
	931	while (cp < cLast) {
	932	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	933	cp++;
	934	bp++;
	935	}
	936	break;
	937	default:
	938	break;
	939	}
	940	#endregion
	941	System.Threading.Interlocked.Decrement(ref workerCount);
	942	};
	943
	944	#region do the work
	945	fixed (fcomplex* arrAP = arrA)
	946	fixed (fcomplex* arrBP = arrB)
	947	fixed (byte* retArrP = retArr) {
	948	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	949	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	950	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	951	workItemLength = outLen / workItemCount;
	952	} else {
	953	workItemLength = outLen / 2;
	954	workItemCount = 2;
	955	}
	956	} else {
	957	workItemLength = outLen;
	958	workItemCount = 1;
	959	}
	960
	961	for (; i < workItemCount - 1; i++) {
	962	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	963	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	964	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	965	System.Threading.Interlocked.Increment(ref workerCount);
	966	ILThreadPool.QueueUserWorkItem(i, worker, range);
	967	}
	968	// the last (or may the only) chunk is done right here
	969	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	970	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	971	ILThreadPool.Wait4Workers(ref workerCount);
	972	}
	973
	974	#endregion
	975	return new ILRetLogical(retArr, outDims);
	976	}
	977	}
	978
	979	private static unsafe ILRetLogical orEx(ILInArray<fcomplex> A, ILInArray<fcomplex> B) {
	980	using (ILScope.Enter(A, B)) {
	981	#region parameter checking
	982	if (isnull(A) \|\| isnull(B))
	983	return new ILRetLogical(ILSize.Empty00);
	984	if (A.IsEmpty) {
	985	return new ILRetLogical(B.S);
	986	} else if (B.IsEmpty) {
	987	return new ILRetLogical(A.S);
	988	}
	989	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	990	// return add(A,B);
	991	int dim = -1;
	992	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	993	if (A.S[l] != B.S[l]) {
	994	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	995	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	996	}
	997	dim = l;
	998	}
	999	}
	1000	if (dim > 1)
	1001	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	1002	#endregion
	1003
	1004	#region parameter preparation
	1005	byte[] retArr;
	1006
	1007	fcomplex[] arrA = A.GetArrayForRead();
	1008
	1009	fcomplex[] arrB = B.GetArrayForRead();
	1010	ILSize outDims;
	1011	BinOptItExMode mode;
	1012	int arrInc = 0;
	1013	int arrStepInc = 0;
	1014	int dimLen = 0;
	1015	if (A.IsVector) {
	1016	outDims = B.S;
	1017	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1018	mode = BinOptItExMode.VAN;
	1019	dimLen = A.Length;
	1020	} else if (B.IsVector) {
	1021	outDims = A.S;
	1022	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1023	mode = BinOptItExMode.AVN;
	1024	dimLen = B.Length;
	1025	} else {
	1026	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1027	}
	1028	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	1029	arrStepInc = outDims.SequentialIndexDistance(dim);
	1030	#endregion
	1031
	1032	#region worker loops definition
	1033	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	1034	int workerCount = 1;
	1035	Action<object> worker = data => {
	1036	// expects: iStart, iLen, ap, bp, cp
	1037	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	1038	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	1039
	1040	fcomplex* ap;
	1041
	1042	fcomplex* bp;
	1043	byte* cp;
	1044	switch (mode) {
	1045	case BinOptItExMode.VAN:
	1046	for (int s = 0; s < range.Item2; s++) {
	1047	ap = (fcomplex*)range.Item3;
	1048	bp = (fcomplex)range.Item4 + range.Item1 + s arrStepInc; ;
	1049	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1050	for (int l = 0; l < dimLen; l++) {
	1051	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1052	ap++;
	1053	bp += arrInc;
	1054	cp += arrInc;
	1055	}
	1056	}
	1057	break;
	1058	case BinOptItExMode.AVN:
	1059	for (int s = 0; s < range.Item2; s++) {
	1060	ap = (fcomplex)range.Item3 + range.Item1 + s arrStepInc;
	1061	bp = (fcomplex*)range.Item4;
	1062	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1063	for (int l = 0; l < dimLen; l++) {
	1064	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1065	ap += arrInc;
	1066	bp++;
	1067	cp += arrInc;
	1068	}
	1069	}
	1070	break;
	1071	}
	1072	System.Threading.Interlocked.Decrement(ref workerCount);
	1073	};
	1074	#endregion
	1075
	1076	#region work distribution
	1077	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1078	int outLen = outDims.NumberOfElements;
	1079	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	1080	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1081	workItemLength = outLen / dimLen / workItemCount;
	1082	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	1083	} else {
	1084	workItemLength = outLen / dimLen / 2;
	1085	workItemCount = 2;
	1086	}
	1087	} else {
	1088	workItemLength = outLen / dimLen;
	1089	workItemCount = 1;
	1090	}
	1091
	1092	fixed ( fcomplex* arrAP = arrA)
	1093	fixed ( fcomplex* arrBP = arrB)
	1094	fixed (byte* retArrP = retArr) {
	1095
	1096	for (; i < workItemCount - 1; i++) {
	1097	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	1098	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1099	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	1100	System.Threading.Interlocked.Increment(ref workerCount);
	1101	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1102	}
	1103	// the last (or may the only) chunk is done right here
	1104	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	1105	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1106	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	1107
	1108	ILThreadPool.Wait4Workers(ref workerCount);
	1109	}
	1110	#endregion
	1111
	1112	return new ILRetLogical(retStorage);
	1113	}
	1114	}
	1115
	1116	/// <summary>Elementwise logical 'or' operator</summary>
	1117	/// <param name="A">Input array A</param>
	1118	/// <param name="B">Input array B</param>
	1119	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	1120	/// <remarks><para>On empty input an empty array will be returned.</para>
	1121	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	1122	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	1123	public unsafe static ILRetLogical or(ILInArray<complex> A, ILInArray<complex> B) {
	1124	using (ILScope.Enter(A, B)) {
	1125	int outLen;
	1126	BinOpItMode mode;
	1127	byte[] retArr;
	1128
	1129	complex[] arrA = A.GetArrayForRead();
	1130
	1131	complex[] arrB = B.GetArrayForRead();
	1132	ILSize outDims;
	1133	if (A.IsScalar) {
	1134	if (B.IsScalar) {
	1135	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	1136	} else if (B.IsEmpty) {
	1137	return new ILRetLogical(B.Size);
	1138	} else {
	1139	outLen = B.S.NumberOfElements;
	1140	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1141	mode = BinOpItMode.SAN;
	1142	}
	1143	outDims = B.Size;
	1144	} else {
	1145	outDims = A.Size;
	1146	if (B.IsScalar) {
	1147	if (A.IsEmpty) {
	1148	return new ILRetLogical(A.Size);
	1149	}
	1150	outLen = A.S.NumberOfElements;
	1151	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1152	mode = BinOpItMode.ASN;
	1153	} else {
	1154	// array + array
	1155	if (!A.Size.IsSameSize(B.Size)) {
	1156	return orEx(A,B);
	1157	}
	1158	outLen = A.S.NumberOfElements;
	1159	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1160	mode = BinOpItMode.AAN;
	1161	}
	1162	}
	1163	int workerCount = 1;
	1164	Action<object> worker = data => {
	1165	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	1166	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	1167	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	1168	complex scalar;
	1169	cLast = cp + range.Item2;
	1170	#region loops
	1171	switch (mode) {
	1172	case BinOpItMode.AAN:
	1173	complex* ap = ((complex*)range.Item3 + range.Item1);
	1174	complex* bp = ((complex*)range.Item4 + range.Item1);
	1175	while (cp < cLast) {
	1176	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1177	cp++;
	1178	ap++;
	1179	bp++;
	1180	}
	1181	break;
	1182	case BinOpItMode.ASN:
	1183	ap = ((complex*)range.Item3 + range.Item1);
	1184	scalar = ((complex)range.Item4);
	1185	while (cp < cLast) {
	1186	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	1187	cp++;
	1188	ap++;
	1189
	1190	}
	1191	break;
	1192	case BinOpItMode.SAN:
	1193	scalar = ((complex)range.Item3);
	1194	bp = ((complex*)range.Item4 + range.Item1);
	1195	while (cp < cLast) {
	1196	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	1197	cp++;
	1198	bp++;
	1199	}
	1200	break;
	1201	default:
	1202	break;
	1203	}
	1204	#endregion
	1205	System.Threading.Interlocked.Decrement(ref workerCount);
	1206	};
	1207
	1208	#region do the work
	1209	fixed (complex* arrAP = arrA)
	1210	fixed (complex* arrBP = arrB)
	1211	fixed (byte* retArrP = retArr) {
	1212	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1213	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	1214	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1215	workItemLength = outLen / workItemCount;
	1216	} else {
	1217	workItemLength = outLen / 2;
	1218	workItemCount = 2;
	1219	}
	1220	} else {
	1221	workItemLength = outLen;
	1222	workItemCount = 1;
	1223	}
	1224
	1225	for (; i < workItemCount - 1; i++) {
	1226	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	1227	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	1228	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	1229	System.Threading.Interlocked.Increment(ref workerCount);
	1230	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1231	}
	1232	// the last (or may the only) chunk is done right here
	1233	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	1234	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	1235	ILThreadPool.Wait4Workers(ref workerCount);
	1236	}
	1237
	1238	#endregion
	1239	return new ILRetLogical(retArr, outDims);
	1240	}
	1241	}
	1242
	1243	private static unsafe ILRetLogical orEx(ILInArray<complex> A, ILInArray<complex> B) {
	1244	using (ILScope.Enter(A, B)) {
	1245	#region parameter checking
	1246	if (isnull(A) \|\| isnull(B))
	1247	return new ILRetLogical(ILSize.Empty00);
	1248	if (A.IsEmpty) {
	1249	return new ILRetLogical(B.S);
	1250	} else if (B.IsEmpty) {
	1251	return new ILRetLogical(A.S);
	1252	}
	1253	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	1254	// return add(A,B);
	1255	int dim = -1;
	1256	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	1257	if (A.S[l] != B.S[l]) {
	1258	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	1259	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1260	}
	1261	dim = l;
	1262	}
	1263	}
	1264	if (dim > 1)
	1265	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	1266	#endregion
	1267
	1268	#region parameter preparation
	1269	byte[] retArr;
	1270
	1271	complex[] arrA = A.GetArrayForRead();
	1272
	1273	complex[] arrB = B.GetArrayForRead();
	1274	ILSize outDims;
	1275	BinOptItExMode mode;
	1276	int arrInc = 0;
	1277	int arrStepInc = 0;
	1278	int dimLen = 0;
	1279	if (A.IsVector) {
	1280	outDims = B.S;
	1281	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1282	mode = BinOptItExMode.VAN;
	1283	dimLen = A.Length;
	1284	} else if (B.IsVector) {
	1285	outDims = A.S;
	1286	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1287	mode = BinOptItExMode.AVN;
	1288	dimLen = B.Length;
	1289	} else {
	1290	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1291	}
	1292	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	1293	arrStepInc = outDims.SequentialIndexDistance(dim);
	1294	#endregion
	1295
	1296	#region worker loops definition
	1297	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	1298	int workerCount = 1;
	1299	Action<object> worker = data => {
	1300	// expects: iStart, iLen, ap, bp, cp
	1301	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	1302	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	1303
	1304	complex* ap;
	1305
	1306	complex* bp;
	1307	byte* cp;
	1308	switch (mode) {
	1309	case BinOptItExMode.VAN:
	1310	for (int s = 0; s < range.Item2; s++) {
	1311	ap = (complex*)range.Item3;
	1312	bp = (complex)range.Item4 + range.Item1 + s arrStepInc; ;
	1313	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1314	for (int l = 0; l < dimLen; l++) {
	1315	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1316	ap++;
	1317	bp += arrInc;
	1318	cp += arrInc;
	1319	}
	1320	}
	1321	break;
	1322	case BinOptItExMode.AVN:
	1323	for (int s = 0; s < range.Item2; s++) {
	1324	ap = (complex)range.Item3 + range.Item1 + s arrStepInc;
	1325	bp = (complex*)range.Item4;
	1326	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1327	for (int l = 0; l < dimLen; l++) {
	1328	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1329	ap += arrInc;
	1330	bp++;
	1331	cp += arrInc;
	1332	}
	1333	}
	1334	break;
	1335	}
	1336	System.Threading.Interlocked.Decrement(ref workerCount);
	1337	};
	1338	#endregion
	1339
	1340	#region work distribution
	1341	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1342	int outLen = outDims.NumberOfElements;
	1343	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	1344	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1345	workItemLength = outLen / dimLen / workItemCount;
	1346	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	1347	} else {
	1348	workItemLength = outLen / dimLen / 2;
	1349	workItemCount = 2;
	1350	}
	1351	} else {
	1352	workItemLength = outLen / dimLen;
	1353	workItemCount = 1;
	1354	}
	1355
	1356	fixed ( complex* arrAP = arrA)
	1357	fixed ( complex* arrBP = arrB)
	1358	fixed (byte* retArrP = retArr) {
	1359
	1360	for (; i < workItemCount - 1; i++) {
	1361	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	1362	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1363	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	1364	System.Threading.Interlocked.Increment(ref workerCount);
	1365	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1366	}
	1367	// the last (or may the only) chunk is done right here
	1368	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	1369	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1370	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	1371
	1372	ILThreadPool.Wait4Workers(ref workerCount);
	1373	}
	1374	#endregion
	1375
	1376	return new ILRetLogical(retStorage);
	1377	}
	1378	}
	1379
	1380	/// <summary>Elementwise logical 'or' operator</summary>
	1381	/// <param name="A">Input array A</param>
	1382	/// <param name="B">Input array B</param>
	1383	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	1384	/// <remarks><para>On empty input an empty array will be returned.</para>
	1385	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	1386	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	1387	public unsafe static ILRetLogical or(ILInArray<byte> A, ILInArray<byte> B) {
	1388	using (ILScope.Enter(A, B)) {
	1389	int outLen;
	1390	BinOpItMode mode;
	1391	byte[] retArr;
	1392
	1393	byte[] arrA = A.GetArrayForRead();
	1394
	1395	byte[] arrB = B.GetArrayForRead();
	1396	ILSize outDims;
	1397	if (A.IsScalar) {
	1398	if (B.IsScalar) {
	1399	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	1400	} else if (B.IsEmpty) {
	1401	return new ILRetLogical(B.Size);
	1402	} else {
	1403	outLen = B.S.NumberOfElements;
	1404	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1405	mode = BinOpItMode.SAN;
	1406	}
	1407	outDims = B.Size;
	1408	} else {
	1409	outDims = A.Size;
	1410	if (B.IsScalar) {
	1411	if (A.IsEmpty) {
	1412	return new ILRetLogical(A.Size);
	1413	}
	1414	outLen = A.S.NumberOfElements;
	1415	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1416	mode = BinOpItMode.ASN;
	1417	} else {
	1418	// array + array
	1419	if (!A.Size.IsSameSize(B.Size)) {
	1420	return orEx(A,B);
	1421	}
	1422	outLen = A.S.NumberOfElements;
	1423	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1424	mode = BinOpItMode.AAN;
	1425	}
	1426	}
	1427	int workerCount = 1;
	1428	Action<object> worker = data => {
	1429	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	1430	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	1431	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	1432	byte scalar;
	1433	cLast = cp + range.Item2;
	1434	#region loops
	1435	switch (mode) {
	1436	case BinOpItMode.AAN:
	1437	byte* ap = ((byte*)range.Item3 + range.Item1);
	1438	byte* bp = ((byte*)range.Item4 + range.Item1);
	1439	while (cp < cLast) {
	1440	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1441	cp++;
	1442	ap++;
	1443	bp++;
	1444	}
	1445	break;
	1446	case BinOpItMode.ASN:
	1447	ap = ((byte*)range.Item3 + range.Item1);
	1448	scalar = ((byte)range.Item4);
	1449	while (cp < cLast) {
	1450	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	1451	cp++;
	1452	ap++;
	1453
	1454	}
	1455	break;
	1456	case BinOpItMode.SAN:
	1457	scalar = ((byte)range.Item3);
	1458	bp = ((byte*)range.Item4 + range.Item1);
	1459	while (cp < cLast) {
	1460	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	1461	cp++;
	1462	bp++;
	1463	}
	1464	break;
	1465	default:
	1466	break;
	1467	}
	1468	#endregion
	1469	System.Threading.Interlocked.Decrement(ref workerCount);
	1470	};
	1471
	1472	#region do the work
	1473	fixed (byte* arrAP = arrA)
	1474	fixed (byte* arrBP = arrB)
	1475	fixed (byte* retArrP = retArr) {
	1476	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1477	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	1478	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1479	workItemLength = outLen / workItemCount;
	1480	} else {
	1481	workItemLength = outLen / 2;
	1482	workItemCount = 2;
	1483	}
	1484	} else {
	1485	workItemLength = outLen;
	1486	workItemCount = 1;
	1487	}
	1488
	1489	for (; i < workItemCount - 1; i++) {
	1490	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	1491	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	1492	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	1493	System.Threading.Interlocked.Increment(ref workerCount);
	1494	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1495	}
	1496	// the last (or may the only) chunk is done right here
	1497	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	1498	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	1499	ILThreadPool.Wait4Workers(ref workerCount);
	1500	}
	1501
	1502	#endregion
	1503	return new ILRetLogical(retArr, outDims);
	1504	}
	1505	}
	1506
	1507	private static unsafe ILRetLogical orEx(ILInArray<byte> A, ILInArray<byte> B) {
	1508	using (ILScope.Enter(A, B)) {
	1509	#region parameter checking
	1510	if (isnull(A) \|\| isnull(B))
	1511	return new ILRetLogical(ILSize.Empty00);
	1512	if (A.IsEmpty) {
	1513	return new ILRetLogical(B.S);
	1514	} else if (B.IsEmpty) {
	1515	return new ILRetLogical(A.S);
	1516	}
	1517	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	1518	// return add(A,B);
	1519	int dim = -1;
	1520	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	1521	if (A.S[l] != B.S[l]) {
	1522	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	1523	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1524	}
	1525	dim = l;
	1526	}
	1527	}
	1528	if (dim > 1)
	1529	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	1530	#endregion
	1531
	1532	#region parameter preparation
	1533	byte[] retArr;
	1534
	1535	byte[] arrA = A.GetArrayForRead();
	1536
	1537	byte[] arrB = B.GetArrayForRead();
	1538	ILSize outDims;
	1539	BinOptItExMode mode;
	1540	int arrInc = 0;
	1541	int arrStepInc = 0;
	1542	int dimLen = 0;
	1543	if (A.IsVector) {
	1544	outDims = B.S;
	1545	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1546	mode = BinOptItExMode.VAN;
	1547	dimLen = A.Length;
	1548	} else if (B.IsVector) {
	1549	outDims = A.S;
	1550	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1551	mode = BinOptItExMode.AVN;
	1552	dimLen = B.Length;
	1553	} else {
	1554	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1555	}
	1556	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	1557	arrStepInc = outDims.SequentialIndexDistance(dim);
	1558	#endregion
	1559
	1560	#region worker loops definition
	1561	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	1562	int workerCount = 1;
	1563	Action<object> worker = data => {
	1564	// expects: iStart, iLen, ap, bp, cp
	1565	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	1566	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	1567
	1568	byte* ap;
	1569
	1570	byte* bp;
	1571	byte* cp;
	1572	switch (mode) {
	1573	case BinOptItExMode.VAN:
	1574	for (int s = 0; s < range.Item2; s++) {
	1575	ap = (byte*)range.Item3;
	1576	bp = (byte)range.Item4 + range.Item1 + s arrStepInc; ;
	1577	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1578	for (int l = 0; l < dimLen; l++) {
	1579	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1580	ap++;
	1581	bp += arrInc;
	1582	cp += arrInc;
	1583	}
	1584	}
	1585	break;
	1586	case BinOptItExMode.AVN:
	1587	for (int s = 0; s < range.Item2; s++) {
	1588	ap = (byte)range.Item3 + range.Item1 + s arrStepInc;
	1589	bp = (byte*)range.Item4;
	1590	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1591	for (int l = 0; l < dimLen; l++) {
	1592	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1593	ap += arrInc;
	1594	bp++;
	1595	cp += arrInc;
	1596	}
	1597	}
	1598	break;
	1599	}
	1600	System.Threading.Interlocked.Decrement(ref workerCount);
	1601	};
	1602	#endregion
	1603
	1604	#region work distribution
	1605	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1606	int outLen = outDims.NumberOfElements;
	1607	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	1608	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1609	workItemLength = outLen / dimLen / workItemCount;
	1610	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	1611	} else {
	1612	workItemLength = outLen / dimLen / 2;
	1613	workItemCount = 2;
	1614	}
	1615	} else {
	1616	workItemLength = outLen / dimLen;
	1617	workItemCount = 1;
	1618	}
	1619
	1620	fixed ( byte* arrAP = arrA)
	1621	fixed ( byte* arrBP = arrB)
	1622	fixed (byte* retArrP = retArr) {
	1623
	1624	for (; i < workItemCount - 1; i++) {
	1625	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	1626	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1627	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	1628	System.Threading.Interlocked.Increment(ref workerCount);
	1629	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1630	}
	1631	// the last (or may the only) chunk is done right here
	1632	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	1633	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1634	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	1635
	1636	ILThreadPool.Wait4Workers(ref workerCount);
	1637	}
	1638	#endregion
	1639
	1640	return new ILRetLogical(retStorage);
	1641	}
	1642	}
	1643
	1644	/// <summary>Elementwise logical 'or' operator</summary>
	1645	/// <param name="A">Input array A</param>
	1646	/// <param name="B">Input array B</param>
	1647	/// <returns>Logical array with the elementwise result of logical 'or' for all elements in A and B</returns>
	1648	/// <remarks><para>On empty input an empty array will be returned.</para>
	1649	/// <para>A and/or B may be scalar. The scalar value will be applied on all elements of the other array.</para></remarks>
	1650	/// <exception cref="ILNumerics.Exceptions.ILDimensionMismatchException">If neither A nor B is scalar or empty, the dimensions of both arrays must match.</exception>
	1651	public unsafe static ILRetLogical or(ILInArray<double> A, ILInArray<double> B) {
	1652	using (ILScope.Enter(A, B)) {
	1653	int outLen;
	1654	BinOpItMode mode;
	1655	byte[] retArr;
	1656
	1657	double[] arrA = A.GetArrayForRead();
	1658
	1659	double[] arrB = B.GetArrayForRead();
	1660	ILSize outDims;
	1661	if (A.IsScalar) {
	1662	if (B.IsScalar) {
	1663	return new ILRetLogical(new byte[1] { (A.GetValue(0) != 0 \|\| B.GetValue(0) != 0) ? (byte)1 : (byte)0 });
	1664	} else if (B.IsEmpty) {
	1665	return new ILRetLogical(B.Size);
	1666	} else {
	1667	outLen = B.S.NumberOfElements;
	1668	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1669	mode = BinOpItMode.SAN;
	1670	}
	1671	outDims = B.Size;
	1672	} else {
	1673	outDims = A.Size;
	1674	if (B.IsScalar) {
	1675	if (A.IsEmpty) {
	1676	return new ILRetLogical(A.Size);
	1677	}
	1678	outLen = A.S.NumberOfElements;
	1679	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1680	mode = BinOpItMode.ASN;
	1681	} else {
	1682	// array + array
	1683	if (!A.Size.IsSameSize(B.Size)) {
	1684	return orEx(A,B);
	1685	}
	1686	outLen = A.S.NumberOfElements;
	1687	retArr = ILMemoryPool.Pool.New<byte>(outLen);
	1688	mode = BinOpItMode.AAN;
	1689	}
	1690	}
	1691	int workerCount = 1;
	1692	Action<object> worker = data => {
	1693	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	1694	= (Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>)data;
	1695	byte* cLast, cp = (byte*)range.Item5 + range.Item1;
	1696	double scalar;
	1697	cLast = cp + range.Item2;
	1698	#region loops
	1699	switch (mode) {
	1700	case BinOpItMode.AAN:
	1701	double* ap = ((double*)range.Item3 + range.Item1);
	1702	double* bp = ((double*)range.Item4 + range.Item1);
	1703	while (cp < cLast) {
	1704	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1705	cp++;
	1706	ap++;
	1707	bp++;
	1708	}
	1709	break;
	1710	case BinOpItMode.ASN:
	1711	ap = ((double*)range.Item3 + range.Item1);
	1712	scalar = ((double)range.Item4);
	1713	while (cp < cLast) {
	1714	cp = (scalar != 0 \|\| ap != 0) ? (byte)1 : (byte)0;
	1715	cp++;
	1716	ap++;
	1717
	1718	}
	1719	break;
	1720	case BinOpItMode.SAN:
	1721	scalar = ((double)range.Item3);
	1722	bp = ((double*)range.Item4 + range.Item1);
	1723	while (cp < cLast) {
	1724	cp = (scalar != 0 \|\| bp != 0) ? (byte)1 : (byte)0;
	1725	cp++;
	1726	bp++;
	1727	}
	1728	break;
	1729	default:
	1730	break;
	1731	}
	1732	#endregion
	1733	System.Threading.Interlocked.Decrement(ref workerCount);
	1734	};
	1735
	1736	#region do the work
	1737	fixed (double* arrAP = arrA)
	1738	fixed (double* arrBP = arrB)
	1739	fixed (byte* retArrP = retArr) {
	1740	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1741	if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
	1742	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1743	workItemLength = outLen / workItemCount;
	1744	} else {
	1745	workItemLength = outLen / 2;
	1746	workItemCount = 2;
	1747	}
	1748	} else {
	1749	workItemLength = outLen;
	1750	workItemCount = 1;
	1751	}
	1752
	1753	for (; i < workItemCount - 1; i++) {
	1754	Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode> range
	1755	= new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	1756	(i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode);
	1757	System.Threading.Interlocked.Increment(ref workerCount);
	1758	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1759	}
	1760	// the last (or may the only) chunk is done right here
	1761	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr, BinOpItMode>
	1762	(i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP, mode));
	1763	ILThreadPool.Wait4Workers(ref workerCount);
	1764	}
	1765
	1766	#endregion
	1767	return new ILRetLogical(retArr, outDims);
	1768	}
	1769	}
	1770
	1771	private static unsafe ILRetLogical orEx(ILInArray<double> A, ILInArray<double> B) {
	1772	using (ILScope.Enter(A, B)) {
	1773	#region parameter checking
	1774	if (isnull(A) \|\| isnull(B))
	1775	return new ILRetLogical(ILSize.Empty00);
	1776	if (A.IsEmpty) {
	1777	return new ILRetLogical(B.S);
	1778	} else if (B.IsEmpty) {
	1779	return new ILRetLogical(A.S);
	1780	}
	1781	//if (A.IsScalar \|\| B.IsScalar \|\| A.D.IsSameSize(B.D))
	1782	// return add(A,B);
	1783	int dim = -1;
	1784	for (int l = 0; l < Math.Max(A.S.NumberOfDimensions, B.S.NumberOfDimensions); l++) {
	1785	if (A.S[l] != B.S[l]) {
	1786	if (dim >= 0 \|\| (A.S[l] != 1 && B.S[l] != 1)) {
	1787	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1788	}
	1789	dim = l;
	1790	}
	1791	}
	1792	if (dim > 1)
	1793	throw new ILArgumentException("singleton dimension expansion currently is only supported for colum- and row vectors");
	1794	#endregion
	1795
	1796	#region parameter preparation
	1797	byte[] retArr;
	1798
	1799	double[] arrA = A.GetArrayForRead();
	1800
	1801	double[] arrB = B.GetArrayForRead();
	1802	ILSize outDims;
	1803	BinOptItExMode mode;
	1804	int arrInc = 0;
	1805	int arrStepInc = 0;
	1806	int dimLen = 0;
	1807	if (A.IsVector) {
	1808	outDims = B.S;
	1809	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1810	mode = BinOptItExMode.VAN;
	1811	dimLen = A.Length;
	1812	} else if (B.IsVector) {
	1813	outDims = A.S;
	1814	retArr = ILMemoryPool.Pool.New<byte>(outDims.NumberOfElements);
	1815	mode = BinOptItExMode.AVN;
	1816	dimLen = B.Length;
	1817	} else {
	1818	throw new ILArgumentException("A and B must have the same size except for one simgleton dimension in A or B");
	1819	}
	1820	arrInc = (dim == 0) ? outDims.SequentialIndexDistance(1) : outDims.SequentialIndexDistance(0);
	1821	arrStepInc = outDims.SequentialIndexDistance(dim);
	1822	#endregion
	1823
	1824	#region worker loops definition
	1825	ILLogicalStorage retStorage = new ILLogicalStorage(retArr, outDims);
	1826	int workerCount = 1;
	1827	Action<object> worker = data => {
	1828	// expects: iStart, iLen, ap, bp, cp
	1829	Tuple<int, int, IntPtr, IntPtr, IntPtr> range =
	1830	(Tuple<int, int, IntPtr, IntPtr, IntPtr>)data;
	1831
	1832	double* ap;
	1833
	1834	double* bp;
	1835	byte* cp;
	1836	switch (mode) {
	1837	case BinOptItExMode.VAN:
	1838	for (int s = 0; s < range.Item2; s++) {
	1839	ap = (double*)range.Item3;
	1840	bp = (double)range.Item4 + range.Item1 + s arrStepInc; ;
	1841	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1842	for (int l = 0; l < dimLen; l++) {
	1843	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1844	ap++;
	1845	bp += arrInc;
	1846	cp += arrInc;
	1847	}
	1848	}
	1849	break;
	1850	case BinOptItExMode.AVN:
	1851	for (int s = 0; s < range.Item2; s++) {
	1852	ap = (double)range.Item3 + range.Item1 + s arrStepInc;
	1853	bp = (double*)range.Item4;
	1854	cp = (byte)range.Item5 + range.Item1 + s arrStepInc;
	1855	for (int l = 0; l < dimLen; l++) {
	1856	cp = (ap != 0 \|\| *bp != 0) ? (byte)1 : (byte)0;
	1857	ap += arrInc;
	1858	bp++;
	1859	cp += arrInc;
	1860	}
	1861	}
	1862	break;
	1863	}
	1864	System.Threading.Interlocked.Decrement(ref workerCount);
	1865	};
	1866	#endregion
	1867
	1868	#region work distribution
	1869	int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength;
	1870	int outLen = outDims.NumberOfElements;
	1871	if (Settings.s_maxNumberThreads > 1 && outLen / 2 >= Settings.s_minParallelElement1Count) {
	1872	if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
	1873	workItemLength = outLen / dimLen / workItemCount;
	1874	//workItemLength = (int)((double)outLen / workItemCount * 1.05);
	1875	} else {
	1876	workItemLength = outLen / dimLen / 2;
	1877	workItemCount = 2;
	1878	}
	1879	} else {
	1880	workItemLength = outLen / dimLen;
	1881	workItemCount = 1;
	1882	}
	1883
	1884	fixed ( double* arrAP = arrA)
	1885	fixed ( double* arrBP = arrB)
	1886	fixed (byte* retArrP = retArr) {
	1887
	1888	for (; i < workItemCount - 1; i++) {
	1889	Tuple<int, int, IntPtr, IntPtr, IntPtr> range
	1890	= new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1891	(i * workItemLength * arrStepInc, workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP);
	1892	System.Threading.Interlocked.Increment(ref workerCount);
	1893	ILThreadPool.QueueUserWorkItem(i, worker, range);
	1894	}
	1895	// the last (or may the only) chunk is done right here
	1896	//System.Threading.Interlocked.Increment(ref retStorage.PendingTasks);
	1897	worker(new Tuple<int, int, IntPtr, IntPtr, IntPtr>
	1898	(i * workItemLength * arrStepInc, (outLen / dimLen) - i * workItemLength, (IntPtr)arrAP, (IntPtr)arrBP, (IntPtr)retArrP));
	1899
	1900	ILThreadPool.Wait4Workers(ref workerCount);
	1901	}
	1902	#endregion
	1903
	1904	return new ILRetLogical(retStorage);
	1905	}
	1906	}
	1907
	1908
	1909	#endregion HYCALPER AUTO GENERATED CODE
	1910
	1911	}
	1912	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Update cookies preferences