def interface(inp, xdata, ydata):
p2 = _gauss2.gaussjac(xdata, ydata, inp)
#p = scipy.zeros((len(inp)+1,))
#p[0] = compare(inp,xdata,ydata)
#p[1:] = _approx_fprime_helper(inp,compare,1e-8,args=(xdata,ydata))
#print('new')
#print(inp)
#print(p)
#print(p2)
if scipy.any(scipy.isnan(p2)):
print(inp, p2)
return p2[0], p2[1:]
def interface(inp, xdata, ydata):
""" splits the c function output to two variables, the RMSE and the derivative of RMSE with respect to the parameters"""
p = _gauss.gaussjac(xdata,ydata,inp)
return p[0],p[1:]
def interface3(inp, xdata, ydata):
p = _gauss2.gaussjac(xdata, ydata, inp)
return p[0]