private void RatesButton_Click(object sender, EventArgs e)
{
if (openFileDialog2.ShowDialog() == DialogResult.OK)
{
int numFiles = openFileDialog2.FileNames.Length;
for (int i = 0; i < numFiles; i++)
{
//string msg = "";
//string latex = "";
textBox1.Text += "Loading file " + openFileDialog2.FileNames[i] + "\r\n";
JOS.Multiplet EmiMult = new JOS.Multiplet();
EmiMult.GetOmegas(exp);
EmiMult.LoadEmiFromFile(openFileDialog2.FileNames[i]);
EmiMult.CalculateRates();
EmiMult.DumpEmidata(out string msg);
textBox1.Text += msg;
EmiMult.ReportRates(out msg, out string latex);
textBox1.Text += msg + "\r\n";
textBox2.Text += latex + "\r\n";
Console.WriteLine(msg);
}
string uniq = DateTime.Now.ToString("yyyyMMddHHmmss");
System.IO.File.WriteAllText("emi" + uniq + exp.absofilename + ".log", textBox1.Text);
}
}
static void Main(string[] args)
{
Console.WriteLine(Globals.qe.ToString("g4", System.Globalization.CultureInfo.InvariantCulture));
JOS.Solver testowy = new JOS.Solver();
testowy.dumpVersion();
Console.WriteLine(testowy.f(0.1, 0.1, 0.1, 450e-7, 1.85, 9, 1e-20, 1e-20, 1e-20).ToString("g4", System.Globalization.CultureInfo.InvariantCulture));
JOS.Multiplet exp = new JOS.Multiplet();
//
exp.o2 = 1.0e-20;
exp.o4 = 1.1e-20;
exp.o6 = 1.2e-20;
exp.LoadFromFile(@"r:\testszarp.txt");
//Console.WriteLine( testowy.chi2(exp));
//MathNet.Numerics.LinearAlgebra.Matrix<double> Hess ;
//MathNet.Numerics.LinearAlgebra.Matrix<double> Grad;
//testowy.CalculateHessian(exp, out Hess,out Grad);
//Console.Write(Hess.ToString());
string msg, lat;
msg = "";
testowy.FitLM(exp, out msg, out lat);
//Console.WriteLine(msg);
exp.CalculateRates();
exp.ReportRates(out msg);
//Console.WriteLine(exp.o2);
Console.WriteLine("Report");
Console.WriteLine(msg);
}
public void Fit_Works()
{
string s1, s2;
JOS.Solver tests = new JOS.Solver();
JOS.Multiplet testm = new JOS.Multiplet();
testm.u2 = new System.Collections.Generic.List <double>();
testm.u4 = new System.Collections.Generic.List <double>();
testm.u6 = new System.Collections.Generic.List <double>();
testm.fexp = new System.Collections.Generic.List <double>();
testm.lambda0 = new System.Collections.Generic.List <double>();
testm.n = 1.95;
testm.TwoJPlusOne = 13;
testm.u2.Add(0.001); testm.u2.Add(0.002); testm.u2.Add(0.004); testm.u2.Add(0.301);
testm.u4.Add(0.001); testm.u4.Add(0.343); testm.u4.Add(0.054); testm.u4.Add(0.101);
testm.u6.Add(0.001); testm.u6.Add(0.212); testm.u6.Add(0.404); testm.u6.Add(0.501);
testm.fexp.Add(1e-6); testm.fexp.Add(3.56e-6); testm.fexp.Add(2.34e-6); testm.fexp.Add(0.34e-6);
testm.lambda0.Add(345e-7); testm.lambda0.Add(645e-7); testm.lambda0.Add(349e-7); testm.lambda0.Add(475e-7);
tests.FitLM(testm, out s1, out s2);
double testchi2 = tests.Chi2(testm);
Assert.IsTrue((testchi2 < 1e-7), "Słaba zbieżność");
}
public Form1()
{
InitializeComponent();
exp = new JOS.Multiplet();
solver = new JOS.Solver();
}
