///
///    This file is part of ILNumerics Community Edition.
///
///    ILNumerics Community Edition - high performance computing for applications.
///    Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net
///
///    ILNumerics Community Edition is free software: you can redistribute it and/or modify
///    it under the terms of the GNU General Public License version 3 as published by
///    the Free Software Foundation.
///
///    ILNumerics Community Edition is distributed in the hope that it will be useful,
///    but WITHOUT ANY WARRANTY; without even the implied warranty of
///    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
///    GNU General Public License for more details.
///
///    You should have received a copy of the GNU General Public License
///    along with ILNumerics Community Edition. See the file License.txt in the root
///    of your distribution package. If not, see <http://www.gnu.org/licenses/>.
///
///    In addition this software uses the following components and/or licenses: 
///
///    =================================================================================
///    The Open Toolkit Library License
///    
///    Copyright (c) 2006 - 2009 the Open Toolkit library.
///    
///    Permission is hereby granted, free of charge, to any person obtaining a copy
///    of this software and associated documentation files (the "Software"), to deal
///    in the Software without restriction, including without limitation the rights to 
///    use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
///    the Software, and to permit persons to whom the Software is furnished to do
///    so, subject to the following conditions:
///
///    The above copyright notice and this permission notice shall be included in all
///    copies or substantial portions of the Software.
///
///    =================================================================================
///    

using System;
using System.Collections.Generic;
using System.Text;
using ILNumerics;
using ILNumerics.Exceptions;
using ILNumerics.Storage;
using ILNumerics.Misc;


namespace ILNumerics {
    public partial class ILMath {



#region HYCALPER AUTO GENERATED CODE

        /// <summary>Square root of array elements - complex output</summary>
        /// <param name="A">Input array - positive and negative values allowed</param>
        /// <returns>Square root of elements of A - complex output</returns>
        /// <remarks><para>If the input array is empty, an empty array will be returned.</para>
        /// <para>The array returned will be a dense array.</para></remarks>
        public unsafe static  ILRetArray<fcomplex>  sqrtc (ILInArray< float > A) {
            using (ILScope.Enter(A)) {
                if (A.IsEmpty)
                    return new  ILRetArray<fcomplex>(A.Size); 
                ILSize inDim = A.Size;
                float[] arrA = A.GetArrayForRead(); 
                fcomplex [] retArr;
                int outLen = inDim.NumberOfElements;
                bool inplace = true;
               
                if (true) {
                    retArr = ILMemoryPool.Pool.New<fcomplex>(outLen);
                    inplace = false; 
                }
                int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength, workerCount = 1;
                if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
                    if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
                        workItemLength = outLen / workItemCount;
                        //workItemLength = (int)((double)outLen / workItemCount * 1.05);
                    } else {
                        workItemLength = outLen / 2;
                        workItemCount = 2;
                    }
                } else {
                    workItemLength = outLen;
                    workItemCount = 1;
                }
                ILDenseStorage<fcomplex> retStorage = new ILDenseStorage<fcomplex>(retArr, inDim);
                
