如何从 python 调用八度函数



我有一个 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" 可以解决问题。

最新更新