我有一个 problem.m 文件,它有一个这样的函数:myfun(p,m)
.
它执行一些计算并返回结果。为了测试此函数的执行,我有一个看起来像这样的 test.m 文件。
myfun(myarg1,myarg2)
如果我以以下方式运行此文件: octave test.m
然后它向我返回正确的结果,如下所示:0.38007
现在,问题是当使用python调用这个函数myfun(p,m)时。我尝试使用python库:oct2py
python代码如下所示:
import sys
import subprocess
import oct2py
import requests
from oct2py import octave
def myfun(p,m):
octave.addpath('/mypath');
oc = oct2py.Oct2Py()
#res = octave.myfun(p,m nout=1);#this didn't work, hence comment
#res = oc.myfun(p, m) #this didn't work, hence comment
x = oc.feval("myfun",p,m);
print(x);
if __name__ == '__main__':
myfun(sys.argv[1],sys.argv[2])
当我以:python FileName.py arg1 arg2(与我在test.m中使用的参数相同)运行此代码时,它给了我一条警告消息和一个空列表,如下所示:
警告:返回值列表中的某些元素未定义 []
我不确定该怎么办。由于该函数在使用八度音阶时似乎以正确的格式返回正确的结果。但由于某种原因,Oct2py无法正常工作。
八度代码:
function result = test(problem,method)
Methods={'MC','COS','RBF-AD','RBF-MLT'};
if problem == 1
S=[90,100,110]; K=100; T=1.0; r=0.03; sig=0.15;
U=[2.758443856146076 7.485087593912603 14.702019669720769];
rootpath=pwd;
filepathsBSeuCallUI=getfilenames('./','BSeuCallUI_*.m');
par={S,K,T,r,sig};
[timeBSeuCallUI,relerrBSeuCallUI] =
executor(rootpath,filepathsBSeuCallUI,U,par)
tBSeuCallUI=NaN(numel(Methods),1); rBSeuCallUI=tBSeuCallUI;
for ii=1:numel(Methods)
for jj=1:numel(filepathsBSeuCallUI)
a=filepathsBSeuCallUI{jj}(3:3+numel(Methods{ii}));
b=[Methods{ii},'/'];
if strcmp(a,b)
tBSeuCallUI(ii)=timeBSeuCallUI(jj);
rBSeuCallUI(ii)=relerrBSeuCallUI(jj);
end
end
end
cd(rootpath);
for el=1:numel(Methods)
if strcmp((Methods{el}),method)
result = tBSeuCallUI(el);
end
end
end
蟒蛇代码:
import sys
import subprocess
import oct2py
import requests
def benchop(problem,method):
oc = oct2py.Oct2Py()
res = oc.test(problem,method)
return res
if __name__ == '__main__':
benchop(sys.argv[1],sys.argv[2])
八度代码中的问题就像检查问题时一样,它说 problem==1,但 python 需要一个字符串,所以列表未定义的警告没有结果值。将八度代码从问题 == 1 更改为问题 == "1" 可以解决问题。