                Action<object> worker = data => {
                    Tuple<int, int, IntPtr, IntPtr, bool> range = (Tuple<int, int, IntPtr, IntPtr, bool>)data;
                   
                    fcomplex* cp = ((fcomplex*)range.Item4 + range.Item1);
                    int len = range.Item2; 
                    if (range.Item5) { 
                        // inplace
                        while (len > 20) {
                            cp[0] =  fcomplex.Sqrt(cp[0]  )  /*dummy*/;
                            cp[1] =  fcomplex.Sqrt(cp[1]  )  /*dummy*/;
                            cp[2] =  fcomplex.Sqrt(cp[2]  )  /*dummy*/;
                            cp[3] =  fcomplex.Sqrt(cp[3]  )  /*dummy*/;
                            cp[4] =  fcomplex.Sqrt(cp[4]  )  /*dummy*/;
                            cp[5] =  fcomplex.Sqrt(cp[5]  )  /*dummy*/;
                            cp[6] =  fcomplex.Sqrt(cp[6]  )  /*dummy*/;
                            cp[7] =  fcomplex.Sqrt(cp[7]  )  /*dummy*/;
                            cp[8] =  fcomplex.Sqrt(cp[8]  )  /*dummy*/;
                            cp[9] =  fcomplex.Sqrt(cp[9]  )  /*dummy*/;
                            cp[10] =  fcomplex.Sqrt(cp[10]  )  /*dummy*/;
                            cp[11] =  fcomplex.Sqrt(cp[11]  )  /*dummy*/;
                            cp[12] =  fcomplex.Sqrt(cp[12]  )  /*dummy*/;
                            cp[13] =  fcomplex.Sqrt(cp[13]  )  /*dummy*/;
                            cp[14] =  fcomplex.Sqrt(cp[14]  )  /*dummy*/;
                            cp[15] =  fcomplex.Sqrt(cp[15]  )  /*dummy*/;
                            cp[16] =  fcomplex.Sqrt(cp[16]  )  /*dummy*/;
                            cp[17] =  fcomplex.Sqrt(cp[17]  )  /*dummy*/;
                            cp[18] =  fcomplex.Sqrt(cp[18]  )  /*dummy*/;
                            cp[19] =  fcomplex.Sqrt(cp[19]  )  /*dummy*/;
                            cp[20] =  fcomplex.Sqrt(cp[20]  )  /*dummy*/;
                            cp+=21; len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  fcomplex.Sqrt(*cp  )  /*dummy*/;
                            cp++;
                        }
                    } else {
                        float* ap = ((float*)range.Item3 + range.Item1);
                        while (len > 20) {
                            cp[0] =  fcomplex.Sqrt(ap[0]  )  /*dummy*/;
                            cp[1] =  fcomplex.Sqrt(ap[1]  )  /*dummy*/;
                            cp[2] =  fcomplex.Sqrt(ap[2]  )  /*dummy*/;
                            cp[3] =  fcomplex.Sqrt(ap[3]  )  /*dummy*/;
                            cp[4] =  fcomplex.Sqrt(ap[4]  )  /*dummy*/;
                            cp[5] =  fcomplex.Sqrt(ap[5]  )  /*dummy*/;
                            cp[6] =  fcomplex.Sqrt(ap[6]  )  /*dummy*/;
                            cp[7] =  fcomplex.Sqrt(ap[7]  )  /*dummy*/;
                            cp[8] =  fcomplex.Sqrt(ap[8]  )  /*dummy*/;
                            cp[9] =  fcomplex.Sqrt(ap[9]  )  /*dummy*/;
                            cp[10] =  fcomplex.Sqrt(ap[10]  )  /*dummy*/;
                            cp[11] =  fcomplex.Sqrt(ap[11]  )  /*dummy*/;
                            cp[12] =  fcomplex.Sqrt(ap[12]  )  /*dummy*/;
                            cp[13] =  fcomplex.Sqrt(ap[13]  )  /*dummy*/;
                            cp[14] =  fcomplex.Sqrt(ap[14]  )  /*dummy*/;
                            cp[15] =  fcomplex.Sqrt(ap[15]  )  /*dummy*/;
                            cp[16] =  fcomplex.Sqrt(ap[16]  )  /*dummy*/;
                            cp[17] =  fcomplex.Sqrt(ap[17]  )  /*dummy*/;
                            cp[18] =  fcomplex.Sqrt(ap[18]  )  /*dummy*/;
                            cp[19] =  fcomplex.Sqrt(ap[19]  )  /*dummy*/;
                            cp[20] =  fcomplex.Sqrt(ap[20]  )  /*dummy*/;
                            ap += 21;
                            cp += 21;
                            len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  fcomplex.Sqrt(*ap  )  /*dummy*/;
                            ap++;
                            cp++;
                        }
                    }
                    System.Threading.Interlocked.Decrement(ref workerCount);
                };

                fixed ( float* arrAP = arrA)
                fixed ( fcomplex* retArrP = retArr) {
                    for (; i < workItemCount - 1; i++) {
                        Tuple<int, int, IntPtr, IntPtr, bool> range
                            = new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace);
                        System.Threading.Interlocked.Increment(ref workerCount);
                        ILThreadPool.QueueUserWorkItem(i,worker, range);
                    }
                    // the last (or may the only) chunk is done right here
                    worker(new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace));

