WCF Link - Analysis Functions

// <copyright file="WCFFurtherAnalysis.cs" company="rhd instruments GmbH and Co. KG">
// Copyright (c) rhd instruments GmbH and Co. KG. All rights reserved.
// Licensed under the MIT No Attribution (MIT-0) license. See section 'License' in the 'SDK Examples / Tutorials' topic for full license information.
// </copyright>

namespace RelaxIS_SDK_Examples.WCFInterface
{
    using System;
    using System.Linq;
    using RelaxIS_SDK.Common;
    using RelaxIS_SDK.RelaxISWCFLink;
    using RelaxIS_SDK.RelaxISWCFLink.DRT;
    using RelaxIS_SDK.RelaxISWCFLink.KramersKronig;

    /// <summary>
    /// Implements examples of further analysis functions such as DRT and Kramers-Kronig.
    /// </summary>
    public static class WCFFurtherAnalysis
    {
        /// <summary>
        /// Implements an example of performing a DRT analysis via the WCF Link.
        /// </summary>
        public static void PerformDRT()
        {
            // Get the interface object
            using (var wcf = RelaxISWCFLinkProxy.GetDefaultProxy())
            {
                // Create a spectrum that contains the data
                var spectrum = new ImpedanceSpectrum(CommonFunctions.GetExampleData());
                var settings = new DRTSettings()
                {
                    AutomaticNumberOfPoints = true,
                    DataType = DataType.Complex,
                    DerivativeType = DerivativeType.FirstOrder,
                    DiscretizationFunction = DiscretizationFunction.Gaussian,
                    DoInterpolate = false,
                    InterpolationCount = 0,
                    FitInductivity = true,
                    FrequencyExtrapolationFactor = 10,
                    Lambda = 0.1,
                    NumberOfPoints = 100,
                    PolarizationRemoval = PolarizationRemoval.DontRemove,
                    ShapeFactor = 0.5,
                };

                var drt = wcf.GetDRT(spectrum, settings);
                var reproduction = wcf.GetDRTReproduction(spectrum.Data.Select(d => d.Impedance.Frequency).ToList(), drt.Taus, drt.Gammas, drt.Resistance, drt.Inductivity);

                Console.WriteLine("DRT Result:");
                Console.WriteLine("Tau\tGamma");
                for (var i = 0; i < drt.Gammas.Count; i++)
                {
                    Console.WriteLine(string.Format("{0}\t{1}", drt.Taus[i], drt.Gammas[i]));
                }

                Console.WriteLine();
                Console.WriteLine("Reproduction:");
                Console.WriteLine("Frequency\tZ'\tZ''");
                for (var i = 0; i < reproduction[0].Length; i++)
                {
                    Console.WriteLine(string.Format("{0}\t{1}\t{2}", reproduction[0][i], reproduction[1][i], reproduction[2][i]));
                }
            }
        }

        /// <summary>
        /// Implements an example of performing a KK transformation via the WCF link.
        /// </summary>
        public static void KramersKronig()
        {
            // Get the interface object
            using (var wcf = RelaxISWCFLinkProxy.GetDefaultProxy())
            {
                // Create a spectrum that contains the data
                var spectrum = new ImpedanceSpectrum(CommonFunctions.GetExampleData());
                var settings = new KramersKronigSettings()
                {
                    Algorithm = KKAlgorithmType.Complex,
                    FitCapacitance = true,
                    FitInductivity = true,
                    NoTimeConstants = spectrum.Data.Count / 2,
                    UseAdmittance = false,
                };

                var kk = wcf.GetKramersKronigTransformation(spectrum, settings);
                Console.WriteLine("Kramers Kronig Result:");
                Console.WriteLine("Frequency\tZ'\tZ''");
                for (var i = 0; i < kk.ResultData.Count; i++)
                {
                    var p = kk.ResultData[i].Impedance;
                    Console.WriteLine(string.Format("{0}\t{1}\t{2}", p.Frequency, p.Real, p.Imaginary));
                }
            }
        }
    }
}