WCF Link - Fitting

// <copyright file="WCFFitting.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.Collections.Generic;
    using System.Linq;
    using RelaxIS_SDK.Common;
    using RelaxIS_SDK.RelaxISWCFLink;

    /// <summary>
    /// Implements examples of using the WCF fit functions.
    /// </summary>
    public static class WCFFitting
    {
        /// <summary>
        /// An example of performing a single spectrum fit via the WCF Link.
        /// </summary>
        public static void PerformSingleFit()
        {
            // Get the interface object
            using (var wcf = RelaxISWCFLinkProxy.GetDefaultProxy())
            {
                // Create a spectrum that contains the data
                var spectrum = new ImpedanceSpectrum(CommonFunctions.GetExampleData());

                // Define fit model
                spectrum.Analysis.Model = "R-(R)(C)";

                // Get parameters for this model.
                var parameters = wcf.GetDefaultParameters(spectrum.Analysis.Model);

                // Set initial parameter values
                parameters[0].Value = 40;
                parameters[1].Value = 70;
                parameters[2].Value = 2e-6;

                spectrum.Analysis.Parameters.AddRange(parameters);

                // Empty list to signal that all parameters are active. This overrides the fixed indizes in the parameter list.
                // If null is used as the argument in the GetFit function, the spectrum parameter's fixed property is checked instead.
                var fixedIndices = new List<int>();

                // Check the Internal Name of the available transfer functions and weight modes in the RelaxIS Setup -> Plugins list.
                var transferFunction = "Impedance";
                var weightMode = "Proportional Weighting";

                var fitResult = wcf.GetFit(spectrum, transferFunction, weightMode, fixedIndices);
                OutputFitResult(wcf, fitResult, spectrum, parameters, transferFunction);
            }
        }

        /// <summary>
        ///  An example of performing an Auto Fit using the WCF interface.
        /// </summary>
        public static void PerformAutoFit()
        {
            // Get the interface object
            using (var wcf = RelaxISWCFLinkProxy.GetDefaultProxy())
            {
                // Create a spectrum that contains the data
                var spectrum = new ImpedanceSpectrum(CommonFunctions.GetExampleData());

                // Define fit model
                spectrum.Analysis.Model = "R-(R)(C)";

                // Get parameters for this model.
                var parameters = wcf.GetDefaultParameters(spectrum.Analysis.Model);

                // Set search ranges. Leave default to allow RelaxIS to find optimized limits automatically.
                parameters[0].Limits.SetLimits(1, 1000);

                spectrum.Analysis.Parameters.AddRange(parameters);

                // Check the Internal Name of the available transfer functions and weight modes in the RelaxIS Setup -> Plugins list.
                var transferFunction = "Impedance";
                var weightMode = "Proportional Weighting";

                var fitResult = wcf.GetAutoFit(spectrum, transferFunction, weightMode, 20, 50);
                OutputFitResult(wcf, fitResult, spectrum, parameters, transferFunction);
            }
        }

        private static void OutputFitResult(IRelaxISWCFLink wcf, RelaxISFitResult fitResult, ImpedanceSpectrum spectrum, IList<Fitparameter> parameters, string transferFunction)
        {
            if (fitResult.Status == RelaxISFitResult.FitStatus.FitOK)
            {
                // Output the SSR
                Console.WriteLine(string.Format("SSR: {0}", fitResult.ESS));

                // Output the parameter values
                for (var i = 0; i < fitResult.ParameterValues.Length; i++)
                {
                    Console.WriteLine(string.Format("{0} => {1}", parameters[i].Name, fitResult.ParameterValues[i]));
                }

                // Get plot data for the input frequencies and the result parameters
                var freqs = spectrum.Data.Select(d => d.Impedance.Frequency).ToArray();
                var plotData = wcf.GetPlotDataFixedFreqs(spectrum.Analysis.Model, fitResult.ParameterValues, transferFunction, freqs);
                Console.WriteLine("Plot Data:");
                for (int i = 0; i < freqs.Length; i++)
                {
                    Console.WriteLine(string.Format("{0}\t{1}\t{2}", plotData[0][i], plotData[1][i], plotData[2][i]));
                }
            }
            else
            {
                Console.WriteLine(string.Format("Fit not successful: {0}", fitResult.Status));
            }
        }
    }
}