                    ILThreadPool.Wait4Workers(ref workerCount);
                }
                return new  ILRetArray<fcomplex>(retStorage);
            }
        }
        /// <summary>Square root of array elements - complex output</summary>
        /// <param name="A">Input array - positive and negative values allowed</param>
        /// <returns>Square root of elements of A - complex output</returns>
        /// <remarks><para>If the input array is empty, an empty array will be returned.</para>
        /// <para>The array returned will be a dense array.</para></remarks>
        public unsafe static  ILRetArray<complex>  sqrtc (ILInArray< double > A) {
            using (ILScope.Enter(A)) {
                if (A.IsEmpty)
                    return new  ILRetArray<complex>(A.Size); 
                ILSize inDim = A.Size;
                double[] arrA = A.GetArrayForRead(); 
                complex [] retArr;
                int outLen = inDim.NumberOfElements;
                bool inplace = true;
               
                if (true) {
                    retArr = ILMemoryPool.Pool.New<complex>(outLen);
                    inplace = false; 
                }
                int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength, workerCount = 1;
                if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
                    if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
                        workItemLength = outLen / workItemCount;
                        //workItemLength = (int)((double)outLen / workItemCount * 1.05);
                    } else {
                        workItemLength = outLen / 2;
                        workItemCount = 2;
                    }
                } else {
                    workItemLength = outLen;
                    workItemCount = 1;
                }
                ILDenseStorage<complex> retStorage = new ILDenseStorage<complex>(retArr, inDim);
                
                Action<object> worker = data => {
                    Tuple<int, int, IntPtr, IntPtr, bool> range = (Tuple<int, int, IntPtr, IntPtr, bool>)data;
                   
                    complex* cp = ((complex*)range.Item4 + range.Item1);
                    int len = range.Item2; 
                    if (range.Item5) { 
                        // inplace
                        while (len > 20) {
                            cp[0] =  complex.Sqrt(cp[0]  )  /*dummy*/;
                            cp[1] =  complex.Sqrt(cp[1]  )  /*dummy*/;
                            cp[2] =  complex.Sqrt(cp[2]  )  /*dummy*/;
                            cp[3] =  complex.Sqrt(cp[3]  )  /*dummy*/;
                            cp[4] =  complex.Sqrt(cp[4]  )  /*dummy*/;
                            cp[5] =  complex.Sqrt(cp[5]  )  /*dummy*/;
                            cp[6] =  complex.Sqrt(cp[6]  )  /*dummy*/;
                            cp[7] =  complex.Sqrt(cp[7]  )  /*dummy*/;
                            cp[8] =  complex.Sqrt(cp[8]  )  /*dummy*/;
                            cp[9] =  complex.Sqrt(cp[9]  )  /*dummy*/;
                            cp[10] =  complex.Sqrt(cp[10]  )  /*dummy*/;
                            cp[11] =  complex.Sqrt(cp[11]  )  /*dummy*/;
                            cp[12] =  complex.Sqrt(cp[12]  )  /*dummy*/;
                            cp[13] =  complex.Sqrt(cp[13]  )  /*dummy*/;
                            cp[14] =  complex.Sqrt(cp[14]  )  /*dummy*/;
                            cp[15] =  complex.Sqrt(cp[15]  )  /*dummy*/;
                            cp[16] =  complex.Sqrt(cp[16]  )  /*dummy*/;
                            cp[17] =  complex.Sqrt(cp[17]  )  /*dummy*/;
                            cp[18] =  complex.Sqrt(cp[18]  )  /*dummy*/;
                            cp[19] =  complex.Sqrt(cp[19]  )  /*dummy*/;
                            cp[20] =  complex.Sqrt(cp[20]  )  /*dummy*/;
                            cp+=21; len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  complex.Sqrt(*cp  )  /*dummy*/;
                            cp++;
                        }
                    } else {
                        double* ap = ((double*)range.Item3 + range.Item1);
                        while (len > 20) {
                            cp[0] =  complex.Sqrt(ap[0]  )  /*dummy*/;
                            cp[1] =  complex.Sqrt(ap[1]  )  /*dummy*/;
                            cp[2] =  complex.Sqrt(ap[2]  )  /*dummy*/;
                            cp[3] =  complex.Sqrt(ap[3]  )  /*dummy*/;
                            cp[4] =  complex.Sqrt(ap[4]  )  /*dummy*/;
                            cp[5] =  complex.Sqrt(ap[5]  )  /*dummy*/;
                            cp[6] =  complex.Sqrt(ap[6]  )  /*dummy*/;
                            cp[7] =  complex.Sqrt(ap[7]  )  /*dummy*/;
                            cp[8] =  complex.Sqrt(ap[8]  )  /*dummy*/;
                            cp[9] =  complex.Sqrt(ap[9]  )  /*dummy*/;
                            cp[10] =  complex.Sqrt(ap[10]  )  /*dummy*/;
                            cp[11] =  complex.Sqrt(ap[11]  )  /*dummy*/;
                            cp[12] =  complex.Sqrt(ap[12]  )  /*dummy*/;
                            cp[13] =  complex.Sqrt(ap[13]  )  /*dummy*/;
                            cp[14] =  complex.Sqrt(ap[14]  )  /*dummy*/;
                            cp[15] =  complex.Sqrt(ap[15]  )  /*dummy*/;
                            cp[16] =  complex.Sqrt(ap[16]  )  /*dummy*/;
                            cp[17] =  complex.Sqrt(ap[17]  )  /*dummy*/;
                            cp[18] =  complex.Sqrt(ap[18]  )  /*dummy*/;
                            cp[19] =  complex.Sqrt(ap[19]  )  /*dummy*/;
                            cp[20] =  complex.Sqrt(ap[20]  )  /*dummy*/;
                            ap += 21;
                            cp += 21;
                            len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  complex.Sqrt(*ap  )  /*dummy*/;
                            ap++;
                            cp++;
                        }
                    }
                    System.Threading.Interlocked.Decrement(ref workerCount);
                };

                fixed ( double* arrAP = arrA)
                fixed ( complex* retArrP = retArr) {
                    for (; i < workItemCount - 1; i++) {
                        Tuple<int, int, IntPtr, IntPtr, bool> range
                            = new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace);
                        System.Threading.Interlocked.Increment(ref workerCount);
                        ILThreadPool.QueueUserWorkItem(i,worker, range);
                    }
                    // the last (or may the only) chunk is done right here
                    worker(new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace));

                    ILThreadPool.Wait4Workers(ref workerCount);
                }
                return new  ILRetArray<complex>(retStorage);
            }
        }
        /// <summary>Square root of array elements</summary>
        /// <param name="A">Input array</param>
        /// <returns>Square root of elements of A </returns>
        /// <remarks><para>If the input array is empty, an empty array will be returned.</para>
        /// <para>The array returned will be a dense array.</para></remarks>
        public unsafe static  ILRetArray<complex>  sqrt (ILInArray< complex > A) {
            using (ILScope.Enter(A)) {
                if (A.IsEmpty)
                    return new  ILRetArray<complex>(A.Size); 
                ILSize inDim = A.Size;
                complex[] arrA = A.GetArrayForRead(); 
                complex [] retArr;
                int outLen = inDim.NumberOfElements;
                bool inplace = true;
               
                if (true) {
                    retArr = ILMemoryPool.Pool.New<complex>(outLen);
                    inplace = false; 
                }
                int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength, workerCount = 1;
                if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
                    if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
                        workItemLength = outLen / workItemCount;
                        //workItemLength = (int)((double)outLen / workItemCount * 1.05);
                    } else {
                        workItemLength = outLen / 2;
                        workItemCount = 2;
                    }
                } else {
                    workItemLength = outLen;
                    workItemCount = 1;
                }
                ILDenseStorage<complex> retStorage = new ILDenseStorage<complex>(retArr, inDim);
                
                Action<object> worker = data => {
                    Tuple<int, int, IntPtr, IntPtr, bool> range = (Tuple<int, int, IntPtr, IntPtr, bool>)data;
                   
                    complex* cp = ((complex*)range.Item4 + range.Item1);
                    int len = range.Item2; 
                    if (range.Item5) { 
                        // inplace
                        while (len > 20) {
                            cp[0] =  complex.Sqrt(cp[0]  )  /*dummy*/;
                            cp[1] =  complex.Sqrt(cp[1]  )  /*dummy*/;
                            cp[2] =  complex.Sqrt(cp[2]  )  /*dummy*/;
                            cp[3] =  complex.Sqrt(cp[3]  )  /*dummy*/;
                            cp[4] =  complex.Sqrt(cp[4]  )  /*dummy*/;
                            cp[5] =  complex.Sqrt(cp[5]  )  /*dummy*/;
                            cp[6] =  complex.Sqrt(cp[6]  )  /*dummy*/;
                            cp[7] =  complex.Sqrt(cp[7]  )  /*dummy*/;
                            cp[8] =  complex.Sqrt(cp[8]  )  /*dummy*/;
                            cp[9] =  complex.Sqrt(cp[9]  )  /*dummy*/;
                            cp[10] =  complex.Sqrt(cp[10]  )  /*dummy*/;
                            cp[11] =  complex.Sqrt(cp[11]  )  /*dummy*/;
                            cp[12] =  complex.Sqrt(cp[12]  )  /*dummy*/;
                            cp[13] =  complex.Sqrt(cp[13]  )  /*dummy*/;
                            cp[14] =  complex.Sqrt(cp[14]  )  /*dummy*/;
                            cp[15] =  complex.Sqrt(cp[15]  )  /*dummy*/;
                            cp[16] =  complex.Sqrt(cp[16]  )  /*dummy*/;
                            cp[17] =  complex.Sqrt(cp[17]  )  /*dummy*/;
                            cp[18] =  complex.Sqrt(cp[18]  )  /*dummy*/;
                            cp[19] =  complex.Sqrt(cp[19]  )  /*dummy*/;
                            cp[20] =  complex.Sqrt(cp[20]  )  /*dummy*/;
                            cp+=21; len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  complex.Sqrt(*cp  )  /*dummy*/;
                            cp++;
                        }
                    } else {
                        complex* ap = ((complex*)range.Item3 + range.Item1);
                        while (len > 20) {
                            cp[0] =  complex.Sqrt(ap[0]  )  /*dummy*/;
                            cp[1] =  complex.Sqrt(ap[1]  )  /*dummy*/;
                            cp[2] =  complex.Sqrt(ap[2]  )  /*dummy*/;
                            cp[3] =  complex.Sqrt(ap[3]  )  /*dummy*/;
                            cp[4] =  complex.Sqrt(ap[4]  )  /*dummy*/;
                            cp[5] =  complex.Sqrt(ap[5]  )  /*dummy*/;
                            cp[6] =  complex.Sqrt(ap[6]  )  /*dummy*/;
                            cp[7] =  complex.Sqrt(ap[7]  )  /*dummy*/;
                            cp[8] =  complex.Sqrt(ap[8]  )  /*dummy*/;
                            cp[9] =  complex.Sqrt(ap[9]  )  /*dummy*/;
                            cp[10] =  complex.Sqrt(ap[10]  )  /*dummy*/;
                            cp[11] =  complex.Sqrt(ap[11]  )  /*dummy*/;
                            cp[12] =  complex.Sqrt(ap[12]  )  /*dummy*/;
                            cp[13] =  complex.Sqrt(ap[13]  )  /*dummy*/;
                            cp[14] =  complex.Sqrt(ap[14]  )  /*dummy*/;
                            cp[15] =  complex.Sqrt(ap[15]  )  /*dummy*/;
                            cp[16] =  complex.Sqrt(ap[16]  )  /*dummy*/;
                            cp[17] =  complex.Sqrt(ap[17]  )  /*dummy*/;
                            cp[18] =  complex.Sqrt(ap[18]  )  /*dummy*/;
                            cp[19] =  complex.Sqrt(ap[19]  )  /*dummy*/;
                            cp[20] =  complex.Sqrt(ap[20]  )  /*dummy*/;
                            ap += 21;
                            cp += 21;
                            len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  complex.Sqrt(*ap  )  /*dummy*/;
                            ap++;
                            cp++;
                        }
                    }
                    System.Threading.Interlocked.Decrement(ref workerCount);
                };

                fixed ( complex* arrAP = arrA)
                fixed ( complex* retArrP = retArr) {
                    for (; i < workItemCount - 1; i++) {
                        Tuple<int, int, IntPtr, IntPtr, bool> range
                            = new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace);
                        System.Threading.Interlocked.Increment(ref workerCount);
                        ILThreadPool.QueueUserWorkItem(i,worker, range);
                    }
                    // the last (or may the only) chunk is done right here
                    worker(new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace));

                    ILThreadPool.Wait4Workers(ref workerCount);
                }
                return new  ILRetArray<complex>(retStorage);
            }
        }
        /// <summary>Square root of array elements</summary>
        /// <param name="A">Input array</param>
        /// <returns>Square root of elements of A </returns>
        /// <remarks><para>If the input array is empty, an empty array will be returned.</para>
        /// <para>The array returned will be a dense array.</para></remarks>
        public unsafe static  ILRetArray<fcomplex>  sqrt (ILInArray< fcomplex > A) {
            using (ILScope.Enter(A)) {
                if (A.IsEmpty)
                    return new  ILRetArray<fcomplex>(A.Size); 
                ILSize inDim = A.Size;
                fcomplex[] arrA = A.GetArrayForRead(); 
                fcomplex [] retArr;
                int outLen = inDim.NumberOfElements;
                bool inplace = true;
               
                if (true) {
                    retArr = ILMemoryPool.Pool.New<fcomplex>(outLen);
                    inplace = false; 
                }
                int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength, workerCount = 1;
                if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
                    if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
                        workItemLength = outLen / workItemCount;
                        //workItemLength = (int)((double)outLen / workItemCount * 1.05);
                    } else {
                        workItemLength = outLen / 2;
                        workItemCount = 2;
                    }
                } else {
                    workItemLength = outLen;
                    workItemCount = 1;
                }
                ILDenseStorage<fcomplex> retStorage = new ILDenseStorage<fcomplex>(retArr, inDim);
                
                Action<object> worker = data => {
                    Tuple<int, int, IntPtr, IntPtr, bool> range = (Tuple<int, int, IntPtr, IntPtr, bool>)data;
                   
                    fcomplex* cp = ((fcomplex*)range.Item4 + range.Item1);
                    int len = range.Item2; 
                    if (range.Item5) { 
                        // inplace
                        while (len > 20) {
                            cp[0] =  fcomplex.Sqrt(cp[0]  )  /*dummy*/;
                            cp[1] =  fcomplex.Sqrt(cp[1]  )  /*dummy*/;
                            cp[2] =  fcomplex.Sqrt(cp[2]  )  /*dummy*/;
                            cp[3] =  fcomplex.Sqrt(cp[3]  )  /*dummy*/;
                            cp[4] =  fcomplex.Sqrt(cp[4]  )  /*dummy*/;
                            cp[5] =  fcomplex.Sqrt(cp[5]  )  /*dummy*/;
                            cp[6] =  fcomplex.Sqrt(cp[6]  )  /*dummy*/;
                            cp[7] =  fcomplex.Sqrt(cp[7]  )  /*dummy*/;
                            cp[8] =  fcomplex.Sqrt(cp[8]  )  /*dummy*/;
                            cp[9] =  fcomplex.Sqrt(cp[9]  )  /*dummy*/;
                            cp[10] =  fcomplex.Sqrt(cp[10]  )  /*dummy*/;
                            cp[11] =  fcomplex.Sqrt(cp[11]  )  /*dummy*/;
                            cp[12] =  fcomplex.Sqrt(cp[12]  )  /*dummy*/;
                            cp[13] =  fcomplex.Sqrt(cp[13]  )  /*dummy*/;
                            cp[14] =  fcomplex.Sqrt(cp[14]  )  /*dummy*/;
                            cp[15] =  fcomplex.Sqrt(cp[15]  )  /*dummy*/;
                            cp[16] =  fcomplex.Sqrt(cp[16]  )  /*dummy*/;
                            cp[17] =  fcomplex.Sqrt(cp[17]  )  /*dummy*/;
                            cp[18] =  fcomplex.Sqrt(cp[18]  )  /*dummy*/;
                            cp[19] =  fcomplex.Sqrt(cp[19]  )  /*dummy*/;
                            cp[20] =  fcomplex.Sqrt(cp[20]  )  /*dummy*/;
                            cp+=21; len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  fcomplex.Sqrt(*cp  )  /*dummy*/;
                            cp++;
                        }
                    } else {
                        fcomplex* ap = ((fcomplex*)range.Item3 + range.Item1);
                        while (len > 20) {
                            cp[0] =  fcomplex.Sqrt(ap[0]  )  /*dummy*/;
                            cp[1] =  fcomplex.Sqrt(ap[1]  )  /*dummy*/;
                            cp[2] =  fcomplex.Sqrt(ap[2]  )  /*dummy*/;
                            cp[3] =  fcomplex.Sqrt(ap[3]  )  /*dummy*/;
                            cp[4] =  fcomplex.Sqrt(ap[4]  )  /*dummy*/;
                            cp[5] =  fcomplex.Sqrt(ap[5]  )  /*dummy*/;
                            cp[6] =  fcomplex.Sqrt(ap[6]  )  /*dummy*/;
                            cp[7] =  fcomplex.Sqrt(ap[7]  )  /*dummy*/;
                            cp[8] =  fcomplex.Sqrt(ap[8]  )  /*dummy*/;
                            cp[9] =  fcomplex.Sqrt(ap[9]  )  /*dummy*/;
                            cp[10] =  fcomplex.Sqrt(ap[10]  )  /*dummy*/;
                            cp[11] =  fcomplex.Sqrt(ap[11]  )  /*dummy*/;
                            cp[12] =  fcomplex.Sqrt(ap[12]  )  /*dummy*/;
                            cp[13] =  fcomplex.Sqrt(ap[13]  )  /*dummy*/;
                            cp[14] =  fcomplex.Sqrt(ap[14]  )  /*dummy*/;
                            cp[15] =  fcomplex.Sqrt(ap[15]  )  /*dummy*/;
                            cp[16] =  fcomplex.Sqrt(ap[16]  )  /*dummy*/;
                            cp[17] =  fcomplex.Sqrt(ap[17]  )  /*dummy*/;
                            cp[18] =  fcomplex.Sqrt(ap[18]  )  /*dummy*/;
                            cp[19] =  fcomplex.Sqrt(ap[19]  )  /*dummy*/;
                            cp[20] =  fcomplex.Sqrt(ap[20]  )  /*dummy*/;
                            ap += 21;
                            cp += 21;
                            len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  fcomplex.Sqrt(*ap  )  /*dummy*/;
                            ap++;
                            cp++;
                        }
                    }
                    System.Threading.Interlocked.Decrement(ref workerCount);
                };

                fixed ( fcomplex* arrAP = arrA)
                fixed ( fcomplex* retArrP = retArr) {
                    for (; i < workItemCount - 1; i++) {
                        Tuple<int, int, IntPtr, IntPtr, bool> range
                            = new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace);
                        System.Threading.Interlocked.Increment(ref workerCount);
                        ILThreadPool.QueueUserWorkItem(i,worker, range);
                    }
                    // the last (or may the only) chunk is done right here
                    worker(new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace));

                    ILThreadPool.Wait4Workers(ref workerCount);
                }
                return new  ILRetArray<fcomplex>(retStorage);
            }
        }
        /// <summary>Square root of array elements - real output</summary>
        /// <param name="A">Input array - only positive values are allowed.</param>
        /// <returns>Square root of elements of A - real output</returns>
        /// <remarks><para>If the input array is empty, an empty array will be returned.</para>
        /// <para>The array returned will be a dense array.</para></remarks>
        public unsafe static  ILRetArray<float>  sqrt (ILInArray< float > A) {
            using (ILScope.Enter(A)) {
                if (A.IsEmpty)
                    return new  ILRetArray<float>(A.Size); 
                ILSize inDim = A.Size;
                float[] arrA = A.GetArrayForRead(); 
                float [] retArr;
                int outLen = inDim.NumberOfElements;
                bool inplace = true;
               
                if (true) {
                    retArr = ILMemoryPool.Pool.New<float>(outLen);
                    inplace = false; 
                }
                int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength, workerCount = 1;
                if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
                    if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
                        workItemLength = outLen / workItemCount;
                        //workItemLength = (int)((double)outLen / workItemCount * 1.05);
                    } else {
                        workItemLength = outLen / 2;
                        workItemCount = 2;
                    }
                } else {
                    workItemLength = outLen;
                    workItemCount = 1;
                }
                ILDenseStorage<float> retStorage = new ILDenseStorage<float>(retArr, inDim);
                
                Action<object> worker = data => {
                    Tuple<int, int, IntPtr, IntPtr, bool> range = (Tuple<int, int, IntPtr, IntPtr, bool>)data;
                   
                    float* cp = ((float*)range.Item4 + range.Item1);
                    int len = range.Item2; 
                    if (range.Item5) { 
                        // inplace
                        while (len > 20) {
                            cp[0] =  (float)Math.Sqrt(cp[0]  )  /*dummy*/;
                            cp[1] =  (float)Math.Sqrt(cp[1]  )  /*dummy*/;
                            cp[2] =  (float)Math.Sqrt(cp[2]  )  /*dummy*/;
                            cp[3] =  (float)Math.Sqrt(cp[3]  )  /*dummy*/;
                            cp[4] =  (float)Math.Sqrt(cp[4]  )  /*dummy*/;
                            cp[5] =  (float)Math.Sqrt(cp[5]  )  /*dummy*/;
                            cp[6] =  (float)Math.Sqrt(cp[6]  )  /*dummy*/;
                            cp[7] =  (float)Math.Sqrt(cp[7]  )  /*dummy*/;
                            cp[8] =  (float)Math.Sqrt(cp[8]  )  /*dummy*/;
                            cp[9] =  (float)Math.Sqrt(cp[9]  )  /*dummy*/;
                            cp[10] =  (float)Math.Sqrt(cp[10]  )  /*dummy*/;
                            cp[11] =  (float)Math.Sqrt(cp[11]  )  /*dummy*/;
                            cp[12] =  (float)Math.Sqrt(cp[12]  )  /*dummy*/;
                            cp[13] =  (float)Math.Sqrt(cp[13]  )  /*dummy*/;
                            cp[14] =  (float)Math.Sqrt(cp[14]  )  /*dummy*/;
                            cp[15] =  (float)Math.Sqrt(cp[15]  )  /*dummy*/;
                            cp[16] =  (float)Math.Sqrt(cp[16]  )  /*dummy*/;
                            cp[17] =  (float)Math.Sqrt(cp[17]  )  /*dummy*/;
                            cp[18] =  (float)Math.Sqrt(cp[18]  )  /*dummy*/;
                            cp[19] =  (float)Math.Sqrt(cp[19]  )  /*dummy*/;
                            cp[20] =  (float)Math.Sqrt(cp[20]  )  /*dummy*/;
                            cp+=21; len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  (float)Math.Sqrt(*cp  )  /*dummy*/;
                            cp++;
                        }
                    } else {
                        float* ap = ((float*)range.Item3 + range.Item1);
                        while (len > 20) {
                            cp[0] =  (float)Math.Sqrt(ap[0]  )  /*dummy*/;
                            cp[1] =  (float)Math.Sqrt(ap[1]  )  /*dummy*/;
                            cp[2] =  (float)Math.Sqrt(ap[2]  )  /*dummy*/;
                            cp[3] =  (float)Math.Sqrt(ap[3]  )  /*dummy*/;
                            cp[4] =  (float)Math.Sqrt(ap[4]  )  /*dummy*/;
                            cp[5] =  (float)Math.Sqrt(ap[5]  )  /*dummy*/;
                            cp[6] =  (float)Math.Sqrt(ap[6]  )  /*dummy*/;
                            cp[7] =  (float)Math.Sqrt(ap[7]  )  /*dummy*/;
                            cp[8] =  (float)Math.Sqrt(ap[8]  )  /*dummy*/;
                            cp[9] =  (float)Math.Sqrt(ap[9]  )  /*dummy*/;
                            cp[10] =  (float)Math.Sqrt(ap[10]  )  /*dummy*/;
                            cp[11] =  (float)Math.Sqrt(ap[11]  )  /*dummy*/;
                            cp[12] =  (float)Math.Sqrt(ap[12]  )  /*dummy*/;
                            cp[13] =  (float)Math.Sqrt(ap[13]  )  /*dummy*/;
                            cp[14] =  (float)Math.Sqrt(ap[14]  )  /*dummy*/;
                            cp[15] =  (float)Math.Sqrt(ap[15]  )  /*dummy*/;
                            cp[16] =  (float)Math.Sqrt(ap[16]  )  /*dummy*/;
                            cp[17] =  (float)Math.Sqrt(ap[17]  )  /*dummy*/;
                            cp[18] =  (float)Math.Sqrt(ap[18]  )  /*dummy*/;
                            cp[19] =  (float)Math.Sqrt(ap[19]  )  /*dummy*/;
                            cp[20] =  (float)Math.Sqrt(ap[20]  )  /*dummy*/;
                            ap += 21;
                            cp += 21;
                            len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  (float)Math.Sqrt(*ap  )  /*dummy*/;
                            ap++;
                            cp++;
                        }
                    }
                    System.Threading.Interlocked.Decrement(ref workerCount);
                };

                fixed ( float* arrAP = arrA)
                fixed ( float* retArrP = retArr) {
                    for (; i < workItemCount - 1; i++) {
                        Tuple<int, int, IntPtr, IntPtr, bool> range
                            = new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace);
                        System.Threading.Interlocked.Increment(ref workerCount);
                        ILThreadPool.QueueUserWorkItem(i,worker, range);
                    }
                    // the last (or may the only) chunk is done right here
                    worker(new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace));

                    ILThreadPool.Wait4Workers(ref workerCount);
                }
                return new  ILRetArray<float>(retStorage);
            }
        }
        /// <summary>Square root of array elements - real output</summary>
        /// <param name="A">Input array - only positive values are allowed.</param>
        /// <returns>Square root of elements of A - real output</returns>
        /// <remarks><para>If the input array is empty, an empty array will be returned.</para>
        /// <para>The array returned will be a dense array.</para></remarks>
        public unsafe static  ILRetArray<double>  sqrt (ILInArray< double > A) {
            using (ILScope.Enter(A)) {
                if (A.IsEmpty)
                    return new  ILRetArray<double>(A.Size); 
                ILSize inDim = A.Size;
                double[] arrA = A.GetArrayForRead(); 
                double [] retArr;
                int outLen = inDim.NumberOfElements;
                bool inplace = true;
               
                if (true) {
                    retArr = ILMemoryPool.Pool.New<double>(outLen);
                    inplace = false; 
                }
                int i = 0, workItemCount = Settings.s_maxNumberThreads, workItemLength, workerCount = 1;
                if (Settings.s_maxNumberThreads > 1 && outLen / 2 > Settings.s_minParallelElement1Count) {
                    if (outLen / workItemCount > Settings.s_minParallelElement1Count) {
                        workItemLength = outLen / workItemCount;
                        //workItemLength = (int)((double)outLen / workItemCount * 1.05);
                    } else {
                        workItemLength = outLen / 2;
                        workItemCount = 2;
                    }
                } else {
                    workItemLength = outLen;
                    workItemCount = 1;
                }
                ILDenseStorage<double> retStorage = new ILDenseStorage<double>(retArr, inDim);
                
                Action<object> worker = data => {
                    Tuple<int, int, IntPtr, IntPtr, bool> range = (Tuple<int, int, IntPtr, IntPtr, bool>)data;
                   
                    double* cp = ((double*)range.Item4 + range.Item1);
                    int len = range.Item2; 
                    if (range.Item5) { 
                        // inplace
                        while (len > 20) {
                            cp[0] =  Math.Sqrt(cp[0]  )  /*dummy*/;
                            cp[1] =  Math.Sqrt(cp[1]  )  /*dummy*/;
                            cp[2] =  Math.Sqrt(cp[2]  )  /*dummy*/;
                            cp[3] =  Math.Sqrt(cp[3]  )  /*dummy*/;
                            cp[4] =  Math.Sqrt(cp[4]  )  /*dummy*/;
                            cp[5] =  Math.Sqrt(cp[5]  )  /*dummy*/;
                            cp[6] =  Math.Sqrt(cp[6]  )  /*dummy*/;
                            cp[7] =  Math.Sqrt(cp[7]  )  /*dummy*/;
                            cp[8] =  Math.Sqrt(cp[8]  )  /*dummy*/;
                            cp[9] =  Math.Sqrt(cp[9]  )  /*dummy*/;
                            cp[10] =  Math.Sqrt(cp[10]  )  /*dummy*/;
                            cp[11] =  Math.Sqrt(cp[11]  )  /*dummy*/;
                            cp[12] =  Math.Sqrt(cp[12]  )  /*dummy*/;
                            cp[13] =  Math.Sqrt(cp[13]  )  /*dummy*/;
                            cp[14] =  Math.Sqrt(cp[14]  )  /*dummy*/;
                            cp[15] =  Math.Sqrt(cp[15]  )  /*dummy*/;
                            cp[16] =  Math.Sqrt(cp[16]  )  /*dummy*/;
                            cp[17] =  Math.Sqrt(cp[17]  )  /*dummy*/;
                            cp[18] =  Math.Sqrt(cp[18]  )  /*dummy*/;
                            cp[19] =  Math.Sqrt(cp[19]  )  /*dummy*/;
                            cp[20] =  Math.Sqrt(cp[20]  )  /*dummy*/;
                            cp+=21; len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  Math.Sqrt(*cp  )  /*dummy*/;
                            cp++;
                        }
                    } else {
                        double* ap = ((double*)range.Item3 + range.Item1);
                        while (len > 20) {
                            cp[0] =  Math.Sqrt(ap[0]  )  /*dummy*/;
                            cp[1] =  Math.Sqrt(ap[1]  )  /*dummy*/;
                            cp[2] =  Math.Sqrt(ap[2]  )  /*dummy*/;
                            cp[3] =  Math.Sqrt(ap[3]  )  /*dummy*/;
                            cp[4] =  Math.Sqrt(ap[4]  )  /*dummy*/;
                            cp[5] =  Math.Sqrt(ap[5]  )  /*dummy*/;
                            cp[6] =  Math.Sqrt(ap[6]  )  /*dummy*/;
                            cp[7] =  Math.Sqrt(ap[7]  )  /*dummy*/;
                            cp[8] =  Math.Sqrt(ap[8]  )  /*dummy*/;
                            cp[9] =  Math.Sqrt(ap[9]  )  /*dummy*/;
                            cp[10] =  Math.Sqrt(ap[10]  )  /*dummy*/;
                            cp[11] =  Math.Sqrt(ap[11]  )  /*dummy*/;
                            cp[12] =  Math.Sqrt(ap[12]  )  /*dummy*/;
                            cp[13] =  Math.Sqrt(ap[13]  )  /*dummy*/;
                            cp[14] =  Math.Sqrt(ap[14]  )  /*dummy*/;
                            cp[15] =  Math.Sqrt(ap[15]  )  /*dummy*/;
                            cp[16] =  Math.Sqrt(ap[16]  )  /*dummy*/;
                            cp[17] =  Math.Sqrt(ap[17]  )  /*dummy*/;
                            cp[18] =  Math.Sqrt(ap[18]  )  /*dummy*/;
                            cp[19] =  Math.Sqrt(ap[19]  )  /*dummy*/;
                            cp[20] =  Math.Sqrt(ap[20]  )  /*dummy*/;
                            ap += 21;
                            cp += 21;
                            len -= 21; 
                        }
                        while (len-- > 0) {
                            *cp =  Math.Sqrt(*ap  )  /*dummy*/;
                            ap++;
                            cp++;
                        }
                    }
                    System.Threading.Interlocked.Decrement(ref workerCount);
                };

                fixed ( double* arrAP = arrA)
                fixed ( double* retArrP = retArr) {
                    for (; i < workItemCount - 1; i++) {
                        Tuple<int, int, IntPtr, IntPtr, bool> range
                            = new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace);
                        System.Threading.Interlocked.Increment(ref workerCount);
                        ILThreadPool.QueueUserWorkItem(i,worker, range);
                    }
                    // the last (or may the only) chunk is done right here
                    worker(new Tuple<int, int, IntPtr, IntPtr, bool>
                                (i * workItemLength, outLen - i * workItemLength, (IntPtr)arrAP, (IntPtr)retArrP, inplace));

                    ILThreadPool.Wait4Workers(ref workerCount);
                }
                return new  ILRetArray<double>(retStorage);
            }
        }

#endregion HYCALPER AUTO GENERATED CODE

    }
}