Estou tentando usarfsolve em combinação com ompmath pacote. No entanto, recebo o errocannot create mpf from array([mpf('1.0')], dtype=object).

Aqui está um exemplo mínimo de reprodução do erro. Neste exemplo, tecnicamente não preciso dompmath pacote, mas minha função real contém funções hiperconfluentes que o fazem.

from scipy.optimize import fsolve
#from mpmath import hyp1f1 as hyp1f1mp
#from mpmath import gamma as gammamp
import mpmath as mp
#import numpy as np

mp.dps = 250; mp.pretty = True


def cosFunc(p):
   vn = p
   output = mp.sin(vn)
   return output

estimate = mp.mpf(1)
value = fsolve(cosFunc,estimate)
print value

Encontrei uma pergunta semelhante sugerindo usarnp.frompyfunc (Como mpf uma matriz?), mas informa que a função não pode ser chamada (quando eu a aplico emvn)