private void StartField() { chart1.Invoke(new MethodInvoker(delegate() { chart1.Series[0].Points.Clear(); })); chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.Clear(); })); // chart1.Series[0].Points.Clear(); chart2.Series[0].Points.Clear(); CloseMe = false; DateTime T1; int Ns = Convert.ToInt32(tbNs.Text); int Nt = Convert.ToInt32(tbNt.Text); SU3 Field = new SU3(Ns,Nt, Convert.ToDouble(tbBeta.Text),Convert.ToDouble(tbflux.Text), checkBox1.Checked); for (int i = 0; i < 100000; i++) { T1 = DateTime.Now; //for (int or = 0; or < 5; or++) Field.Sweep(false); //SU(2) ONLY!!! this.Text = Field.Sweep(false).ToString(); chart1.Invoke(new MethodInvoker(delegate() { chart1.Series[0].Points.AddY(Field.MeasureS()); })); chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.AddY(Field.MeasurePL()); })); lblTime1.Invoke(new MethodInvoker(delegate() { lblTime1.Text = Math.Round((DateTime.Now - T1).TotalMilliseconds).ToString() ; })); if (CloseMe) break; // if (cbDraw.Checked) DrawPicture(Field); Application.DoEvents(); } CloseMe = false; }
private void GetSample() { double Beta = Convert.ToDouble(tbBeta.Text); int Ns = Convert.ToInt32(tbNs.Text); int Nt = Convert.ToInt32(tbNt.Text); int Tsteps = Convert.ToInt32(tbTsteps.Text); int Msteps = Convert.ToInt32(tbMsteps.Text); // double[] Sample = new double[N]; SU3 Field = new SU3(Ns, Nt, Beta, Convert.ToDouble(tbflux.Text),checkBox1.Checked); double sum = 0, sum2 = 0, t = 0; int n = 0; int Mark = Convert.ToInt32(tbMark.Text); InvokeMe(() => lblProcess.Text = "Thermalizing"); for (int i = 0; i < Tsteps; i++) { Field.Sweep(false); lblProcess.Invoke(new MethodInvoker(delegate() { lblProcess.Text = "Thermalizing " + i.ToString(); })); } while (true) { Field.Sweep(false); Field.Sweep(false); //Field.Sweep(false); Field.Sweep(false); t = Field.MeasureS(); // Sample[n] = t; n++; sum += t; sum2 += t * t; InvokeMe(() => lblProcess.Text = "Measuring " + n.ToString()); // if (n % Mark == 0) { InvokeMe(() => chart2.Series[0].Points.AddY(sum / n)); InvokeMe(() => lblValue.Text = (sum / n).ToString()); } } sum = sum / n; //finding error /*double sigma2=0, Sn=0, Ssample=0; for (int i=0;i<N;i++) { Sn = 1 / (N - 1.0) * (N*sum - Sample[i]); Ssample += Sn; sigma2 += (Sn - sum) * (Sn - sum); chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.AddY(Sample[i]); })); } Ssample = Ssample / N; sigma2 = Math.Sqrt((N-1.0)/N*sigma2); Ssample = sum - (N - 1f) * (Ssample - sum); */ // int k = 10; // int d = N / k; sum2 = 0; // for (int i = 0; i < d; i++) // { // sum = 0; // for (int j = 0; j < k; j++) sum += Sample[i * k + j]; // sum /= k; // chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[1].Points.AddY(sum); })); // sum2 = (sum2 * i + sum) / ((i + 1)); // chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.AddY(sum2); })); // } // lblSampleValue.Invoke(new MethodInvoker(delegate() { lblSampleValue.Text = Ssample.ToString(); })); // lblSampleError.Invoke(new MethodInvoker(delegate() { lblSampleError.Text = sigma2.ToString(); })); //average and dispersion sum2 = Math.Sqrt(sum2 / n - sum * sum); InvokeMe(() => lblValue.Text = sum.ToString()); InvokeMe(() => lblError.Text = sum2.ToString()); }
private void button10_Click(object sender, EventArgs e) { double beta = 10; SU3 Field = new SU3(int.Parse(tbNs.Text), int.Parse(tbNt.Text), beta, 0, true);//4х2 field beta=6 flux=0 for (int i = 0; i < 10; i++) Field.Sweep(false); Dirac Dt = new Dirac(Field); Complex res = Complex.Zero; for (int n = 0; n < Dirac.N; n++) for (byte alpha = 0; alpha < 4; alpha++) for (byte a = 0; a < 3; a++) { DateTime t = DateTime.Now; Dirac D = new Dirac(Field); D.phi[n, alpha, a] = Complex.One; //double philen = Dirac.V1hermV2(D.phi, D.phi).Module(); D.Prepare(); for (int i = 0; i < 5; i++) { D.Iter(); //Vector phinew = D.MulD(D.x); //double r = 0;// D.V1hermV2(phinew, phinew).Module(); //D.V1hermV2(D.ri, D.ri).Module(); //Text = "i = " + i.ToString() + " philen=" + philen.ToString() + " =? " + r.ToString() + " " + (DateTime.Now - t).ToString(); //t = DateTime.Now; } res += D.x[n, alpha, a]; Text = n.ToString() + " " + (res / new Complex((n + 1) * 12, 0)).ToString(); Application.DoEvents(); if (CloseMe) break; } Text = (res / new Complex(Dirac.N, 0)).ToString(); //for (int exp = 0; exp < 1; exp++) //{ // double beta = 6;// 1 + 10 * exp; // SU3 Field = new SU3(4, 2, beta, 0, true);//4х2 field beta=6 flux=0 // for (int i = 0; i < 10; i++) Field.Sweep(false); // DateTime t = DateTime.Now; // int K = 100; // Vector[] chi = new Vector[K]; // Complex tr = new Complex(); // for (int k = 0; k < K; k++) // { // Dirac D = new Dirac(Field); // chi[k] = new Vector(Dirac.N, rand.Next()); // chi[k].FillGaussRandom(); // D.phi = chi[k]; // //double philen = Dirac.V1hermV2(D.phi, D.phi).Module(); // D.Prepare(); // for (int i = 0; i < 4; i++) // { // D.Iter(); // //Vector v = D.MulD(D.x); // //Text = Dirac.V1hermV2(D.ri, D.ri).ToString();//"Dx = " + Dirac.V1hermV2(v, v).Module() + " phi = " + philen.ToString(); // //Application.DoEvents(); // } // tr += Dirac.V1hermV2(D.phi, D.x); // Text = beta.ToString() + " " + k.ToString(); // Application.DoEvents(); // } // tr *= new Complex(1 / (double)K / (4 * 4 * 4 * 2 * 12 * 2), 0); // Text = tr.ToString(); // chart1.Series[0].Points.AddXY(beta, tr.x); //} }
private void GetData() { var l = new List<int>(new int[]{1,2,3}); //l.Where( DateTime T, Tfull; double Chi = 0; double BetaStart = Convert.ToDouble(tbBS.Text); double BetaEnd = Convert.ToDouble(tbBE.Text); double BetaStep = Convert.ToDouble(tbBSt.Text); int Tsteps = Convert.ToInt16(tbTsteps.Text); int Msteps = Convert.ToInt16(tbMsteps.Text); for (double b = BetaStart; b <= BetaEnd; b = b + BetaStep) { SU3 Field = new SU3(Convert.ToInt32(tbNs.Text), Convert.ToInt32(tbNt.Text), b,Convert.ToDouble(tbflux.Text), checkBox1.Checked); double PL = 0, PL2 = 0, PLt; Chi = 0; int n = 0; Tfull = DateTime.Now; InvokeMe(() => Text = "1"); InvokeMe(() => lblProcess.Text = "Thermalizing"); for (int i = 0; i < Tsteps; i++) { //T = DateTime.Now; Field.Sweep(false); InvokeMe(() => lblProcess.Text = "Thermalizing " + i.ToString()); // lblTime1.Invoke(new MethodInvoker(delegate() { lblTime1.Text = (DateTime.Now - T).TotalMilliseconds.ToString(); })); } //chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.AddXY(b, 0); })); for (int i = 0; i < Msteps; i++) { Field.Sweep(false); // for (int or = 0; or < 5; or++) Field.Sweep(false); if (rbS.Checked) PLt = Field.MeasureS(); else PLt = Field.MeasurePL(); PL += PLt; n++; PL2 += PLt * PLt; Chi = PL2 / n - (PL / n) * (PL / n); // chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.Last().YValues = new double[1] { Chi }; })); InvokeMe(() => chart1.Series[0].Points.Add(Chi)); InvokeMe(() => chart2.Series[0].Points.Add(PLt)); InvokeMe(() => lblProcess.Text = "Measuring " + i.ToString()); } // Chi = PL2 / n - (PL / n) * (PL / n); // chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.Last().YValues = new double[1] { Chi }; })); // chart2.Invoke(new MethodInvoker(delegate() { chart2.Series[0].Points.AddXY(b, PL/n); })); double t = (DateTime.Now - Tfull).TotalMilliseconds * (BetaEnd - b) / BetaStep; InvokeMe(() => lblRemain.Text = (DateTime.Now.AddMilliseconds(t)).ToString()); InvokeMe(() => chart2.Update()); } }
private void button14_Click(object sender, EventArgs e) { double beta = 6; SU3 Field = new SU3(int.Parse(tbNs.Text), int.Parse(tbNt.Text), beta, 0, true);//4х2 field beta=6 flux=0 DateTime t = DateTime.Now; for (int i = 0; i < 4; i++) { Field.Sweep(false); MessageBox.Show(Field.MeasureS().ToString()); Text = (DateTime.Now - t).TotalMilliseconds.ToString(); Application.DoEvents(); } }