src/p/y/pysces-HEAD/trunk/pysces_pysundials/pysundials/ida.py pysces(Download)
def IDAGetWorkSpace(ida_mem):
lenrw = ctypes.c_long(0)
leniw = ctypes.c_long(0)
ret = ida.IDAGetWorkSpace(ida_mem.obj, ctypes.byref(lenrw), ctypes.byref(leniw))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetWorkSpace() failed with flag %i"%(ret))
return (lenrw.value, leniw.value)
ida.IDAGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDAGetWorkSpace.restype = ctypes.c_int
def IDAGetNumSteps(ida_mem):
nsteps = ctypes.c_long(0)
def IDAGetNumResEvals(ida_mem):
nrevals = ctypes.c_long(0)
ret = ida.IDAGetNumResEvals(ida_mem.obj, ctypes.byref(nrevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumResEvals() failed with flag %i"%(ret))
return nrevals.value
ida.IDAGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumResEvals.restype = ctypes.c_int
def IDAGetNumLinSolvSetups(ida_mem):
nlinsetups = ctypes.c_long(0)
def IDAGetNumErrTestFails(ida_mem):
netfails = ctypes.c_long(0)
ret = ida.IDAGetNumErrTestFails(ida_mem.obj, ctypes.byref(netfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumErrTestFails() failed with flag %i"%(ret))
return netfails.value
ida.IDAGetNumErrTestFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumErrTestFails.restype = ctypes.c_int
def IDAGetNumBacktrackOps(ida_mem):
nbacktr = ctypes.c_long(0)
def IDAGetNumGEvals(ida_mem):
ngevals = ctypes.c_long(0)
ret = ida.IDAGetNumGEvals(ida_mem.obj, ctypes.byref(ngevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumGEvals() failed with flag %i"%(ret))
return ngevals.value
ida.IDAGetNumGEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDAGetIntegratorStats(ida_mem): nsteps = ctypes.c_long(0) nrevals = ctypes.c_long(0) nlinsetups = ctypes.c_long(0) netfails = ctypes.c_long(0) qlast = ctypes.c_int(0) qcur = ctypes.c_int(0)
def IDAGetNumNonlinSolvIters(ida_mem):
nniters = ctypes.c_long(0)
ret = ida.IDAGetNumNonlinSolvIters(ida_mem.obj, ctypes.byref(nniters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumNonlinSolvIters() failed with flag %i"%(ret))
return nniters.value
ida.IDAGetNumNonlinSolvIters.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumNonlinSolvIters.restype = ctypes.c_int
def IDAGetNumNonlinSolvConvFails(ida_mem):
nncfails = ctypes.c_long(0)
def IDAGetNonlinSolvStats(ida_mem):
nniters = ctypes.c_long(0)
nncfails = ctypes.c_long(0)
ret = ida.IDAGetNonlinSolvStats(ida_mem.obj, ctypes.byref(nniters), ctypes.byref(nncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNonlinSolvStats() failed with flag %i"%(ret))
return (nniters.value, nncfils.value)
def IDABandGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDABandGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABandGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDABandGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDABandGetWorkSpace.restype = ctypes.c_int
def IDABandGetNumJacEvals(ida_mem):
njevals = ctypes.c_long(0)
def IDABandGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDABandGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABandGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDABandGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDABBDPrecGetWorkSpace(bbd_data):
lenrwBBDP = ctypes.c_long(0)
leniwBBDP = ctypes.c_long(0)
ret = ida.IDABBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwBBDP), ctypes.byref(leniwBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABBDPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwBBDP.value, leniwBBDP.value)
ida.IDABBDPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDABBDPrecGetWorkSpace.restype = ctypes.c_int
def IDABBDPrecGetNumGfnEvals(bbd_data):
ngevalsBBDP = ctypes.c_long(0)
def IDADenseGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDADenseGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDADenseGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDADenseGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDADenseGetWorkSpace.restype = ctypes.c_int
def IDADenseGetNumJacEvals(ida_mem):
njevals = ctypes.c_long(0)
def IDADenseGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDADenseGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDADenseGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDADenseGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDASpilsGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDASpilsGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDASpilsGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetWorkSpace.restype = ctypes.c_int
def IDASpilsGetNumPrecEvals(ida_mem):
npevals = ctypes.c_long(0)
def IDASpilsGetNumPrecSolves(ida_mem):
npsolves = ctypes.c_long(0)
ret = ida.IDASpilsGetNumPrecSolves(ida_mem.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumPrecSolves() failed with flag %i"%(ret))
return npsolves.value
ida.IDASpilsGetNumPrecSolves.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetNumPrecSolves.restype = ctypes.c_int
def IDASpilsGetNumLinIters(ida_mem):
nliters = ctypes.c_long(0)
def IDASpilsGetNumConvFails(ida_mem):
nlcfails = ctypes.c_long(0)
ret = ida.IDASpilsGetNumConvFails(ida_mem.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumConvFails() failed with flag %i"%(ret))
return nlcfails.value
ida.IDASpilsGetNumConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetNumConvFails.restype = ctypes.c_int
def IDASpilsGetNumJtimesEvals(ida_mem):
njvevals = ctypes.c_long(0)
def IDASpilsGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDASpilsGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDASpilsGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
src/p/y/pysces-HEAD/trunk/pysces_pysundials/pysundials/cvodes.py pysces(Download)
def CVodeGetWorkSpace(cvodememobj):
"""CVodeGetWorkSpace returns the CVODES real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrw = ctypes.c_long()
leniw = ctypes.c_long()
ret = cvodes.CVodeGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrw), ctypes.byref(leniw))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetWorkSpace() failed with flag %i"%(ret))
return (lenrw, leniw)
cvodes.CVodeGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetWorkSpace.restype = ctypes.c_int
def CVodeGetNumSteps(cvodememobj):
"""CVodeGetNumSteps returns the cumulative number of internal steps taken by the solver\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
def CVodeGetNumRhsEvals(cvodememobj):
"""CVodeGetNumRhsEvals returns the number of calls to the user's f function\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvals() failed with flag %i"%(ret))
return retval.value
cvodes.CVodeGetNumRhsEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumRhsEvals.restype = ctypes.c_int
def CVodeGetNumLinSolvSetups(cvodememobj):
"""CVodeGetNumLinSolvSetups returns the number of calls made to the linear solver's setup routine\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
def CVodeGetNumErrTestFails(cvodememobj):
"""CVodeGetNumErrTestFails returns the number of local error test failures that have occured\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumErrTestFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumErrTestFails() failed with flag %i"%(ret))
return retval.value
def CVodeGetNumStabLimOrderReds(cvodememobj):
"""CVodeGetNumStabLimOrderReds returns the number of order reductions due to stability limit detection\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumStabLimOrderReds(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStabLimOrderReds() failed with flag %i"%(ret))
return retval.value
def CVodeGetNumGEvals(cvodememobj):
"""CVodeGetNumGEvals returns the number of calls to the user's g function (for rootfinding)\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumGEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumGEvals() failed with flag %i"%(ret))
return retval.value
def CVodeGetIntegratorStats(cvodememobj): """A convenience function which returns a tuple of all available integrator stats.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()\n\nReturn tuple is of he format:\n \n\t(nsteps, nfevals, nlinsetups, netfails, qlast, qcur, hinused, hlast, hcur, tcur) where: \n\t\tnsteps: cumulative number of internal steps taken by the solver \n\t\tnfevals: number of rhs evaluations \n\t\tnlinsetups: number of calls to the linear setup function \n\t\tnetfails: number of local error test failures \n\t\tqlast: order used during ast internal step \n\t\tqcur: order to be used on next internal step \n\t\thinused: actual initial step size ised by CVODES \n\t\thlast: last step size used by CVODES \n\t\thcur: next step size to be used by CVODES \n\t\ttcur = current time""" nsteps = ctypes.c_long() nfevals = ctypes.c_long() nlinsetups = ctypes.c_long() netfails = ctypes.c_long() qlast = ctypes.c_int()
def CVodeGetNumNonlinSolvIters(cvodememobj):
"""CVodeGetNumNonlinSolvIters returns the number of nonlinear solver iterations performed.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumNonlinSolvIters(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvIters() failed with flag %i"%(ret))
return retval.value
cvodes.CVodeGetNumNonlinSolvIters.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumNonlinSolvIters.restype = ctypes.c_int
def CVodeGetNumNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumNonlinSolvConvFails returns the number of nonlinear convergence failures.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
def CVodeGetNonlinSolvStats(cvodememobj):
"""A convenience function that provides the nonlinear solver optional outputs in a tuple (NumNonlinSolvIters, NumNonlinsolvConvFails).\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nniters = ctypes.c_long()
nncfails = ctypes.c_long()
ret = cvodes.CVodeGetNonlinSolvStats(cvodememobj.obj, ctypes.byref(nniters), ctypes.byref(nncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNonlinSolvStats() failed with flag %i"%(ret))
def CVodeGetQuadNumRhsEvals(cvodememobj):
nfQevals = ctypes.c_long(0)
ret = cvodes.CVodeGetQuadNumRhsEvals(cvodememobj.obj, ctypes.byref(nfQevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetQuadNumRhsEvals() failed with flag %i"%(ret))
return nfQevals.value
cvodes.CVodeGetQuadNumRhsEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetQuadNumRhsEvals.restype = ctypes.c_int
def CVodeGetQuadNumErrTestFails(cvodememobj):
nQetfails = ctypes.c_long(0)
def CVodeGetQuadStats(cvodememobj):
nfQevals = ctypes.c_long(0)
nQetfails = ctypes.c_long(0)
ret = cvodes.CVodeGetQuadStats(cvodememobj.obj, ctypes.byref(nfQevals), ctypes.byref(nQetfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetQuadStats() failed with flag %i"%(ret))
return (nfQevals.value, nQetfails.value)
def CVodeGetNumSensRhsEvals(cvodememobj):
nfSevals = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensRhsEvals(cvodememobj.obj, ctypes.byref(nfSevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensRhsEvals() failed with flag %i"%(ret))
return nfSevals.value
cvodes.CVodeGetNumSensRhsEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumSensRhsEvals.restype = ctypes.c_int
def CVodeGetNumRhsEvalsSens(cvodememobj):
nfevalsS = ctypes.c_long(0)
def CVodeGetNumSensErrTestFails(cvodememobj):
nSetfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensErrTestFails(cvodememobj.obj, ctypes.byref(nSetfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensErrTestFails() failed with flag %i"%(ret))
return nSetfails.value
cvodes.CVodeGetNumSensErrTestFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumSensErrTestFails.restype = ctypes.c_int
def CVodeGetNumSensLinSolvSetups(cvodememobj):
nlinsetupsS = ctypes.c_long(0)
def CVodeGetSensStats(cvodememobj):
nfSevals = ctypes.c_long(0)
nfevalsS = ctypes.c_long(0)
nSetfails = ctypes.c_long(0)
nlinsetupS = ctypes.c_long(0)
ret = cvodes.CVodeGetSensStats(cvodememobj.obj, ctypes.byref(nfSevals), ctypes.byref(nfevalsS), ctypes.byref(nSetfails), ctypes.byref(nlinsetupsS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetSensStats() failed with flag %i"%(ret))
return (nfSevals.value, nfevalsS.value, nSetfails.value, nlinsetupS.value)
cvodes.CVodeGetSensStats.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetSensStats.restype = ctypes.c_int
def CVodeGetNumSensNonlinSolvIters(cvodememobj):
nSniters = ctypes.c_long(0)
def CVodeGetNumSensNonlinSolvConvFails(cvodememobj):
nSncfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(nSncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensNonlinSolvConvFails() failed with flag %i"%(ret))
return nSncfails.value
cvodes.CVodeGetNumSensNonlinSolvConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumSensNonlinSolvConvFails.restype = ctypes.c_int
def CVodeGetNumStgrSensNonlinSolvIters(cvodememobj):
nSTGR1niters = ctypes.c_long(0)
def CVodeGetNumStgrSensNonlinSolvConvFails(cvodememobj):
nSTGR1ncfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumStgrSensNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(nSTGR1ncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStgrSensNonlinSolvConvFails() failed with flag %i"%(ret))
return nSTGR1ncfails.value
cvodes.CVodeGetNumStgrSensNonlinSolvConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumStgrSensNonlinSolvConvFails.restype = ctypes.c_int
def CVodeGetSensNonlinSolvStats(cvodememobj):
nSniters = ctypes.c_long(0)
nSncfails = ctypes.c_long(0)
def CVDiagGetWorkSpace(cvodememobj):
"""CVDiagGetWorkSpace returns the CVDIAG real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVDiagGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDiagGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVDiagGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVDiagGetWorkSpace.restype = ctypes.c_int
def CVDiagGetNumRhsEvals(cvodememobj):
"""CVDiagGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian evaluation. Note: The number of diagonal approximate Jacobians formed is equal to the number of CVDiagSetup calls. This number is available through CVodeGetNumLinSolvSetups."""
retval = ctypes.c_long(0)
def CVBBDPrecGetWorkSpace(bbd_data):
"""CVBBDPrecGetWorkSpace returns the CVBBDPrec real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVBBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBBDPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVBBDPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVBBDPrecGetWorkSpace.restype = ctypes.c_int
def CVBBDPrecGetNumGfnEvals(bbd_data):
"""CVBBDPrecGetNumGfnEvals returns the number of calls to gfn."""
ngevalsBBDP = ctypes.c_long(0)
def CVBandPrecGetWorkSpace(bp_data):
"""CVBandPrecGetWorkSpace returns the CVBandPrec real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVBandPrecGetWorkSpace(bp_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVBandPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVBandPrecGetWorkSpace.restype = ctypes.c_int
def CVBandPrecGetNumRhsEvals(bp_data):
"""CVBandPrecGetNumGfnEvals returns the number of calls made from CVBANDPRE to the user's RHS."""
nfevalsBP = ctypes.c_long(0)
def CVBandGetWorkSpace(cvodememobj):
"""CVBandGetWorkSpace returns the CVBand real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long()
leniwLS = ctypes.c_long()
ret = cvodes.CVBandGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVBandGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVBandGetWorkSpace.restype = ctypes.c_int
def CVBandGetNumJacEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's jacobian function."""
retval = ctypes.c_long(0)
def CVBandGetNumRhsEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's RHS function."""
retval = ctypes.c_long(0)
ret = cvodes.CVBandGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumRhsEvals() failed with flag %i"%(ret))
return retval.value
def CVDenseGetWorkSpace(cvodememobj):
"""CVDenseGetWorkSpace returns the CVDense real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long()
leniwLS = ctypes.c_long()
ret = cvodes.CVDenseGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVDenseGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVDenseGetWorkSpace.restype = ctypes.c_int
def CVDenseGetNumJacEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's Jacobian function."""
retval = ctypes.c_long(0)
def CVDenseGetNumRhsEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's RHS function."""
retval = ctypes.c_long(0)
ret = cvodes.CVDenseGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumRhsEvals() failed with flag %i"%(ret))
return retval.value
def CVSpilsGetWorkSpace(cvodememobj):
"""CVSpilsGetWorkSpace returns the CVSpils real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVSpilsGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVSpilsGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVSpilsGetWorkSpace.restype = ctypes.c_int
def CVSpilsGetNumPrecEvals(cvodememobj):
"""CVSpilsGetNumPrecEvals returns the number of preconditioner evaluations, i.e. the number of calls made to PrecSetup with jok==FALSE.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
npevals = ctypes.c_long(0)
def CVSpilsGetNumPrecSolves(cvodememobj):
"""CVSpilsGetNumPrecSolves returns the number of calls made to PrecSolve.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
npsolves = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumPrecSolves(cvodememobj.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecSolves() failed with flag %i"%(ret))
return npsolves.value
cvodes.CVSpilsGetNumPrecSolves.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVSpilsGetNumPrecSolves.restype = ctypes.c_int
def CVSpilsGetNumLinIters(cvodememobj):
"""CVSpilsGetNumLinIters returns the number of linear iterations.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nliters = ctypes.c_long(0)
def CVSpilsGetNumConvFails(cvodememobj):
"""CVSpilsGetNumConvFails returns the number of linear convergence failures.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nlcfails = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumConvFails(cvodememobj.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumConvFails() failed with flag %i"%(ret))
return nlcfails.value
cvodes.CVSpilsGetNumConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVSpilsGetNumConvFails.restype = ctypes.c_int
def CVSpilsGetNumJtimesEvals(cvodememobj):
"""CVSpilsGetNumJtimesEvals returns the number of calls to jtimes.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
njvevals = ctypes.c_long(0)
def CVSpilsGetNumRhsEvals(cvodememobj):
"""CVSpilsGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian times vector evaluation.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nfevalsLS = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumRhsEvals(cvodememobj.obj, ctypes.byref(nfevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumRhsEvals() failed with flag %i"%(ret))
return nfevalsLS.value
src/p/y/pysces-HEAD/pysces_pysundials/pysundials/ida.py pysces(Download)
def IDAGetWorkSpace(ida_mem):
lenrw = ctypes.c_long(0)
leniw = ctypes.c_long(0)
ret = ida.IDAGetWorkSpace(ida_mem.obj, ctypes.byref(lenrw), ctypes.byref(leniw))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetWorkSpace() failed with flag %i"%(ret))
return (lenrw.value, leniw.value)
ida.IDAGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDAGetWorkSpace.restype = ctypes.c_int
def IDAGetNumSteps(ida_mem):
nsteps = ctypes.c_long(0)
def IDAGetNumResEvals(ida_mem):
nrevals = ctypes.c_long(0)
ret = ida.IDAGetNumResEvals(ida_mem.obj, ctypes.byref(nrevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumResEvals() failed with flag %i"%(ret))
return nrevals.value
ida.IDAGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumResEvals.restype = ctypes.c_int
def IDAGetNumLinSolvSetups(ida_mem):
nlinsetups = ctypes.c_long(0)
def IDAGetNumErrTestFails(ida_mem):
netfails = ctypes.c_long(0)
ret = ida.IDAGetNumErrTestFails(ida_mem.obj, ctypes.byref(netfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumErrTestFails() failed with flag %i"%(ret))
return netfails.value
ida.IDAGetNumErrTestFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumErrTestFails.restype = ctypes.c_int
def IDAGetNumBacktrackOps(ida_mem):
nbacktr = ctypes.c_long(0)
def IDAGetNumGEvals(ida_mem):
ngevals = ctypes.c_long(0)
ret = ida.IDAGetNumGEvals(ida_mem.obj, ctypes.byref(ngevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumGEvals() failed with flag %i"%(ret))
return ngevals.value
ida.IDAGetNumGEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDAGetIntegratorStats(ida_mem): nsteps = ctypes.c_long(0) nrevals = ctypes.c_long(0) nlinsetups = ctypes.c_long(0) netfails = ctypes.c_long(0) qlast = ctypes.c_int(0) qcur = ctypes.c_int(0)
def IDAGetNumNonlinSolvIters(ida_mem):
nniters = ctypes.c_long(0)
ret = ida.IDAGetNumNonlinSolvIters(ida_mem.obj, ctypes.byref(nniters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumNonlinSolvIters() failed with flag %i"%(ret))
return nniters.value
ida.IDAGetNumNonlinSolvIters.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumNonlinSolvIters.restype = ctypes.c_int
def IDAGetNumNonlinSolvConvFails(ida_mem):
nncfails = ctypes.c_long(0)
def IDAGetNonlinSolvStats(ida_mem):
nniters = ctypes.c_long(0)
nncfails = ctypes.c_long(0)
ret = ida.IDAGetNonlinSolvStats(ida_mem.obj, ctypes.byref(nniters), ctypes.byref(nncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNonlinSolvStats() failed with flag %i"%(ret))
return (nniters.value, nncfils.value)
def IDABandGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDABandGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABandGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDABandGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDABandGetWorkSpace.restype = ctypes.c_int
def IDABandGetNumJacEvals(ida_mem):
njevals = ctypes.c_long(0)
def IDABandGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDABandGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABandGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDABandGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDABBDPrecGetWorkSpace(bbd_data):
lenrwBBDP = ctypes.c_long(0)
leniwBBDP = ctypes.c_long(0)
ret = ida.IDABBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwBBDP), ctypes.byref(leniwBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABBDPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwBBDP.value, leniwBBDP.value)
ida.IDABBDPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDABBDPrecGetWorkSpace.restype = ctypes.c_int
def IDABBDPrecGetNumGfnEvals(bbd_data):
ngevalsBBDP = ctypes.c_long(0)
def IDADenseGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDADenseGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDADenseGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDADenseGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDADenseGetWorkSpace.restype = ctypes.c_int
def IDADenseGetNumJacEvals(ida_mem):
njevals = ctypes.c_long(0)
def IDADenseGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDADenseGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDADenseGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDADenseGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDASpilsGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDASpilsGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDASpilsGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetWorkSpace.restype = ctypes.c_int
def IDASpilsGetNumPrecEvals(ida_mem):
npevals = ctypes.c_long(0)
def IDASpilsGetNumPrecSolves(ida_mem):
npsolves = ctypes.c_long(0)
ret = ida.IDASpilsGetNumPrecSolves(ida_mem.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumPrecSolves() failed with flag %i"%(ret))
return npsolves.value
ida.IDASpilsGetNumPrecSolves.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetNumPrecSolves.restype = ctypes.c_int
def IDASpilsGetNumLinIters(ida_mem):
nliters = ctypes.c_long(0)
def IDASpilsGetNumConvFails(ida_mem):
nlcfails = ctypes.c_long(0)
ret = ida.IDASpilsGetNumConvFails(ida_mem.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumConvFails() failed with flag %i"%(ret))
return nlcfails.value
ida.IDASpilsGetNumConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetNumConvFails.restype = ctypes.c_int
def IDASpilsGetNumJtimesEvals(ida_mem):
njvevals = ctypes.c_long(0)
def IDASpilsGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDASpilsGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDASpilsGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
src/p/y/pysces-HEAD/pysces_pysundials/pysundials/cvodes.py pysces(Download)
def CVodeGetWorkSpace(cvodememobj):
"""CVodeGetWorkSpace returns the CVODES real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrw = ctypes.c_long()
leniw = ctypes.c_long()
ret = cvodes.CVodeGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrw), ctypes.byref(leniw))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetWorkSpace() failed with flag %i"%(ret))
return (lenrw, leniw)
cvodes.CVodeGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetWorkSpace.restype = ctypes.c_int
def CVodeGetNumSteps(cvodememobj):
"""CVodeGetNumSteps returns the cumulative number of internal steps taken by the solver\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
def CVodeGetNumRhsEvals(cvodememobj):
"""CVodeGetNumRhsEvals returns the number of calls to the user's f function\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvals() failed with flag %i"%(ret))
return retval.value
cvodes.CVodeGetNumRhsEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumRhsEvals.restype = ctypes.c_int
def CVodeGetNumLinSolvSetups(cvodememobj):
"""CVodeGetNumLinSolvSetups returns the number of calls made to the linear solver's setup routine\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
def CVodeGetNumErrTestFails(cvodememobj):
"""CVodeGetNumErrTestFails returns the number of local error test failures that have occured\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumErrTestFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumErrTestFails() failed with flag %i"%(ret))
return retval.value
def CVodeGetNumStabLimOrderReds(cvodememobj):
"""CVodeGetNumStabLimOrderReds returns the number of order reductions due to stability limit detection\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumStabLimOrderReds(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStabLimOrderReds() failed with flag %i"%(ret))
return retval.value
def CVodeGetNumGEvals(cvodememobj):
"""CVodeGetNumGEvals returns the number of calls to the user's g function (for rootfinding)\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumGEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumGEvals() failed with flag %i"%(ret))
return retval.value
def CVodeGetIntegratorStats(cvodememobj): """A convenience function which returns a tuple of all available integrator stats.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()\n\nReturn tuple is of he format:\n \n\t(nsteps, nfevals, nlinsetups, netfails, qlast, qcur, hinused, hlast, hcur, tcur) where: \n\t\tnsteps: cumulative number of internal steps taken by the solver \n\t\tnfevals: number of rhs evaluations \n\t\tnlinsetups: number of calls to the linear setup function \n\t\tnetfails: number of local error test failures \n\t\tqlast: order used during ast internal step \n\t\tqcur: order to be used on next internal step \n\t\thinused: actual initial step size ised by CVODES \n\t\thlast: last step size used by CVODES \n\t\thcur: next step size to be used by CVODES \n\t\ttcur = current time""" nsteps = ctypes.c_long() nfevals = ctypes.c_long() nlinsetups = ctypes.c_long() netfails = ctypes.c_long() qlast = ctypes.c_int()
def CVodeGetNumNonlinSolvIters(cvodememobj):
"""CVodeGetNumNonlinSolvIters returns the number of nonlinear solver iterations performed.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumNonlinSolvIters(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvIters() failed with flag %i"%(ret))
return retval.value
cvodes.CVodeGetNumNonlinSolvIters.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumNonlinSolvIters.restype = ctypes.c_int
def CVodeGetNumNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumNonlinSolvConvFails returns the number of nonlinear convergence failures.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
def CVodeGetNonlinSolvStats(cvodememobj):
"""A convenience function that provides the nonlinear solver optional outputs in a tuple (NumNonlinSolvIters, NumNonlinsolvConvFails).\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nniters = ctypes.c_long()
nncfails = ctypes.c_long()
ret = cvodes.CVodeGetNonlinSolvStats(cvodememobj.obj, ctypes.byref(nniters), ctypes.byref(nncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNonlinSolvStats() failed with flag %i"%(ret))
def CVodeGetQuadNumRhsEvals(cvodememobj):
nfQevals = ctypes.c_long(0)
ret = cvodes.CVodeGetQuadNumRhsEvals(cvodememobj.obj, ctypes.byref(nfQevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetQuadNumRhsEvals() failed with flag %i"%(ret))
return nfQevals.value
cvodes.CVodeGetQuadNumRhsEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetQuadNumRhsEvals.restype = ctypes.c_int
def CVodeGetQuadNumErrTestFails(cvodememobj):
nQetfails = ctypes.c_long(0)
def CVodeGetQuadStats(cvodememobj):
nfQevals = ctypes.c_long(0)
nQetfails = ctypes.c_long(0)
ret = cvodes.CVodeGetQuadStats(cvodememobj.obj, ctypes.byref(nfQevals), ctypes.byref(nQetfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetQuadStats() failed with flag %i"%(ret))
return (nfQevals.value, nQetfails.value)
def CVodeGetNumSensRhsEvals(cvodememobj):
nfSevals = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensRhsEvals(cvodememobj.obj, ctypes.byref(nfSevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensRhsEvals() failed with flag %i"%(ret))
return nfSevals.value
cvodes.CVodeGetNumSensRhsEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumSensRhsEvals.restype = ctypes.c_int
def CVodeGetNumRhsEvalsSens(cvodememobj):
nfevalsS = ctypes.c_long(0)
def CVodeGetNumSensErrTestFails(cvodememobj):
nSetfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensErrTestFails(cvodememobj.obj, ctypes.byref(nSetfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensErrTestFails() failed with flag %i"%(ret))
return nSetfails.value
cvodes.CVodeGetNumSensErrTestFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumSensErrTestFails.restype = ctypes.c_int
def CVodeGetNumSensLinSolvSetups(cvodememobj):
nlinsetupsS = ctypes.c_long(0)
def CVodeGetSensStats(cvodememobj):
nfSevals = ctypes.c_long(0)
nfevalsS = ctypes.c_long(0)
nSetfails = ctypes.c_long(0)
nlinsetupS = ctypes.c_long(0)
ret = cvodes.CVodeGetSensStats(cvodememobj.obj, ctypes.byref(nfSevals), ctypes.byref(nfevalsS), ctypes.byref(nSetfails), ctypes.byref(nlinsetupsS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetSensStats() failed with flag %i"%(ret))
return (nfSevals.value, nfevalsS.value, nSetfails.value, nlinsetupS.value)
cvodes.CVodeGetSensStats.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetSensStats.restype = ctypes.c_int
def CVodeGetNumSensNonlinSolvIters(cvodememobj):
nSniters = ctypes.c_long(0)
def CVodeGetNumSensNonlinSolvConvFails(cvodememobj):
nSncfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(nSncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensNonlinSolvConvFails() failed with flag %i"%(ret))
return nSncfails.value
cvodes.CVodeGetNumSensNonlinSolvConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumSensNonlinSolvConvFails.restype = ctypes.c_int
def CVodeGetNumStgrSensNonlinSolvIters(cvodememobj):
nSTGR1niters = ctypes.c_long(0)
def CVodeGetNumStgrSensNonlinSolvConvFails(cvodememobj):
nSTGR1ncfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumStgrSensNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(nSTGR1ncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStgrSensNonlinSolvConvFails() failed with flag %i"%(ret))
return nSTGR1ncfails.value
cvodes.CVodeGetNumStgrSensNonlinSolvConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVodeGetNumStgrSensNonlinSolvConvFails.restype = ctypes.c_int
def CVodeGetSensNonlinSolvStats(cvodememobj):
nSniters = ctypes.c_long(0)
nSncfails = ctypes.c_long(0)
def CVDiagGetWorkSpace(cvodememobj):
"""CVDiagGetWorkSpace returns the CVDIAG real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVDiagGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDiagGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVDiagGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVDiagGetWorkSpace.restype = ctypes.c_int
def CVDiagGetNumRhsEvals(cvodememobj):
"""CVDiagGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian evaluation. Note: The number of diagonal approximate Jacobians formed is equal to the number of CVDiagSetup calls. This number is available through CVodeGetNumLinSolvSetups."""
retval = ctypes.c_long(0)
def CVBBDPrecGetWorkSpace(bbd_data):
"""CVBBDPrecGetWorkSpace returns the CVBBDPrec real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVBBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBBDPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVBBDPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVBBDPrecGetWorkSpace.restype = ctypes.c_int
def CVBBDPrecGetNumGfnEvals(bbd_data):
"""CVBBDPrecGetNumGfnEvals returns the number of calls to gfn."""
ngevalsBBDP = ctypes.c_long(0)
def CVBandPrecGetWorkSpace(bp_data):
"""CVBandPrecGetWorkSpace returns the CVBandPrec real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVBandPrecGetWorkSpace(bp_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVBandPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVBandPrecGetWorkSpace.restype = ctypes.c_int
def CVBandPrecGetNumRhsEvals(bp_data):
"""CVBandPrecGetNumGfnEvals returns the number of calls made from CVBANDPRE to the user's RHS."""
nfevalsBP = ctypes.c_long(0)
def CVBandGetWorkSpace(cvodememobj):
"""CVBandGetWorkSpace returns the CVBand real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long()
leniwLS = ctypes.c_long()
ret = cvodes.CVBandGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVBandGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVBandGetWorkSpace.restype = ctypes.c_int
def CVBandGetNumJacEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's jacobian function."""
retval = ctypes.c_long(0)
def CVBandGetNumRhsEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's RHS function."""
retval = ctypes.c_long(0)
ret = cvodes.CVBandGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumRhsEvals() failed with flag %i"%(ret))
return retval.value
def CVDenseGetWorkSpace(cvodememobj):
"""CVDenseGetWorkSpace returns the CVDense real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long()
leniwLS = ctypes.c_long()
ret = cvodes.CVDenseGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVDenseGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVDenseGetWorkSpace.restype = ctypes.c_int
def CVDenseGetNumJacEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's Jacobian function."""
retval = ctypes.c_long(0)
def CVDenseGetNumRhsEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's RHS function."""
retval = ctypes.c_long(0)
ret = cvodes.CVDenseGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumRhsEvals() failed with flag %i"%(ret))
return retval.value
def CVSpilsGetWorkSpace(cvodememobj):
"""CVSpilsGetWorkSpace returns the CVSpils real and integer workspaces as a tuple (integer, real)\n\tcvodememobj\ta CVodeMemObj as returned by CVodeCreate()"""
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = cvodes.CVSpilsGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
cvodes.CVSpilsGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
cvodes.CVSpilsGetWorkSpace.restype = ctypes.c_int
def CVSpilsGetNumPrecEvals(cvodememobj):
"""CVSpilsGetNumPrecEvals returns the number of preconditioner evaluations, i.e. the number of calls made to PrecSetup with jok==FALSE.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
npevals = ctypes.c_long(0)
def CVSpilsGetNumPrecSolves(cvodememobj):
"""CVSpilsGetNumPrecSolves returns the number of calls made to PrecSolve.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
npsolves = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumPrecSolves(cvodememobj.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecSolves() failed with flag %i"%(ret))
return npsolves.value
cvodes.CVSpilsGetNumPrecSolves.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVSpilsGetNumPrecSolves.restype = ctypes.c_int
def CVSpilsGetNumLinIters(cvodememobj):
"""CVSpilsGetNumLinIters returns the number of linear iterations.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nliters = ctypes.c_long(0)
def CVSpilsGetNumConvFails(cvodememobj):
"""CVSpilsGetNumConvFails returns the number of linear convergence failures.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nlcfails = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumConvFails(cvodememobj.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumConvFails() failed with flag %i"%(ret))
return nlcfails.value
cvodes.CVSpilsGetNumConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvodes.CVSpilsGetNumConvFails.restype = ctypes.c_int
def CVSpilsGetNumJtimesEvals(cvodememobj):
"""CVSpilsGetNumJtimesEvals returns the number of calls to jtimes.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
njvevals = ctypes.c_long(0)
def CVSpilsGetNumRhsEvals(cvodememobj):
"""CVSpilsGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian times vector evaluation.\n\tcvodememobj\ta CvodeMemObj as returned by CvodeCreate()"""
nfevalsLS = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumRhsEvals(cvodememobj.obj, ctypes.byref(nfevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumRhsEvals() failed with flag %i"%(ret))
return nfevalsLS.value
src/p/y/pysces-HEAD/trunk/pysces_pysundials/pysundials/cvode.py pysces(Download)
def CVodeGetWorkSpace(cvodememobj): """CVodeGetWorkSpace returns the CVODE real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrw = ctypes.c_long() leniw = ctypes.c_long() ret = cvode.CVodeGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrw), ctypes.byref(leniw)) if ret < 0:
def CVodeGetNumSteps(cvodememobj):
"""CVodeGetNumSteps returns the cumulative number of internal steps taken by the solver
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumSteps(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSteps() failed with flag %i"%(ret))
def CVodeGetNumRhsEvals(cvodememobj):
"""CVodeGetNumRhsEvals returns the number of calls to the user's f function
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvals() failed with flag %i"%(ret))
def CVodeGetNumLinSolvSetups(cvodememobj):
"""CVodeGetNumLinSolvSetups returns the number of calls made to the linear solver's setup routine
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumLinSolvSetups(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumLinSolvSetups() failed with flag %i"%(ret))
def CVodeGetNumErrTestFails(cvodememobj):
"""CVodeGetNumErrTestFails returns the number of local error test failures that have occured
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumErrTestFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumErrTestFails() failed with flag %i"%(ret))
def CVodeGetNumStabLimOrderReds(cvodememobj):
"""CVodeGetNumStabLimOrderReds returns the number of order reductions due to stability limit detection
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumStabLimOrderReds(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStabLimOrderReds() failed with flag %i"%(ret))
def CVodeGetNumGEvals(cvodememobj):
"""CVodeGetNumGEvals returns the number of calls to the user's g function (for rootfinding)
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumGEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumGEvals() failed with flag %i"%(ret))
hlast last step size used by CVODE hcur next step size to be used by CVODE tcur current time""" nsteps = ctypes.c_long() nfevals = ctypes.c_long() nlinsetups = ctypes.c_long() netfails = ctypes.c_long()
def CVodeGetNumNonlinSolvIters(cvodememobj):
"""CVodeGetNumNonlinSolvIters returns the number of nonlinear solver iterations performed.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumNonlinSolvIters(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvIters() failed with flag %i"%(ret))
def CVodeGetNumNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumNonlinSolvConvFails returns the number of nonlinear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvConvFails() failed with flag %i"%(ret))
def CVodeGetNonlinSolvStats(cvodememobj): """A convenience function that provides the nonlinear solver optional outputs in a tuple (NumNonlinSolvIters, NumNonlinsolvConvFails). cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nniters = ctypes.c_long() nncfails = ctypes.c_long() ret = cvode.CVodeGetNonlinSolvStats(cvodememobj.obj, ctypes.byref(nniters), ctypes.byref(nncfails)) if ret < 0:
def CVDiagGetWorkSpace(cvodememobj): """Returns the CVDIAG real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVDiagGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDiagGetNumRhsEvals(cvodememobj): """Returns the number of calls to the user f routine due to finite difference Jacobian evaluation. cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate() Note: The number of diagonal approximate Jacobians formed is equal to the number of CVDiagSetup calls. This number is available through CVodeGetNumLinSolvSetups.""" retval = ctypes.c_long(0) ret = cvode.CVDiagGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
def CVBBDPrecGetWorkSpace(bbd_data): """CVBBDPrecGetWorkSpace returns the CVBBDPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVBBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBBDPrecGetNumGfnEvals(bbd_data):
"""CVBBDPrecGetNumGfnEvals returns the number of calls to gfn."""
ngevalsBBDP = ctypes.c_long(0)
ret = cvode.CVBBDPrecGetNumGfnEvals(bbd_data, ctypes.byref(ngevalsBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBBDPrecGetNumGfnEvals() failed with flag %i"%(ret))
return ngevalsBBDP.value
def CVBandPrecGetWorkSpace(bp_data): """CVBandPrecGetWorkSpace returns the CVBandPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVBandPrecGetWorkSpace(bp_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandPrecGetNumRhsEvals(bp_data):
"""CVBandPrecGetNumGfnEvals returns the number of calls made from CVBANDPRE to the user's RHS."""
nfevalsBP = ctypes.c_long(0)
ret = cvode.CVBandPrecGetNumRhsEvals(bp_data, ctypes.byref(nfevalsBP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandPrecGetNumRhsEvals() failed with flag %i"%(ret))
return nfevalsBP.value
def CVBandGetWorkSpace(cvodememobj): """CVBandGetWorkSpace returns the CVBand real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvode.CVBandGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandGetNumJacEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's jacobian function."""
retval = ctypes.c_long(0)
ret = cvode.CVBandGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumJacEvals() failed with flag %i"%(ret))
return retval.value
cvode.CVBandGetNumJacEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvode.CVBandGetNumJacEvals.restype = ctypes.c_int
def CVBandGetNumRhsEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's RHS function."""
retval = ctypes.c_long(0)
def CVDenseGetWorkSpace(cvodememobj): """CVDenseGetWorkSpace returns the CVDense real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvode.CVDenseGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDenseGetNumJacEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's Jacobian function."""
retval = ctypes.c_long(0)
ret = cvode.CVDenseGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumJacEvals() failed with flag %i"%(ret))
return retval.value
cvode.CVDenseGetNumJacEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvode.CVDenseGetNumJacEvals.restype = ctypes.c_int
def CVDenseGetNumRhsEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's RHS function."""
retval = ctypes.c_long(0)
def CVSpilsGetWorkSpace(cvodememobj): """CVSpilsGetWorkSpace returns the CVSpils real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVSpilsGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVSpilsGetNumPrecEvals(cvodememobj):
"""CVSpilsGetNumPrecEvals returns the number of preconditioner evaluations, i.e. the number of calls made to PrecSetup with jok==FALSE.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npevals = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumPrecEvals(cvodememobj.obj, ctypes.byref(npevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecEvals() failed with flag %i"%(ret))
def CVSpilsGetNumPrecSolves(cvodememobj):
"""CVSpilsGetNumPrecSolves returns the number of calls made to PrecSolve.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npsolves = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumPrecSolves(cvodememobj.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecSolves() failed with flag %i"%(ret))
def CVSpilsGetNumLinIters(cvodememobj):
"""CVSpilsGetNumLinIters returns the number of linear iterations.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nliters = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumLinIters(cvodememobj.obj, ctypes.byref(nliters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumLinIters() failed with flag %i"%(ret))
def CVSpilsGetNumConvFails(cvodememobj):
"""CVSpilsGetNumConvFails returns the number of linear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nlcfails = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumConvFails(cvodememobj.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumConvFails() failed with flag %i"%(ret))
def CVSpilsGetNumJtimesEvals(cvodememobj):
"""CVSpilsGetNumJtimesEvals returns the number of calls to jtimes.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
njvevals = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumJtimesEvals(cvodememobj.obj, ctypes.byref(njvevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumJtimesEvals() failed with flag %i"%(ret))
def CVSpilsGetNumRhsEvals(cvodememobj):
"""CVSpilsGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian times vector evaluation.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfevalsLS = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumRhsEvals(cvodememobj.obj, ctypes.byref(nfevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumRhsEvals() failed with flag %i"%(ret))
src/p/y/pysundials-HEAD/2.3.0/src/ida.py pysundials(Download)
def IDAGetWorkSpace(ida_mem):
lenrw = ctypes.c_long(0)
leniw = ctypes.c_long(0)
ret = ida.IDAGetWorkSpace(ida_mem.obj, ctypes.byref(lenrw), ctypes.byref(leniw))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetWorkSpace() failed with flag %i"%(ret))
return (lenrw.value, leniw.value)
ida.IDAGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDAGetWorkSpace.restype = ctypes.c_int
def IDAGetNumSteps(ida_mem):
nsteps = ctypes.c_long(0)
def IDAGetNumResEvals(ida_mem):
nrevals = ctypes.c_long(0)
ret = ida.IDAGetNumResEvals(ida_mem.obj, ctypes.byref(nrevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumResEvals() failed with flag %i"%(ret))
return nrevals.value
ida.IDAGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumResEvals.restype = ctypes.c_int
def IDAGetNumLinSolvSetups(ida_mem):
nlinsetups = ctypes.c_long(0)
def IDAGetNumErrTestFails(ida_mem):
netfails = ctypes.c_long(0)
ret = ida.IDAGetNumErrTestFails(ida_mem.obj, ctypes.byref(netfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumErrTestFails() failed with flag %i"%(ret))
return netfails.value
ida.IDAGetNumErrTestFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumErrTestFails.restype = ctypes.c_int
def IDAGetNumBacktrackOps(ida_mem):
nbacktr = ctypes.c_long(0)
def IDAGetNumGEvals(ida_mem):
ngevals = ctypes.c_long(0)
ret = ida.IDAGetNumGEvals(ida_mem.obj, ctypes.byref(ngevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumGEvals() failed with flag %i"%(ret))
return ngevals.value
ida.IDAGetNumGEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDAGetIntegratorStats(ida_mem): nsteps = ctypes.c_long(0) nrevals = ctypes.c_long(0) nlinsetups = ctypes.c_long(0) netfails = ctypes.c_long(0) qlast = ctypes.c_int(0) qcur = ctypes.c_int(0)
def IDAGetNumNonlinSolvIters(ida_mem):
nniters = ctypes.c_long(0)
ret = ida.IDAGetNumNonlinSolvIters(ida_mem.obj, ctypes.byref(nniters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNumNonlinSolvIters() failed with flag %i"%(ret))
return nniters.value
ida.IDAGetNumNonlinSolvIters.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDAGetNumNonlinSolvIters.restype = ctypes.c_int
def IDAGetNumNonlinSolvConvFails(ida_mem):
nncfails = ctypes.c_long(0)
def IDAGetNonlinSolvStats(ida_mem):
nniters = ctypes.c_long(0)
nncfails = ctypes.c_long(0)
ret = ida.IDAGetNonlinSolvStats(ida_mem.obj, ctypes.byref(nniters), ctypes.byref(nncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDAGetNonlinSolvStats() failed with flag %i"%(ret))
return (nniters.value, nncfils.value)
def IDABandGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDABandGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABandGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDABandGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDABandGetWorkSpace.restype = ctypes.c_int
def IDABandGetNumJacEvals(ida_mem):
njevals = ctypes.c_long(0)
def IDABandGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDABandGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABandGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDABandGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDABBDPrecGetWorkSpace(bbd_data):
lenrwBBDP = ctypes.c_long(0)
leniwBBDP = ctypes.c_long(0)
ret = ida.IDABBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwBBDP), ctypes.byref(leniwBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDABBDPrecGetWorkSpace() failed with flag %i"%(ret))
return (lenrwBBDP.value, leniwBBDP.value)
ida.IDABBDPrecGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDABBDPrecGetWorkSpace.restype = ctypes.c_int
def IDABBDPrecGetNumGfnEvals(bbd_data):
ngevalsBBDP = ctypes.c_long(0)
def IDADenseGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDADenseGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDADenseGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDADenseGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDADenseGetWorkSpace.restype = ctypes.c_int
def IDADenseGetNumJacEvals(ida_mem):
njevals = ctypes.c_long(0)
def IDADenseGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDADenseGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDADenseGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDADenseGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
def IDASpilsGetWorkSpace(ida_mem):
lenrwLS = ctypes.c_long(0)
leniwLS = ctypes.c_long(0)
ret = ida.IDASpilsGetWorkSpace(ida_mem.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetWorkSpace() failed with flag %i"%(ret))
return (lenrwLS.value, leniwLS.value)
ida.IDASpilsGetWorkSpace.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long), ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetWorkSpace.restype = ctypes.c_int
def IDASpilsGetNumPrecEvals(ida_mem):
npevals = ctypes.c_long(0)
def IDASpilsGetNumPrecSolves(ida_mem):
npsolves = ctypes.c_long(0)
ret = ida.IDASpilsGetNumPrecSolves(ida_mem.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumPrecSolves() failed with flag %i"%(ret))
return npsolves.value
ida.IDASpilsGetNumPrecSolves.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetNumPrecSolves.restype = ctypes.c_int
def IDASpilsGetNumLinIters(ida_mem):
nliters = ctypes.c_long(0)
def IDASpilsGetNumConvFails(ida_mem):
nlcfails = ctypes.c_long(0)
ret = ida.IDASpilsGetNumConvFails(ida_mem.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumConvFails() failed with flag %i"%(ret))
return nlcfails.value
ida.IDASpilsGetNumConvFails.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
ida.IDASpilsGetNumConvFails.restype = ctypes.c_int
def IDASpilsGetNumJtimesEvals(ida_mem):
njvevals = ctypes.c_long(0)
def IDASpilsGetNumResEvals(ida_mem):
nrevalsLS = ctypes.c_long(0)
ret = ida.IDASpilsGetNumResEvals(ida_mem.obj, ctypes.byref(nrevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: IDASpilsGetNumResEvals() failed with flag %i"%(ret))
return nrevalsLS.value
ida.IDASpilsGetNumResEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
src/p/y/pysundials-HEAD/2.3.0/src/cvodes.py pysundials(Download)
def CVodeGetWorkSpace(cvodememobj): """CVodeGetWorkSpace returns the CVODE real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrw = ctypes.c_long() leniw = ctypes.c_long() ret = cvodes.CVodeGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrw), ctypes.byref(leniw)) if ret < 0:
def CVodeGetNumSteps(cvodememobj):
"""CVodeGetNumSteps returns the cumulative number of internal steps taken by the solver
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumSteps(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSteps() failed with flag %i"%(ret))
def CVodeGetNumRhsEvals(cvodememobj):
"""CVodeGetNumRhsEvals returns the number of calls to the user's f function
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvals() failed with flag %i"%(ret))
def CVodeGetNumLinSolvSetups(cvodememobj):
"""CVodeGetNumLinSolvSetups returns the number of calls made to the linear solver's setup routine
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumLinSolvSetups(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumLinSolvSetups() failed with flag %i"%(ret))
def CVodeGetNumErrTestFails(cvodememobj):
"""CVodeGetNumErrTestFails returns the number of local error test failures that have occured
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumErrTestFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumErrTestFails() failed with flag %i"%(ret))
def CVodeGetNumStabLimOrderReds(cvodememobj):
"""CVodeGetNumStabLimOrderReds returns the number of order reductions due to stability limit detection
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvodes.CVodeGetNumStabLimOrderReds(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStabLimOrderReds() failed with flag %i"%(ret))
def CVodeGetNumGEvals(cvodememobj):
"""CVodeGetNumGEvals returns the number of calls to the user's g function (for rootfinding)
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumGEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumGEvals() failed with flag %i"%(ret))
hlast last step size used by CVODE hcur next step size to be used by CVODE tcur current time""" nsteps = ctypes.c_long() nfevals = ctypes.c_long() nlinsetups = ctypes.c_long() netfails = ctypes.c_long()
def CVodeGetNumNonlinSolvIters(cvodememobj):
"""CVodeGetNumNonlinSolvIters returns the number of nonlinear solver iterations performed.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumNonlinSolvIters(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvIters() failed with flag %i"%(ret))
def CVodeGetNumNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumNonlinSolvConvFails returns the number of nonlinear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVodeGetNumNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvConvFails() failed with flag %i"%(ret))
def CVodeGetNonlinSolvStats(cvodememobj): """A convenience function that provides the nonlinear solver optional outputs in a tuple (NumNonlinSolvIters, NumNonlinsolvConvFails). cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nniters = ctypes.c_long() nncfails = ctypes.c_long() ret = cvodes.CVodeGetNonlinSolvStats(cvodememobj.obj, ctypes.byref(nniters), ctypes.byref(nncfails)) if ret < 0:
def CVodeGetQuadNumRhsEvals(cvodememobj):
"""CVodeGetQuadNumRhsEvals returns the number of calls to the user function fQ defining the right hand side of the quadrature variables.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfQevals = ctypes.c_long(0)
ret = cvodes.CVodeGetQuadNumRhsEvals(cvodememobj.obj, ctypes.byref(nfQevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetQuadNumRhsEvals() failed with flag %i"%(ret))
def CVodeGetQuadNumErrTestFails(cvodememobj):
"""CVodeGetQuadNumErrTestFails returns the number of local error test failures for quadrature variables.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nQetfails = ctypes.c_long(0)
ret = cvodes.CVodeGetQuadNumErrTestFails(cvodememobj.obj, ctypes.byref(nQetfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetQuadNumErrTestFails() failed with flag %i"%(ret))
def CVodeGetQuadStats(cvodememobj): """A convenience function return a tuple of results from CVodeGetQuadNumRhsEvals and CVodeGetQuadErrTestFails respectively. cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nfQevals = ctypes.c_long(0) nQetfails = ctypes.c_long(0) ret = cvodes.CVodeGetQuadStats(cvodememobj.obj, ctypes.byref(nfQevals), ctypes.byref(nQetfails)) if ret < 0:
def CVodeGetNumSensRhsEvals(cvodememobj):
"""CVodeGetNumSensRhsEvals returns the number of calls to the sensitivity right hand side routine.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfSevals = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensRhsEvals(cvodememobj.obj, ctypes.byref(nfSevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensRhsEvals() failed with flag %i"%(ret))
def CVodeGetNumRhsEvalsSens(cvodememobj):
"""CVodeGetNumRhsEvalsSens returns the number of calls to the user f routine due to finite difference evaluations of the sensitivity equations.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfevalsS = ctypes.c_long(0)
ret = cvodes.CVodeGetNumRhsEvalsSens(cvodememobj.obj, ctypes.byref(nfevalsS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvalsSens() failed with flag %i"%(ret))
def CVodeGetNumSensErrTestFails(cvodememobj):
"""CVodeGetNumSensErrTestFails returns the number of local error test failures for sensitivity variables.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nSetfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensErrTestFails(cvodememobj.obj, ctypes.byref(nSetfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensErrTestFails() failed with flag %i"%(ret))
def CVodeGetNumSensLinSolvSetups(cvodememobj):
"""CVodeGetNumSensLinSolvSetups returns the number of calls made to the linear solver's setup routine due to sensitivity computations.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nlinsetupsS = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensLinSolvSetups(cvodememobj.obj, ctypes.byref(nlinsetupsS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensLinSolvSetups() failed with flag %i"%(ret))
def CVodeGetSensStats(cvodememobj): """A convenience function that provides the optional outputs in a tuple (NumSensRhsEvals, NumRhsEvalsSens, NumSensErrTestFails, NumSensLinSolvSetups) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nfSevals = ctypes.c_long(0) nfevalsS = ctypes.c_long(0) nSetfails = ctypes.c_long(0) nlinsetupS = ctypes.c_long(0)
def CVodeGetNumSensNonlinSolvIters(cvodememobj):
"""CVodeGetNumSensNonlinSolvIters returns the total number of nonlinear iterations for sensitivity variables.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nSniters = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensNonlinSolvIters(cvodememobj.obj, ctypes.byref(nSniters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensNonlinSolvIters() failed with flag %i"%(ret))
def CVodeGetNumSensNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumSensNonlinSolvConvFails returns the total number of nonlinear convergence failures for sensitivity variables
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nSncfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumSensNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(nSncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSensNonlinSolvConvFails() failed with flag %i"%(ret))
def CVodeGetNumStgrSensNonlinSolvIters(cvodememobj):
"""CVodeGetNumStgrSensNonlinSolvIters returns a vector of Ns nonlinear iteration counters for sensitivity variables in the CV_STAGGERED1 method.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nSTGR1niters = ctypes.c_long(0)
ret = cvodes.CVodeGetNumStgrSensNonlinSolvIters(cvodememobj.obj, ctypes.byref(nSTGR1niters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStgrSensNonlinSolvIters() failed with flag %i"%(ret))
def CVodeGetNumStgrSensNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumStgrSensNonlinSolvConvFails returns a vector of Ns nonlinear solver convergence failure counters for sensitivity variables in the CV_STAGGERED1 method.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nSTGR1ncfails = ctypes.c_long(0)
ret = cvodes.CVodeGetNumStgrSensNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(nSTGR1ncfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStgrSensNonlinSolvConvFails() failed with flag %i"%(ret))
def CVodeGetSensNonlinSolvStats(cvodememobj): """A convenience function that provides the optional outputs in a tuple (NumSensNonlinSolvIters, NumSensNonlinSolvConvFails, NumStgrSensNonlinSolvIters, NumStgrSensNonlinSolvConvFails) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nSniters = ctypes.c_long(0) nSncfails = ctypes.c_long(0) ret = cvodes.CVodeGetSensNonlinSolvStats(cvodememobj.obj, ctypes.byref(nSniters), ctypes.byref(nSncfails)) if ret < 0:
def CVDiagGetWorkSpace(cvodememobj): """Returns the CVDIAG real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvodes.CVDiagGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDiagGetNumRhsEvals(cvodememobj): """Returns the number of calls to the user f routine due to finite difference Jacobian evaluation. cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate() Note: The number of diagonal approximate Jacobians formed is equal to the number of CVDiagSetup calls. This number is available through CVodeGetNumLinSolvSetups.""" retval = ctypes.c_long(0) ret = cvodes.CVDiagGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
def CVBBDPrecGetWorkSpace(bbd_data): """CVBBDPrecGetWorkSpace returns the CVBBDPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvodes.CVBBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBBDPrecGetNumGfnEvals(bbd_data):
"""CVBBDPrecGetNumGfnEvals returns the number of calls to gfn."""
ngevalsBBDP = ctypes.c_long(0)
ret = cvodes.CVBBDPrecGetNumGfnEvals(bbd_data, ctypes.byref(ngevalsBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBBDPrecGetNumGfnEvals() failed with flag %i"%(ret))
return ngevalsBBDP.value
def CVBandPrecGetWorkSpace(bp_data): """CVBandPrecGetWorkSpace returns the CVBandPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvodes.CVBandPrecGetWorkSpace(bp_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandPrecGetNumRhsEvals(bp_data):
"""CVBandPrecGetNumGfnEvals returns the number of calls made from CVBANDPRE to the user's RHS."""
nfevalsBP = ctypes.c_long(0)
ret = cvodes.CVBandPrecGetNumRhsEvals(bp_data, ctypes.byref(nfevalsBP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandPrecGetNumRhsEvals() failed with flag %i"%(ret))
return nfevalsBP.value
def CVBandGetWorkSpace(cvodememobj): """CVBandGetWorkSpace returns the CVBand real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvodes.CVBandGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandGetNumJacEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's jacobian function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVBandGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumJacEvals() failed with flag %i"%(ret))
def CVBandGetNumRhsEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's RHS function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVBandGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumRhsEvals() failed with flag %i"%(ret))
def CVDenseGetWorkSpace(cvodememobj): """CVDenseGetWorkSpace returns the CVDense real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvodes.CVDenseGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDenseGetNumJacEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's Jacobian function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVDenseGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumJacEvals() failed with flag %i"%(ret))
def CVDenseGetNumRhsEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's RHS function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvodes.CVDenseGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumRhsEvals() failed with flag %i"%(ret))
def CVSpilsGetWorkSpace(cvodememobj): """CVSpilsGetWorkSpace returns the CVSpils real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvodes.CVSpilsGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVSpilsGetNumPrecEvals(cvodememobj):
"""CVSpilsGetNumPrecEvals returns the number of preconditioner evaluations, i.e. the number of calls made to PrecSetup with jok==FALSE.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npevals = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumPrecEvals(cvodememobj.obj, ctypes.byref(npevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecEvals() failed with flag %i"%(ret))
def CVSpilsGetNumPrecSolves(cvodememobj):
"""CVSpilsGetNumPrecSolves returns the number of calls made to PrecSolve.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npsolves = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumPrecSolves(cvodememobj.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecSolves() failed with flag %i"%(ret))
def CVSpilsGetNumLinIters(cvodememobj):
"""CVSpilsGetNumLinIters returns the number of linear iterations.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nliters = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumLinIters(cvodememobj.obj, ctypes.byref(nliters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumLinIters() failed with flag %i"%(ret))
def CVSpilsGetNumConvFails(cvodememobj):
"""CVSpilsGetNumConvFails returns the number of linear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nlcfails = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumConvFails(cvodememobj.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumConvFails() failed with flag %i"%(ret))
def CVSpilsGetNumJtimesEvals(cvodememobj):
"""CVSpilsGetNumJtimesEvals returns the number of calls to jtimes.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
njvevals = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumJtimesEvals(cvodememobj.obj, ctypes.byref(njvevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumJtimesEvals() failed with flag %i"%(ret))
def CVSpilsGetNumRhsEvals(cvodememobj):
"""CVSpilsGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian times vector evaluation.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfevalsLS = ctypes.c_long(0)
ret = cvodes.CVSpilsGetNumRhsEvals(cvodememobj.obj, ctypes.byref(nfevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumRhsEvals() failed with flag %i"%(ret))
src/p/y/pysces-HEAD/pysces_pysundials/pysundials/cvode.py pysces(Download)
def CVodeGetWorkSpace(cvodememobj): """CVodeGetWorkSpace returns the CVODE real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrw = ctypes.c_long() leniw = ctypes.c_long() ret = cvode.CVodeGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrw), ctypes.byref(leniw)) if ret < 0:
def CVodeGetNumSteps(cvodememobj):
"""CVodeGetNumSteps returns the cumulative number of internal steps taken by the solver
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumSteps(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSteps() failed with flag %i"%(ret))
def CVodeGetNumRhsEvals(cvodememobj):
"""CVodeGetNumRhsEvals returns the number of calls to the user's f function
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvals() failed with flag %i"%(ret))
def CVodeGetNumLinSolvSetups(cvodememobj):
"""CVodeGetNumLinSolvSetups returns the number of calls made to the linear solver's setup routine
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumLinSolvSetups(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumLinSolvSetups() failed with flag %i"%(ret))
def CVodeGetNumErrTestFails(cvodememobj):
"""CVodeGetNumErrTestFails returns the number of local error test failures that have occured
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumErrTestFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumErrTestFails() failed with flag %i"%(ret))
def CVodeGetNumStabLimOrderReds(cvodememobj):
"""CVodeGetNumStabLimOrderReds returns the number of order reductions due to stability limit detection
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumStabLimOrderReds(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStabLimOrderReds() failed with flag %i"%(ret))
def CVodeGetNumGEvals(cvodememobj):
"""CVodeGetNumGEvals returns the number of calls to the user's g function (for rootfinding)
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumGEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumGEvals() failed with flag %i"%(ret))
hlast last step size used by CVODE hcur next step size to be used by CVODE tcur current time""" nsteps = ctypes.c_long() nfevals = ctypes.c_long() nlinsetups = ctypes.c_long() netfails = ctypes.c_long()
def CVodeGetNumNonlinSolvIters(cvodememobj):
"""CVodeGetNumNonlinSolvIters returns the number of nonlinear solver iterations performed.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumNonlinSolvIters(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvIters() failed with flag %i"%(ret))
def CVodeGetNumNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumNonlinSolvConvFails returns the number of nonlinear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvConvFails() failed with flag %i"%(ret))
def CVodeGetNonlinSolvStats(cvodememobj): """A convenience function that provides the nonlinear solver optional outputs in a tuple (NumNonlinSolvIters, NumNonlinsolvConvFails). cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nniters = ctypes.c_long() nncfails = ctypes.c_long() ret = cvode.CVodeGetNonlinSolvStats(cvodememobj.obj, ctypes.byref(nniters), ctypes.byref(nncfails)) if ret < 0:
def CVDiagGetWorkSpace(cvodememobj): """Returns the CVDIAG real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVDiagGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDiagGetNumRhsEvals(cvodememobj): """Returns the number of calls to the user f routine due to finite difference Jacobian evaluation. cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate() Note: The number of diagonal approximate Jacobians formed is equal to the number of CVDiagSetup calls. This number is available through CVodeGetNumLinSolvSetups.""" retval = ctypes.c_long(0) ret = cvode.CVDiagGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
def CVBBDPrecGetWorkSpace(bbd_data): """CVBBDPrecGetWorkSpace returns the CVBBDPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVBBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBBDPrecGetNumGfnEvals(bbd_data):
"""CVBBDPrecGetNumGfnEvals returns the number of calls to gfn."""
ngevalsBBDP = ctypes.c_long(0)
ret = cvode.CVBBDPrecGetNumGfnEvals(bbd_data, ctypes.byref(ngevalsBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBBDPrecGetNumGfnEvals() failed with flag %i"%(ret))
return ngevalsBBDP.value
def CVBandPrecGetWorkSpace(bp_data): """CVBandPrecGetWorkSpace returns the CVBandPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVBandPrecGetWorkSpace(bp_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandPrecGetNumRhsEvals(bp_data):
"""CVBandPrecGetNumGfnEvals returns the number of calls made from CVBANDPRE to the user's RHS."""
nfevalsBP = ctypes.c_long(0)
ret = cvode.CVBandPrecGetNumRhsEvals(bp_data, ctypes.byref(nfevalsBP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandPrecGetNumRhsEvals() failed with flag %i"%(ret))
return nfevalsBP.value
def CVBandGetWorkSpace(cvodememobj): """CVBandGetWorkSpace returns the CVBand real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvode.CVBandGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandGetNumJacEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's jacobian function."""
retval = ctypes.c_long(0)
ret = cvode.CVBandGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumJacEvals() failed with flag %i"%(ret))
return retval.value
cvode.CVBandGetNumJacEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvode.CVBandGetNumJacEvals.restype = ctypes.c_int
def CVBandGetNumRhsEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's RHS function."""
retval = ctypes.c_long(0)
def CVDenseGetWorkSpace(cvodememobj): """CVDenseGetWorkSpace returns the CVDense real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvode.CVDenseGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDenseGetNumJacEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's Jacobian function."""
retval = ctypes.c_long(0)
ret = cvode.CVDenseGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumJacEvals() failed with flag %i"%(ret))
return retval.value
cvode.CVDenseGetNumJacEvals.argtypes = [ctypes.c_void_p, ctypes.POINTER(ctypes.c_long)]
cvode.CVDenseGetNumJacEvals.restype = ctypes.c_int
def CVDenseGetNumRhsEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's RHS function."""
retval = ctypes.c_long(0)
def CVSpilsGetWorkSpace(cvodememobj): """CVSpilsGetWorkSpace returns the CVSpils real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVSpilsGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVSpilsGetNumPrecEvals(cvodememobj):
"""CVSpilsGetNumPrecEvals returns the number of preconditioner evaluations, i.e. the number of calls made to PrecSetup with jok==FALSE.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npevals = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumPrecEvals(cvodememobj.obj, ctypes.byref(npevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecEvals() failed with flag %i"%(ret))
def CVSpilsGetNumPrecSolves(cvodememobj):
"""CVSpilsGetNumPrecSolves returns the number of calls made to PrecSolve.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npsolves = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumPrecSolves(cvodememobj.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecSolves() failed with flag %i"%(ret))
def CVSpilsGetNumLinIters(cvodememobj):
"""CVSpilsGetNumLinIters returns the number of linear iterations.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nliters = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumLinIters(cvodememobj.obj, ctypes.byref(nliters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumLinIters() failed with flag %i"%(ret))
def CVSpilsGetNumConvFails(cvodememobj):
"""CVSpilsGetNumConvFails returns the number of linear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nlcfails = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumConvFails(cvodememobj.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumConvFails() failed with flag %i"%(ret))
def CVSpilsGetNumJtimesEvals(cvodememobj):
"""CVSpilsGetNumJtimesEvals returns the number of calls to jtimes.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
njvevals = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumJtimesEvals(cvodememobj.obj, ctypes.byref(njvevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumJtimesEvals() failed with flag %i"%(ret))
def CVSpilsGetNumRhsEvals(cvodememobj):
"""CVSpilsGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian times vector evaluation.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfevalsLS = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumRhsEvals(cvodememobj.obj, ctypes.byref(nfevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumRhsEvals() failed with flag %i"%(ret))
src/p/y/pysundials-HEAD/2.3.0/src/cvode.py pysundials(Download)
def CVodeGetWorkSpace(cvodememobj): """CVodeGetWorkSpace returns the CVODE real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrw = ctypes.c_long() leniw = ctypes.c_long() ret = cvode.CVodeGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrw), ctypes.byref(leniw)) if ret < 0:
def CVodeGetNumSteps(cvodememobj):
"""CVodeGetNumSteps returns the cumulative number of internal steps taken by the solver
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumSteps(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumSteps() failed with flag %i"%(ret))
def CVodeGetNumRhsEvals(cvodememobj):
"""CVodeGetNumRhsEvals returns the number of calls to the user's f function
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumRhsEvals() failed with flag %i"%(ret))
def CVodeGetNumLinSolvSetups(cvodememobj):
"""CVodeGetNumLinSolvSetups returns the number of calls made to the linear solver's setup routine
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumLinSolvSetups(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumLinSolvSetups() failed with flag %i"%(ret))
def CVodeGetNumErrTestFails(cvodememobj):
"""CVodeGetNumErrTestFails returns the number of local error test failures that have occured
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumErrTestFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumErrTestFails() failed with flag %i"%(ret))
def CVodeGetNumStabLimOrderReds(cvodememobj):
"""CVodeGetNumStabLimOrderReds returns the number of order reductions due to stability limit detection
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long()
ret = cvode.CVodeGetNumStabLimOrderReds(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumStabLimOrderReds() failed with flag %i"%(ret))
def CVodeGetNumGEvals(cvodememobj):
"""CVodeGetNumGEvals returns the number of calls to the user's g function (for rootfinding)
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumGEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumGEvals() failed with flag %i"%(ret))
hlast last step size used by CVODE hcur next step size to be used by CVODE tcur current time""" nsteps = ctypes.c_long() nfevals = ctypes.c_long() nlinsetups = ctypes.c_long() netfails = ctypes.c_long()
def CVodeGetNumNonlinSolvIters(cvodememobj):
"""CVodeGetNumNonlinSolvIters returns the number of nonlinear solver iterations performed.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumNonlinSolvIters(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvIters() failed with flag %i"%(ret))
def CVodeGetNumNonlinSolvConvFails(cvodememobj):
"""CVodeGetNumNonlinSolvConvFails returns the number of nonlinear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVodeGetNumNonlinSolvConvFails(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVodeGetNumNonlinSolvConvFails() failed with flag %i"%(ret))
def CVodeGetNonlinSolvStats(cvodememobj): """A convenience function that provides the nonlinear solver optional outputs in a tuple (NumNonlinSolvIters, NumNonlinsolvConvFails). cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()""" nniters = ctypes.c_long() nncfails = ctypes.c_long() ret = cvode.CVodeGetNonlinSolvStats(cvodememobj.obj, ctypes.byref(nniters), ctypes.byref(nncfails)) if ret < 0:
def CVDiagGetWorkSpace(cvodememobj): """Returns the CVDIAG real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVDiagGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDiagGetNumRhsEvals(cvodememobj): """Returns the number of calls to the user f routine due to finite difference Jacobian evaluation. cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate() Note: The number of diagonal approximate Jacobians formed is equal to the number of CVDiagSetup calls. This number is available through CVodeGetNumLinSolvSetups.""" retval = ctypes.c_long(0) ret = cvode.CVDiagGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
def CVBBDPrecGetWorkSpace(bbd_data): """CVBBDPrecGetWorkSpace returns the CVBBDPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVBBDPrecGetWorkSpace(bbd_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBBDPrecGetNumGfnEvals(bbd_data):
"""CVBBDPrecGetNumGfnEvals returns the number of calls to gfn."""
ngevalsBBDP = ctypes.c_long(0)
ret = cvode.CVBBDPrecGetNumGfnEvals(bbd_data, ctypes.byref(ngevalsBBDP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBBDPrecGetNumGfnEvals() failed with flag %i"%(ret))
return ngevalsBBDP.value
def CVBandPrecGetWorkSpace(bp_data): """CVBandPrecGetWorkSpace returns the CVBandPrec real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVBandPrecGetWorkSpace(bp_data, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandPrecGetNumRhsEvals(bp_data):
"""CVBandPrecGetNumGfnEvals returns the number of calls made from CVBANDPRE to the user's RHS."""
nfevalsBP = ctypes.c_long(0)
ret = cvode.CVBandPrecGetNumRhsEvals(bp_data, ctypes.byref(nfevalsBP))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandPrecGetNumRhsEvals() failed with flag %i"%(ret))
return nfevalsBP.value
def CVBandGetWorkSpace(cvodememobj): """CVBandGetWorkSpace returns the CVBand real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvode.CVBandGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVBandGetNumJacEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's jacobian function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVBandGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumJacEvals() failed with flag %i"%(ret))
def CVBandGetNumRhsEvals(cvodememobj):
"""CVBandGetNumGfnEvals returns the number of calls made from CVBAND to the user's RHS function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVBandGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVBandGetNumRhsEvals() failed with flag %i"%(ret))
def CVDenseGetWorkSpace(cvodememobj): """CVDenseGetWorkSpace returns the CVDense real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long() leniwLS = ctypes.c_long() ret = cvode.CVDenseGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVDenseGetNumJacEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's Jacobian function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVDenseGetNumJacEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumJacEvals() failed with flag %i"%(ret))
def CVDenseGetNumRhsEvals(cvodememobj):
"""CVDenseGetNumGfnEvals returns the number of calls made from CVDENSE to the user's RHS function.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()"""
retval = ctypes.c_long(0)
ret = cvode.CVDenseGetNumRhsEvals(cvodememobj.obj, ctypes.byref(retval))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVDenseGetNumRhsEvals() failed with flag %i"%(ret))
def CVSpilsGetWorkSpace(cvodememobj): """CVSpilsGetWorkSpace returns the CVSpils real and integer workspaces as a tuple (integer, real) cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CVodeCreate()""" lenrwLS = ctypes.c_long(0) leniwLS = ctypes.c_long(0) ret = cvode.CVSpilsGetWorkSpace(cvodememobj.obj, ctypes.byref(lenrwLS), ctypes.byref(leniwLS)) if ret < 0:
def CVSpilsGetNumPrecEvals(cvodememobj):
"""CVSpilsGetNumPrecEvals returns the number of preconditioner evaluations, i.e. the number of calls made to PrecSetup with jok==FALSE.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npevals = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumPrecEvals(cvodememobj.obj, ctypes.byref(npevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecEvals() failed with flag %i"%(ret))
def CVSpilsGetNumPrecSolves(cvodememobj):
"""CVSpilsGetNumPrecSolves returns the number of calls made to PrecSolve.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
npsolves = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumPrecSolves(cvodememobj.obj, ctypes.byref(npsolves))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumPrecSolves() failed with flag %i"%(ret))
def CVSpilsGetNumLinIters(cvodememobj):
"""CVSpilsGetNumLinIters returns the number of linear iterations.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nliters = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumLinIters(cvodememobj.obj, ctypes.byref(nliters))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumLinIters() failed with flag %i"%(ret))
def CVSpilsGetNumConvFails(cvodememobj):
"""CVSpilsGetNumConvFails returns the number of linear convergence failures.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nlcfails = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumConvFails(cvodememobj.obj, ctypes.byref(nlcfails))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumConvFails() failed with flag %i"%(ret))
def CVSpilsGetNumJtimesEvals(cvodememobj):
"""CVSpilsGetNumJtimesEvals returns the number of calls to jtimes.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
njvevals = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumJtimesEvals(cvodememobj.obj, ctypes.byref(njvevals))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumJtimesEvals() failed with flag %i"%(ret))
def CVSpilsGetNumRhsEvals(cvodememobj):
"""CVSpilsGetNumRhsEvals returns the number of calls to the user f routine due to finite difference Jacobian times vector evaluation.
cvodememobj (CVodeMemObj) a CVodeMemObj as returned by CvodeCreate()"""
nfevalsLS = ctypes.c_long(0)
ret = cvode.CVSpilsGetNumRhsEvals(cvodememobj.obj, ctypes.byref(nfevalsLS))
if ret < 0:
raise AssertionError("SUNDIALS ERROR: CVSpilsGetNumRhsEvals() failed with flag %i"%(ret))
src/l/a/LabJackPython-HEAD/src/u12.py LabJackPython(Download)
idNum = self.id
if ON_WINDOWS:
ljid = ctypes.c_long(idNum)
ad0 = ctypes.c_long(999)
ad1 = ctypes.c_float(999)
idNum = self.id
if ON_WINDOWS:
ljid = ctypes.c_long(idNum)
ecode = staticLib.EAnalogOut(ctypes.byref(ljid), demo, ctypes.c_float(analogOut0), ctypes.c_float(analogOut1))
if ecode != 0: raise U12Exception(ecode)
idNum = self.id
if ON_WINDOWS:
ljid = ctypes.c_long(idNum)
count = ctypes.c_double()
ms = ctypes.c_double()
idNum = self.id
if ON_WINDOWS:
ljid = ctypes.c_long(idNum)
state = ctypes.c_long(999)
ecode = staticLib.EDigitalIn(ctypes.byref(ljid), demo, channel, readD, ctypes.byref(state))
idNum = self.id
if ON_WINDOWS:
ljid = ctypes.c_long(idNum)
ecode = staticLib.EDigitalOut(ctypes.byref(ljid), demo, channel, writeD, state)
# Check id num
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Check to make sure that everything is checked
if not isIterable(channels): raise TypeError("channels must be iterable")
# Convert lists to arrays and create other ctypes
channelsArray = listToCArray(channels, ctypes.c_long)
gainsArray = listToCArray(gains, ctypes.c_long)
overVoltage = ctypes.c_long(999)
longArrayType = (ctypes.c_long * 4)
floatArrayType = (ctypes.c_float * 4)
voltages = floatArrayType(0, 0, 0, 0)
stateIOin = ctypes.c_long(stateIOin)
# Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# check list sizes
if len(channels) < numChannels: raise ValueError("channels must have atleast numChannels elements")
voltages_type = arr4096_type * 4
voltages = voltages_type()
stateIOout = (ctypes.c_long * 4096)()
overVoltage = ctypes.c_long(999)
ecode = staticLib.AIBurst(ctypes.byref(idNum), demo, stateIOin, updateIO, ledOn, numChannels, ctypes.byref(channelsArray), ctypes.byref(gainsArray), ctypes.byref(scanRate), disableCal, triggerIO, triggerState, numScans, timeout, ctypes.byref(voltages), ctypes.byref(stateIOout), ctypes.byref(overVoltage), transferMode)
# Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Convert lists to arrays and create other ctypes
channelsArray = listToCArray(channels, ctypes.c_long)
voltages_type = arr4096_type * 4
voltages = voltages_type()
stateIOout = (ctypes.c_long * 4096)()
reserved = ctypes.c_long(0)
ljScanBacklog = ctypes.c_long(99999)
overVoltage = ctypes.c_long(999)
# Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Check tris and state arguments
if updateDigital > 0:
if trisD is None: raise ValueError("keyword argument trisD must be set")
if trisIO is None: raise ValueError("keyword argument trisIO must be set")
if stateD is None: raise ValueError("keyword argument stateD must be set")
if stateIO is None: raise ValueError("keyword argument stateIO must be set")
# Create ctypes
if stateD is None: stateD = ctypes.c_long(0)
# Create ctypes
if stateD is None: stateD = ctypes.c_long(0)
else: stateD = ctypes.c_long(stateD)
if stateIO is None: stateIO = ctypes.c_long(0)
else: stateIO = ctypes.c_long(stateIO)
count = ctypes.c_ushort(999)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.AsynchConfig(ctypes.byref(idNum), demo, timeoutMult, configA, configB, configTE, fullA, fullB, fullC, halfA, halfB, halfC)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Check size of data
if len(data) > 18: raise ValueError("data can not be larger than 18 elements")
>>> {'bits': 2662}
"""
if ON_WINDOWS:
bits = ctypes.c_long(999)
ecode = staticLib.VoltsToBits(chnum, chgain, ctypes.c_float(volts), ctypes.byref(bits))
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Create ctypes
stateD = ctypes.c_long(999)
stateIO = ctypes.c_long(999)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Check tris and state parameters
if updateDigital > 0:
if trisD is None: raise ValueError("keyword argument trisD must be set")
if trisIO is None: raise ValueError("keyword argument trisIO must be set")
if stateD is None: raise ValueError("keyword argument stateD must be set")
if stateIO is None: raise ValueError("keyword argument stateIO must be set")
# Create ctypes
if trisD is None: trisD = ctypes.c_long(999)
# Create ctypes
if trisD is None: trisD = ctypes.c_long(999)
else:trisD = ctypes.c_long(trisD)
if stateD is None:stateD = ctypes.c_long(999)
else: stateD = ctypes.c_long(stateD)
if stateIO is None: stateIO = ctypes.c_long(0)
else: stateIO = ctypes.c_long(stateIO)
outputD = ctypes.c_long(999)
# Check ID number
if idNum is None: idNum = self.id
idNum = ctypes.c_long(idNum)
staticLib.GetFirmwareVersion.restype = ctypes.c_float
firmware = staticLib.GetFirmwareVersion(ctypes.byref(idNum))
arr127_type = ctypes.c_long * 127
calMatrix_type = arr127_type * 20
calMatrix = calMatrix_type()
reserved = ctypes.c_long()
numberFound = ctypes.c_long()
ecode = staticLib.ListAll(ctypes.byref(productIDList), ctypes.byref(serialnumList), ctypes.byref(localIDList), ctypes.byref(powerList), ctypes.byref(calMatrix), ctypes.byref(numberFound), ctypes.byref(reserved), ctypes.byref(reserved))
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.LocalID(ctypes.byref(idNum), localID)
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.NoThread(ctypes.byref(idNum), noThread)
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.PulseOut(ctypes.byref(idNum), demo, lowFirst, bitSelect, numPulses, timeB1, timeC1, timeB2, timeC2)
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.PulseOutStart(ctypes.byref(idNum), demo, lowFirst, bitSelect, numPulses, timeB1, timeC1, timeB2, timeC2)
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.PulseOutFinish(ctypes.byref(idNum), demo, timeoutMS)
if ecode != 0: raise U12Exception(ecode)
# Create ctypes
frequency = ctypes.c_float(frequency)
timeB = ctypes.c_long(0)
timeC = ctypes.c_long(0)
ecode = staticLib.PulseOutCalc(ctypes.byref(frequency), ctypes.byref(timeB), ctypes.byref(timeC))
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.ReEnum(ctypes.byref(idNum))
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
ecode = staticLib.Reset(ctypes.byref(idNum))
if ecode != 0: raise U12Exception(ecode)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Create ctypes
tempC = ctypes.c_float(0)
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
# Check size of datatx
if len(datatx) != 4: raise ValueError("datatx must have exactly 4 elements")
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
if controlCS > 0 and csLine is None: raise ValueError("csLine must be specified")
#Check id number
if idNum is None:
idNum = self.id
idNum = ctypes.c_long(idNum)
if len(activeDn) is not 3: raise ValueError("activeDn must have 3 elements")
if len(stateDn) is not 3: raise Value("stateDn must have 3 elements")
ad2 = ctypes.c_ulong()
ad3 = ctypes.c_ulong()
ec = staticLib.ReadMem(ctypes.byref(ljid), ctypes.c_long(address), ctypes.byref(ad3), ctypes.byref(ad2), ctypes.byref(ad1), ctypes.byref(ad0))
if ec != 0: raise U12Exception(ec)
addr = [0] * 4
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 Next