using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using CategoryTheory;
using BaseTypes;
using DiagramUI;
using DiagramUI.Interfaces;
using DiagramUI.Helpers;
using DataPerformer.Interfaces;
namespace DataPerformer.Helpers
{
/// <summary>
/// Transformer of differential equatuins
/// </summary>
public class StateTransformer : AbstractDoubleTransformer, IMeasure, ITimeMeasureProvider
{
#region Fields
IDifferentialEquationProcessor processor;
ITimeMeasureProvider provider;
double time;
double step;
#endregion
#region Ctor
protected StateTransformer(IObjectCollection collection, IDifferentialEquationProcessor processor,
ITimeMeasureProvider provider) : base(collection)
{
this.processor = processor;
this.provider = provider;
}
public StateTransformer(IObjectCollection collection,
IDifferentialEquationProcessor processor)
: this(collection, processor, null)
{
provider = this;
processor.TimeProvider = provider;
}
public StateTransformer(IObjectCollection collection)
: this(collection, new InternalRungeProcessor())
{
}
#endregion
#region IMeasure Members
Func<object> IMeasure.Parameter
{
get { return GetTime; }
}
string IMeasure.Name
{
get { return "Time"; }
}
object IMeasure.Type
{
get { return a; }
}
#endregion
#region ITimeMeasureProvider Members
IMeasure ITimeMeasureProvider.TimeMeasure
{
get { return this; }
}
double ITimeMeasureProvider.Time
{
get
{
return time;
}
set
{
time = value;
}
}
double ITimeMeasureProvider.Step
{
get
{
return step;
}
set
{
step = value;
}
}
#endregion
#region Members
protected override void Prepare()
{
if (processor == null)
{
processor = new InternalRungeProcessor();
}
base.Prepare();
processor.Set(collection);
}
/// <summary>
/// Calculation
/// </summary>
/// <param name="input">Input</param>
/// <param name="output">Output</param>
public override void Calculate(object[] input, object[] output)
{
ITimeMeasureProvider old = processor.TimeProvider;
try
{
using (new TimeProviderBackup(collection, provider, DifferentialEquationProcessor.Processor))
{
using (new ComponentCollectionBackup(collection))
{
processor.TimeProvider = provider;
// Input
double[] inp = input[0] as double[];
// Sets state vector
collection.SetStateVector(inp);
double start = provider.Time;
double step = provider.Step;
runtime.StartAll(start);
runtime.UpdateAll();
// Solution of differential equations
processor.Step(start, start + step);
provider.Time = start + step;
collection.GetStateVector(outbuffer);
// Sets final vector
output[0] = outbuffer;
}
}
}
catch (Exception ex)
{
ex.Log();
}
processor.TimeProvider = old;
}
#endregion
object GetTime()
{
return time;
}
#region Internal Runge Processor
class InternalRungeProcessor : RungeProcessor
{
#region Fields
#endregion
internal InternalRungeProcessor()
{
}
}
#endregion
}
}
Submit an article or tip