我已经编写或者尝试使用 ZedGraph 在 C# 中绘制可用空间薛定谔方程,并在运行创建的可执行文件时遇到 System.OutOfMemoryException 错误。
************** Exception Text **************
System.OutOfMemoryException: Exception of type 'System.OutOfMemoryException' was thrown.
at System.Collections.Generic.List`1.set_Capacity(Int32 value)
at System.Collections.Generic.List`1.EnsureCapacity(Int32 min)
at System.Collections.Generic.List`1.Add(T item)
at ZedGraph.PointPairList.Add(Double x, Double y)
at Wavepacket.Form1.CreateGraph(ZedGraphControl schrodinger) in C:\Users\user1748005\Documents\Visual Studio 2010\Projects\Wavepacket\Wavepacket\Form1.cs:line 72
at Wavepacket.Form1.Form1_Load(Object sender, EventArgs e) in C:\Users\user1748005\Documents\Visual Studio 2010\Projects\Wavepacket\Wavepacket\Form1.cs:line 35
at System.Windows.Forms.Form.OnLoad(EventArgs e)
at System.Windows.Forms.Form.OnCreateControl()
at System.Windows.Forms.Control.CreateControl(Boolean fIgnoreVisible)
at System.Windows.Forms.Control.CreateControl()
at System.Windows.Forms.Control.WmShowWindow(Message& m)
at System.Windows.Forms.Control.WndProc(Message& m)
at System.Windows.Forms.ScrollableControl.WndProc(Message& m)
at System.Windows.Forms.Form.WmShowWindow(Message& m)
at System.Windows.Forms.Form.WndProc(Message& m)
at System.Windows.Forms.Control.ControlNativeWindow.OnMessage(Message& m)
at System.Windows.Forms.Control.ControlNativeWindow.WndProc(Message& m)
at System.Windows.Forms.NativeWindow.Callback(IntPtr hWnd, Int32 msg, IntPtr wparam, IntPtr lparam)
我的 Cmplx 类的必要函数/构造函数相对简单,如下所示:
public Cmplx(double r, double i)
{
this.real = r;
this.imag = i;
}
public double Real
{
get { return this.real; }
set { this.real = value; }
}
public double Imag
{
get { return this.imag; }
set { this.imag = value; }
}
public Cmplx Multiply(Cmplx b)
{
return new Cmplx((this.real * b.real - this.imag * b.imag), (b.imag * this.real + this.imag * b.real));
}
public Cmplx Divide(Cmplx b)
{
return new Cmplx((this.real * b.real + this.imag * b.imag) / (b.real * b.real + b.imag * b.imag), (b.real * this.imag - this.real * b.imag) / (b.real * b.real + b.imag * b.imag));
}
public Cmplx Scale(Double b)
{
return new Cmplx(this.real * b, this.imag * b);
}
public double Mod()
{
return Math.Sqrt(this.real * this.real + this.imag * this.imag);
}
public Cmplx Sqrt()
{
double r, theta;
r = this.Mod();
theta = this.Arg();
return new Cmplx(Math.Sqrt(r) * (Math.Cos(theta / 2)), Math.Sqrt(r) * Math.Sin(theta / 2));
}
public Cmplx Exp()
{
double x, y, z;
x = Math.Exp(this.real);
y = Math.Cos(this.imag);
z = Math.Sin(this.imag);
return new Cmplx(x * y, x * z);
}
下面是我的 ZedGraph 方法,其中只有从 Codeproject ZedGraph 教程修改的 CreateGraph 方法:
public Form1()
{
InitializeComponent();
}
private void Form1_Resize(object sender, EventArgs e)
{
SetSize();
}
private void SetSize()
{
zedGraphControl1.Location = new Point(10, 10);
zedGraphControl1.Size = new Size(ClientRectangle.Width - 20,
ClientRectangle.Height - 20);
}
private void Form1_Load(object sender, EventArgs e)
{
CreateGraph(zedGraphControl1);
SetSize();
}
private void CreateGraph(ZedGraphControl schrodinger)
{
GraphPane myPane = schrodinger.GraphPane;
myPane.Title.Text = "Free Space Schrodinger Equation)";
myPane.XAxis.Title.Text = "Time (s)";
myPane.YAxis.Title.Text = "Psi";
double alpha, beta, norm, k0, t;
double L,N;
L = 25;
N = 4096;
alpha = 0.5;
beta = 0.5785;
k0 = 1;
t = 0;
norm = (Math.Exp(-k0 * k0 * alpha)) * (Math.Pow((alpha / (2 * Math.PI)), 0.25));
PointPairList list1 = new PointPairList();
for (double i = -L; i <= L; i -= (-2L) / N)
{
Cmplx gamma = new Cmplx(alpha, beta * t);
Cmplx a = new Cmplx(2 * alpha * k0, 0);
a.Imag = i;
Cmplx b = a.Multiply(a);
Cmplx phi = (((b.Divide(gamma.Scale(4))).Exp()).Scale(norm)).Divide(gamma.Sqrt());
list1.Add(phi.Mod(),i);
}
LineItem myCurve = myPane.AddCurve("t=0",
list1, Color.Red, SymbolType.Default);
schrodinger.AxisChange();
}
我改变了
list1.Add(phi.Mod(), i);
上述代码块的 for 循环中的行到
Console.WriteLine("{0}, {1}", phi.Mod(), i);
这可以按预期工作,尽管在命令行界面中。使用 ZedGraph 库结果的 Windows 窗体应用程序在引发上述异常之前消耗了大约 1.7GB 的 RAM。
我通过绘制大量数据点找到了 ZedGraph 内存不足的实例
但是,我不太倾向于我绘制的 4096 个数据点会导致这种内存消耗,而宁愿指责我的(缺乏)编程技能
注意:这是一项学术作业,因此提供的任何帮助都应考虑到这一点。干杯!
编辑:将 L 和 N 的初始化固定为整数而不是双精度数。如前所述,问题仍然存在。
EDIT2我是个白痴。我在我的测试 CLI 项目中将 L 和 N 的初始化更改为双精度,但未能更新 Windows 窗体应用程序的源代码。感谢 Jeppe Stig Nielsen 为我指明了正确的方向。