using System;
using System.Collections.Generic;
using System.Runtime.InteropServices;
using System.Text;
using System.ComponentModel;
using System.Threading;

namespace C06PKF
{
    class Program
    {
        [StructLayout(LayoutKind.Sequential)]/* Force sequential mapping */
        public struct complex
        {
            public double r;
            public double i;
        };

        [DllImport("FLDLL214Z_nag.dll")]
        public static extern void C06PKF(ref int job, double[,] x, double[,] y,
            ref int n, complex[] work, ref int ifail);

        static void Main(string[] args)
        {
            HiPerfTimer pt = new HiPerfTimer();
            int ifail;
            int n = 5000000;
            double[,] x = new double[n,2];
            double[,] y = new double[n,2];
            complex[] work = new complex[4 * n];
         
            for (int j=0; j<n; j++) 
            {
                  x[j,0]=(double)1.0+Math.Cos(j); 
                  x[j,1]=(double)-0.5+Math.Sin(j);
            }
            

            for (int i = 0; i < n; i++)
            {
                  y[i,0] = (double)0.5 + Math.Sin(i); 
                  y[i,1] = (double)-0.25+Math.Cos(i);
            }
            
            pt.Start();

            ifail = 1;
            int job = 2;
           
            C06PKF(ref job, x, y, ref n, work, ref ifail);
            
            pt.Stop();
            Console.WriteLine("Duration: {0} sec\n", pt.Duration); 
        }
    }

    internal class HiPerfTimer
    {
        [DllImport("Kernel32.dll")]
        private static extern bool QueryPerformanceCounter(
            out long lpPerformanceCount);

        [DllImport("Kernel32.dll")]
        private static extern bool QueryPerformanceFrequency(
            out long lpFrequency);

        private long startTime, stopTime;
        private long freq;

        // Constructor

        public HiPerfTimer()
        {
            startTime = 0;
            stopTime = 0;

            if (QueryPerformanceFrequency(out freq) == false)
            {
                // high-performance counter not supported

                throw new Win32Exception();
            }
        }


        // Start the timer

        public void Start()
        {

            Thread.Sleep(0);

            QueryPerformanceCounter(out startTime);
        }


        // Stop the timer

        public void Stop()
        {
            QueryPerformanceCounter(out stopTime);
        }


        // Returns the duration of the timer (in seconds)

        public double Duration
        {
            get
            {
                return (double)(stopTime - startTime) / (double)freq;
            }
        }
    }
}