我在同一台计算机上安装了 Mathematica 7.01 和 Mathematica 5.2。我希望能够在 Mathematica 7.01 会话中评估 v.5.2 内核中的代码。我的意思是,运行 Mathematica 7.0.1 标准会话时,我希望有一个像 kernel5Evaluate
这样的命令来评估 5.2 内核中的某些代码,并将结果返回到 7.01 内核并链接7.01 FrontEnd Notebook 中的代码将在 7.01 内核中执行。
例如(在标准 Mathematica v.7.01 会话中):
In[1]:= solutionFrom5 = kernel5Evaluate[NDSolve[{(y^\[Prime])[x]==y[x],y[1]==2},y,{x,0,3}]]
Out[1]= {{y -> InterpolatingFunction[{{0., 3.}}, <>]}}
In[2]:= kernel5Evaluate[Plot3D[Sin[x y],{x,-Pi,Pi},{y,-Pi,Pi}]]
During evaluation of In[2]:= GraphicsData["PostScript", "\<\............
Out[2]= -SurfaceGraphics-
在这两种情况下,结果应该就像 v.5.2 内核在 v.7.01 前端中设置为“笔记本内核”一样。当然 solutionFrom5
变量应该设置为 v.5.2 内核返回的真实解决方案。
I have Mathematica 7.01 and Mathematica 5.2 installed on the same machine. I wish to be able to evaluate code in the v.5.2 kernel from within Mathematica 7.01 session. I mean that running Mathematica 7.0.1 standard session I wish to have a command like kernel5Evaluate
to evaluate some code in the 5.2 kernel and return the result into the 7.01 kernel and linked 7.01 FrontEnd notebook in such a way as this code would be executed in the 7.01 kernel.
For example (in the standard Mathematica v.7.01 session):
In[1]:= solutionFrom5 = kernel5Evaluate[NDSolve[{(y^\[Prime])[x]==y[x],y[1]==2},y,{x,0,3}]]
Out[1]= {{y -> InterpolatingFunction[{{0., 3.}}, <>]}}
In[2]:= kernel5Evaluate[Plot3D[Sin[x y],{x,-Pi,Pi},{y,-Pi,Pi}]]
During evaluation of In[2]:= GraphicsData["PostScript", "\<\............
Out[2]= -SurfaceGraphics-
In the both cases the result should be as if the v.5.2 kernel is set to be "Notebook's Kernel" in the v.7.01 FrontEnd. And of course solutionFrom5
variable should be set to the real solution returned by v.5.2 kernel.
发布评论
评论(3)
这是基于 Simon 代码的实现。它仍然需要改进。我不清楚的一件事是如何处理从属(v.5.2)内核中生成的消息。
这是我的代码:
这是测试表达式:
它似乎按预期工作。不过,它可能会更好...
Here is an implementation based on Simon's code. It still requires improvement. The one unclear thing to me is how to handle Messages generated in the slave (v.5.2) kernel.
Here is my code:
Here are test expressions:
It seems to work as expected. However it could be better...
这是我想要的工作实现。我按照 Todd Gayley MathLink 连接的检查“nofollow”>
此处
。现在,即使从属内核以异常方式终止,kernel5Evaluate
也可以可靠地工作。我还大大改进了Message
的解析,并为kernel5Evaluate
添加了一些诊断消息。这是代码:这是一些测试表达式:
Here is working implementation of what I wanted. I have added checking for a dead
MathLink
connection as suggested by Todd Gayleyhere
. Nowkernel5Evaluate
works reliable even if the slave kernel was terminated in unusual way. I also have much improved parsing ofMessage
s and added some diagnostic messages forkernel5Evaluate
. Here is the code:And here are some test expressions:
这是我对您想要的尝试,
首先我定义了
linkEvaluate
,它接受一个活动的Link
并向其传递一个表达式。如果
LinkRead
仍有需要读取的内容,则会读取这些内容,直到没有更多内容为止。然后它写入表达式并等待结果返回。
然后它读取输出,直到没有任何内容可供读取。
通常,它会返回第一个
ReturnExpressionPacket
,除非您已将最后一个可选参数all
设置为True
- 在这种情况下,它会返回所有内容它读到。然后
kernel5Evaluate
首先检查全局$kern5
是否定义为LinkObject
,如果没有则定义它。然后,它只是将工作传递给linkEvaluate
。您必须将“math5”替换为 Mma 5.2 内核的文件名和路径。
Here's my attempt at what you want,
First I define
linkEvaluate
that takes an activeLink
and passes it an expression.If there's things for
LinkRead
still to read, then it reads them until there are no more.Then it writes the expression and waits for the results to come back.
Then it reads the output until there's nothing left to read.
Normally, it then returns the first
ReturnExpressionPacket
unless you have set the final, optional argument,all
, toTrue
- in which case it returns everything it read.Then
kernel5Evaluate
first checks if the global$kern5
is defined as aLinkObject
, if not then it defines it. It then simply passes the work over tolinkEvaluate
.You will have to replace "math5" with the filename and path of your Mma 5.2 kernel.