Generated by Cython 0.15.1 on Wed Jul 18 12:11:31 2012
Raw output: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.c
1:
2: include "interrupt.pxi" # ctrl-c interrupt block support
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":2 * * include "interrupt.pxi" # ctrl-c interrupt block support #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* include "stdsage.pxi" # ctrl-c interrupt block support * */ __pyx_t_2 = PyDict_New(); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 2; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_2)); if (PyObject_SetAttr(__pyx_m, __pyx_n_s____test__, ((PyObject *)__pyx_t_2))<
0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 2; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(((PyObject *)__pyx_t_2)); __pyx_t_2 = 0;
3: include "stdsage.pxi" # ctrl-c interrupt block support
4:
5: include "cdefs.pxi"
6: import numpy as np
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":6 * * include "cdefs.pxi" * import numpy as np #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cimport numpy as np * */ __pyx_t_1 = __Pyx_Import(((PyObject *)__pyx_n_s__numpy), 0, -1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 6; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (PyObject_SetAttr(__pyx_m, __pyx_n_s__np, __pyx_t_1)<
0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 6; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
7: cimport numpy as np
8:
9: cdouble = np.double
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":9 * cimport numpy as np * * cdouble = np.double #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* ctypedef np.double_t cdouble_t * */ __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_2 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__double); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; if (PyObject_SetAttr(__pyx_m, __pyx_n_s__cdouble, __pyx_t_2)<
0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 9; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
10: ctypedef np.double_t cdouble_t
11:
12: cpdef np.ndarray ode_DP(np.ndarray[cdouble_t] x, double t,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":12 * ctypedef np.double_t cdouble_t * * cpdef np.ndarray ode_DP(np.ndarray[cdouble_t] x, double t, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* double a_00, double a_01, double a_10, * double b_00, double b_01, double b_10, */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_ode_DP(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static PyArrayObject *__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_ode_DP(PyArrayObject *__pyx_v_x, double __pyx_v_t, double __pyx_v_a_00, double __pyx_v_a_01, double __pyx_v_a_10, double __pyx_v_b_00, double __pyx_v_b_01, double __pyx_v_b_10, double __pyx_v_c_0001, double __pyx_v_c_0010, double __pyx_v_c_0011, double __pyx_v_c_0111, double __pyx_v_c_1011, double __pyx_v_alpha_01, double __pyx_v_alpha_10, double __pyx_v_alpha_11, double __pyx_v_beta_01, double __pyx_v_beta_10, double __pyx_v_beta_11, double __pyx_v_gamma_0100, double __pyx_v_gamma_1000, double __pyx_v_gamma_1100, double __pyx_v_gamma_1101, double __pyx_v_gamma_1110, int __pyx_skip_dispatch) { double __pyx_v_ks_00; double __pyx_v_ks_01; double __pyx_v_ks_10; double __pyx_v_fs_01; double __pyx_v_fs_10; double __pyx_v_fs_11; double __pyx_v_k; double __pyx_v_f; double __pyx_v_s_00; double __pyx_v_s_01; double __pyx_v_s_10; double __pyx_v_s_11; Py_buffer __pyx_bstruct_x; Py_ssize_t __pyx_bstride_0_x = 0; Py_ssize_t __pyx_bshape_0_x = 0; PyArrayObject *__pyx_r = NULL; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("ode_DP"); __pyx_bstruct_x.buf = NULL; { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_x, (PyObject*)__pyx_v_x, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_x = __pyx_bstruct_x.strides[0]; __pyx_bshape_0_x = __pyx_bstruct_x.shape[0]; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":12 * ctypedef np.double_t cdouble_t * * cpdef np.ndarray ode_DP(np.ndarray[cdouble_t] x, double t, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* double a_00, double a_01, double a_10, * double b_00, double b_01, double b_10, */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_ode_DP(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static char __pyx_doc_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_ode_DP[] = "File: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx (starting at line 12)"; static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_ode_DP(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) { PyArrayObject *__pyx_v_x = 0; double __pyx_v_t; double __pyx_v_a_00; double __pyx_v_a_01; double __pyx_v_a_10; double __pyx_v_b_00; double __pyx_v_b_01; double __pyx_v_b_10; double __pyx_v_c_0001; double __pyx_v_c_0010; double __pyx_v_c_0011; double __pyx_v_c_0111; double __pyx_v_c_1011; double __pyx_v_alpha_01; double __pyx_v_alpha_10; double __pyx_v_alpha_11; double __pyx_v_beta_01; double __pyx_v_beta_10; double __pyx_v_beta_11; double __pyx_v_gamma_0100; double __pyx_v_gamma_1000; double __pyx_v_gamma_1100; double __pyx_v_gamma_1101; double __pyx_v_gamma_1110; Py_buffer __pyx_bstruct_x; Py_ssize_t __pyx_bstride_0_x = 0; Py_ssize_t __pyx_bshape_0_x = 0; PyObject *__pyx_r = NULL; static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__x,&__pyx_n_s__t,&__pyx_n_s__a_00,&__pyx_n_s__a_01,&__pyx_n_s__a_10,&__pyx_n_s__b_00,&__pyx_n_s__b_01,&__pyx_n_s__b_10,&__pyx_n_s__c_0001,&__pyx_n_s__c_0010,&__pyx_n_s__c_0011,&__pyx_n_s__c_0111,&__pyx_n_s__c_1011,&__pyx_n_s__alpha_01,&__pyx_n_s__alpha_10,&__pyx_n_s__alpha_11,&__pyx_n_s__beta_01,&__pyx_n_s__beta_10,&__pyx_n_s__beta_11,&__pyx_n_s__gamma_0100,&__pyx_n_s__gamma_1000,&__pyx_n_s__gamma_1100,&__pyx_n_s__gamma_1101,&__pyx_n_s__gamma_1110,0}; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("ode_DP"); __pyx_self = __pyx_self; { PyObject* values[24] = {0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0}; if (unlikely(__pyx_kwds)) { Py_ssize_t kw_args; switch (PyTuple_GET_SIZE(__pyx_args)) { case 24: values[23] = PyTuple_GET_ITEM(__pyx_args, 23); case 23: values[22] = PyTuple_GET_ITEM(__pyx_args, 22); case 22: values[21] = PyTuple_GET_ITEM(__pyx_args, 21); case 21: values[20] = PyTuple_GET_ITEM(__pyx_args, 20); case 20: values[19] = PyTuple_GET_ITEM(__pyx_args, 19); case 19: values[18] = PyTuple_GET_ITEM(__pyx_args, 18); case 18: values[17] = PyTuple_GET_ITEM(__pyx_args, 17); case 17: values[16] = PyTuple_GET_ITEM(__pyx_args, 16); case 16: values[15] = PyTuple_GET_ITEM(__pyx_args, 15); case 15: values[14] = PyTuple_GET_ITEM(__pyx_args, 14); case 14: values[13] = PyTuple_GET_ITEM(__pyx_args, 13); case 13: values[12] = PyTuple_GET_ITEM(__pyx_args, 12); case 12: values[11] = PyTuple_GET_ITEM(__pyx_args, 11); case 11: values[10] = PyTuple_GET_ITEM(__pyx_args, 10); case 10: values[9] = PyTuple_GET_ITEM(__pyx_args, 9); case 9: values[8] = PyTuple_GET_ITEM(__pyx_args, 8); case 8: values[7] = PyTuple_GET_ITEM(__pyx_args, 7); case 7: values[6] = PyTuple_GET_ITEM(__pyx_args, 6); case 6: values[5] = PyTuple_GET_ITEM(__pyx_args, 5); case 5: values[4] = PyTuple_GET_ITEM(__pyx_args, 4); case 4: values[3] = PyTuple_GET_ITEM(__pyx_args, 3); case 3: values[2] = PyTuple_GET_ITEM(__pyx_args, 2); case 2: values[1] = PyTuple_GET_ITEM(__pyx_args, 1); case 1: values[0] = PyTuple_GET_ITEM(__pyx_args, 0); case 0: break; default: goto __pyx_L5_argtuple_error; } kw_args = PyDict_Size(__pyx_kwds); switch (PyTuple_GET_SIZE(__pyx_args)) { case 0: values[0] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__x); if (likely(values[0])) kw_args--; else goto __pyx_L5_argtuple_error; case 1: values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__t); if (likely(values[1])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 2: values[2] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__a_00); if (likely(values[2])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 3: values[3] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__a_01); if (likely(values[3])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 3); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 4: values[4] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__a_10); if (likely(values[4])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 4); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 5: values[5] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__b_00); if (likely(values[5])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 5); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 6: values[6] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__b_01); if (likely(values[6])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 6); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 7: values[7] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__b_10); if (likely(values[7])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 7); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 8: values[8] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__c_0001); if (likely(values[8])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 8); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 9: values[9] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__c_0010); if (likely(values[9])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 9); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 10: values[10] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__c_0011); if (likely(values[10])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 10); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 11: values[11] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__c_0111); if (likely(values[11])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 11); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 12: values[12] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__c_1011); if (likely(values[12])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 12); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 13: values[13] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__alpha_01); if (likely(values[13])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 13); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 14: values[14] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__alpha_10); if (likely(values[14])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 14); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 15: values[15] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__alpha_11); if (likely(values[15])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 15); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 16: values[16] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__beta_01); if (likely(values[16])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 16); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 17: values[17] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__beta_10); if (likely(values[17])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 17); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 18: values[18] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__beta_11); if (likely(values[18])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 18); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 19: values[19] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__gamma_0100); if (likely(values[19])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 19); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 20: values[20] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__gamma_1000); if (likely(values[20])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 20); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 21: values[21] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__gamma_1100); if (likely(values[21])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 21); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 22: values[22] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__gamma_1101); if (likely(values[22])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 22); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 23: values[23] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__gamma_1110); if (likely(values[23])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, 23); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } } if (unlikely(kw_args > 0)) { if (unlikely(__Pyx_ParseOptionalKeywords(__pyx_kwds, __pyx_pyargnames, 0, values, PyTuple_GET_SIZE(__pyx_args), "ode_DP")<
0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } } else if (PyTuple_GET_SIZE(__pyx_args) != 24) { goto __pyx_L5_argtuple_error; } else { values[0] = PyTuple_GET_ITEM(__pyx_args, 0); values[1] = PyTuple_GET_ITEM(__pyx_args, 1); values[2] = PyTuple_GET_ITEM(__pyx_args, 2); values[3] = PyTuple_GET_ITEM(__pyx_args, 3); values[4] = PyTuple_GET_ITEM(__pyx_args, 4); values[5] = PyTuple_GET_ITEM(__pyx_args, 5); values[6] = PyTuple_GET_ITEM(__pyx_args, 6); values[7] = PyTuple_GET_ITEM(__pyx_args, 7); values[8] = PyTuple_GET_ITEM(__pyx_args, 8); values[9] = PyTuple_GET_ITEM(__pyx_args, 9); values[10] = PyTuple_GET_ITEM(__pyx_args, 10); values[11] = PyTuple_GET_ITEM(__pyx_args, 11); values[12] = PyTuple_GET_ITEM(__pyx_args, 12); values[13] = PyTuple_GET_ITEM(__pyx_args, 13); values[14] = PyTuple_GET_ITEM(__pyx_args, 14); values[15] = PyTuple_GET_ITEM(__pyx_args, 15); values[16] = PyTuple_GET_ITEM(__pyx_args, 16); values[17] = PyTuple_GET_ITEM(__pyx_args, 17); values[18] = PyTuple_GET_ITEM(__pyx_args, 18); values[19] = PyTuple_GET_ITEM(__pyx_args, 19); values[20] = PyTuple_GET_ITEM(__pyx_args, 20); values[21] = PyTuple_GET_ITEM(__pyx_args, 21); values[22] = PyTuple_GET_ITEM(__pyx_args, 22); values[23] = PyTuple_GET_ITEM(__pyx_args, 23); } __pyx_v_x = ((PyArrayObject *)values[0]); __pyx_v_t = __pyx_PyFloat_AsDouble(values[1]); if (unlikely((__pyx_v_t == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_a_00 = __pyx_PyFloat_AsDouble(values[2]); if (unlikely((__pyx_v_a_00 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_a_01 = __pyx_PyFloat_AsDouble(values[3]); if (unlikely((__pyx_v_a_01 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_a_10 = __pyx_PyFloat_AsDouble(values[4]); if (unlikely((__pyx_v_a_10 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_b_00 = __pyx_PyFloat_AsDouble(values[5]); if (unlikely((__pyx_v_b_00 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_b_01 = __pyx_PyFloat_AsDouble(values[6]); if (unlikely((__pyx_v_b_01 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_b_10 = __pyx_PyFloat_AsDouble(values[7]); if (unlikely((__pyx_v_b_10 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_c_0001 = __pyx_PyFloat_AsDouble(values[8]); if (unlikely((__pyx_v_c_0001 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 15; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_c_0010 = __pyx_PyFloat_AsDouble(values[9]); if (unlikely((__pyx_v_c_0010 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 15; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_c_0011 = __pyx_PyFloat_AsDouble(values[10]); if (unlikely((__pyx_v_c_0011 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 15; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_c_0111 = __pyx_PyFloat_AsDouble(values[11]); if (unlikely((__pyx_v_c_0111 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 16; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_c_1011 = __pyx_PyFloat_AsDouble(values[12]); if (unlikely((__pyx_v_c_1011 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 16; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_alpha_01 = __pyx_PyFloat_AsDouble(values[13]); if (unlikely((__pyx_v_alpha_01 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_alpha_10 = __pyx_PyFloat_AsDouble(values[14]); if (unlikely((__pyx_v_alpha_10 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_alpha_11 = __pyx_PyFloat_AsDouble(values[15]); if (unlikely((__pyx_v_alpha_11 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_beta_01 = __pyx_PyFloat_AsDouble(values[16]); if (unlikely((__pyx_v_beta_01 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_beta_10 = __pyx_PyFloat_AsDouble(values[17]); if (unlikely((__pyx_v_beta_10 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_beta_11 = __pyx_PyFloat_AsDouble(values[18]); if (unlikely((__pyx_v_beta_11 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_gamma_0100 = __pyx_PyFloat_AsDouble(values[19]); if (unlikely((__pyx_v_gamma_0100 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 19; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_gamma_1000 = __pyx_PyFloat_AsDouble(values[20]); if (unlikely((__pyx_v_gamma_1000 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 19; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_gamma_1100 = __pyx_PyFloat_AsDouble(values[21]); if (unlikely((__pyx_v_gamma_1100 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 19; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_gamma_1101 = __pyx_PyFloat_AsDouble(values[22]); if (unlikely((__pyx_v_gamma_1101 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_v_gamma_1110 = __pyx_PyFloat_AsDouble(values[23]); if (unlikely((__pyx_v_gamma_1110 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } goto __pyx_L4_argument_unpacking_done; __pyx_L5_argtuple_error:; __Pyx_RaiseArgtupleInvalid("ode_DP", 1, 24, 24, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_L3_error:; __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.ode_DP", __pyx_clineno, __pyx_lineno, __pyx_filename); __Pyx_RefNannyFinishContext(); return NULL; __pyx_L4_argument_unpacking_done:; __pyx_bstruct_x.buf = NULL; if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x), __pyx_ptype_5numpy_ndarray, 1, "x", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;} { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_x, (PyObject*)__pyx_v_x, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_x = __pyx_bstruct_x.strides[0]; __pyx_bshape_0_x = __pyx_bstruct_x.shape[0]; __Pyx_XDECREF(__pyx_r); __pyx_t_1 = ((PyObject *)__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_ode_DP(__pyx_v_x, __pyx_v_t, __pyx_v_a_00, __pyx_v_a_01, __pyx_v_a_10, __pyx_v_b_00, __pyx_v_b_01, __pyx_v_b_10, __pyx_v_c_0001, __pyx_v_c_0010, __pyx_v_c_0011, __pyx_v_c_0111, __pyx_v_c_1011, __pyx_v_alpha_01, __pyx_v_alpha_10, __pyx_v_alpha_11, __pyx_v_beta_01, __pyx_v_beta_10, __pyx_v_beta_11, __pyx_v_gamma_0100, __pyx_v_gamma_1000, __pyx_v_gamma_1100, __pyx_v_gamma_1101, __pyx_v_gamma_1110, 0)); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 12; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_r = __pyx_t_1; __pyx_t_1 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_1); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);} __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.ode_DP", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = NULL; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
13: double a_00, double a_01, double a_10,
14: double b_00, double b_01, double b_10,
15: double c_0001, double c_0010, double c_0011,
16: double c_0111, double c_1011,
17: double alpha_01, double alpha_10, double alpha_11,
18: double beta_01, double beta_10, double beta_11,
19: double gamma_0100, double gamma_1000, double gamma_1100,
20: double gamma_1101, double gamma_1110):
21: cdef double ks_00 = x[ 0]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":21 * double gamma_0100, double gamma_1000, double gamma_1100, * double gamma_1101, double gamma_1110): * cdef double ks_00 = x[ 0] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double ks_01 = x[ 1] * cdef double ks_10 = x[ 2] */ __pyx_t_1 = 0; __pyx_t_2 = -1; if (__pyx_t_1<
0) { __pyx_t_1 += __pyx_bshape_0_x; if (unlikely(__pyx_t_1<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_1 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 21; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_ks_00 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_1, __pyx_bstride_0_x));
22: cdef double ks_01 = x[ 1]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":22 * double gamma_1101, double gamma_1110): * cdef double ks_00 = x[ 0] * cdef double ks_01 = x[ 1] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double ks_10 = x[ 2] * cdef double fs_01 = x[ 3] */ __pyx_t_3 = 1; __pyx_t_2 = -1; if (__pyx_t_3<
0) { __pyx_t_3 += __pyx_bshape_0_x; if (unlikely(__pyx_t_3<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_3 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 22; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_ks_01 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_3, __pyx_bstride_0_x));
23: cdef double ks_10 = x[ 2]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":23 * cdef double ks_00 = x[ 0] * cdef double ks_01 = x[ 1] * cdef double ks_10 = x[ 2] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double fs_01 = x[ 3] * cdef double fs_10 = x[ 4] */ __pyx_t_4 = 2; __pyx_t_2 = -1; if (__pyx_t_4<
0) { __pyx_t_4 += __pyx_bshape_0_x; if (unlikely(__pyx_t_4<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_4 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_ks_10 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_4, __pyx_bstride_0_x));
24: cdef double fs_01 = x[ 3]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":24 * cdef double ks_01 = x[ 1] * cdef double ks_10 = x[ 2] * cdef double fs_01 = x[ 3] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double fs_10 = x[ 4] * cdef double fs_11 = x[ 5] */ __pyx_t_5 = 3; __pyx_t_2 = -1; if (__pyx_t_5<
0) { __pyx_t_5 += __pyx_bshape_0_x; if (unlikely(__pyx_t_5<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_5 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_fs_01 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_5, __pyx_bstride_0_x));
25: cdef double fs_10 = x[ 4]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":25 * cdef double ks_10 = x[ 2] * cdef double fs_01 = x[ 3] * cdef double fs_10 = x[ 4] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double fs_11 = x[ 5] * cdef double k = x[ 6] */ __pyx_t_6 = 4; __pyx_t_2 = -1; if (__pyx_t_6<
0) { __pyx_t_6 += __pyx_bshape_0_x; if (unlikely(__pyx_t_6<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_6 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_fs_10 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_6, __pyx_bstride_0_x));
26: cdef double fs_11 = x[ 5]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":26 * cdef double fs_01 = x[ 3] * cdef double fs_10 = x[ 4] * cdef double fs_11 = x[ 5] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k = x[ 6] * cdef double f = x[ 7] */ __pyx_t_7 = 5; __pyx_t_2 = -1; if (__pyx_t_7<
0) { __pyx_t_7 += __pyx_bshape_0_x; if (unlikely(__pyx_t_7<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_7 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_fs_11 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_7, __pyx_bstride_0_x));
27: cdef double k = x[ 6]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":27 * cdef double fs_10 = x[ 4] * cdef double fs_11 = x[ 5] * cdef double k = x[ 6] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double f = x[ 7] * cdef double s_00 = x[ 8] */ __pyx_t_8 = 6; __pyx_t_2 = -1; if (__pyx_t_8<
0) { __pyx_t_8 += __pyx_bshape_0_x; if (unlikely(__pyx_t_8<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_8 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_8, __pyx_bstride_0_x));
28: cdef double f = x[ 7]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":28 * cdef double fs_11 = x[ 5] * cdef double k = x[ 6] * cdef double f = x[ 7] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double s_00 = x[ 8] * cdef double s_01 = x[ 9] */ __pyx_t_9 = 7; __pyx_t_2 = -1; if (__pyx_t_9<
0) { __pyx_t_9 += __pyx_bshape_0_x; if (unlikely(__pyx_t_9<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_9 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_f = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_9, __pyx_bstride_0_x));
29: cdef double s_00 = x[ 8]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":29 * cdef double k = x[ 6] * cdef double f = x[ 7] * cdef double s_00 = x[ 8] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double s_01 = x[ 9] * cdef double s_10 = x[10] */ __pyx_t_10 = 8; __pyx_t_2 = -1; if (__pyx_t_10<
0) { __pyx_t_10 += __pyx_bshape_0_x; if (unlikely(__pyx_t_10<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_10 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_s_00 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_10, __pyx_bstride_0_x));
30: cdef double s_01 = x[ 9]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":30 * cdef double f = x[ 7] * cdef double s_00 = x[ 8] * cdef double s_01 = x[ 9] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double s_10 = x[10] * cdef double s_11 = x[11] */ __pyx_t_11 = 9; __pyx_t_2 = -1; if (__pyx_t_11<
0) { __pyx_t_11 += __pyx_bshape_0_x; if (unlikely(__pyx_t_11<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_11 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_s_01 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_11, __pyx_bstride_0_x));
31: cdef double s_10 = x[10]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":31 * cdef double s_00 = x[ 8] * cdef double s_01 = x[ 9] * cdef double s_10 = x[10] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double s_11 = x[11] * return np.array( */ __pyx_t_12 = 10; __pyx_t_2 = -1; if (__pyx_t_12<
0) { __pyx_t_12 += __pyx_bshape_0_x; if (unlikely(__pyx_t_12<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_12 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_s_10 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_12, __pyx_bstride_0_x));
32: cdef double s_11 = x[11]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":32 * cdef double s_01 = x[ 9] * cdef double s_10 = x[10] * cdef double s_11 = x[11] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return np.array( * [(-b_00 - c_0001 - c_0010)*ks_00 + a_00*k*s_00, */ __pyx_t_13 = 11; __pyx_t_2 = -1; if (__pyx_t_13<
0) { __pyx_t_13 += __pyx_bshape_0_x; if (unlikely(__pyx_t_13<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_13 >= __pyx_bshape_0_x)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_s_11 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_13, __pyx_bstride_0_x));
33: return np.array(
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":33 * cdef double s_10 = x[10] * cdef double s_11 = x[11] * return np.array( #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* [(-b_00 - c_0001 - c_0010)*ks_00 + a_00*k*s_00, * (-b_01 - c_0111)*ks_01 + a_01*k*s_01, */ __Pyx_XDECREF(((PyObject *)__pyx_r)); __pyx_t_14 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __pyx_t_15 = PyObject_GetAttr(__pyx_t_14, __pyx_n_s__array); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_14); __pyx_t_14 = 0;
34: [(-b_00 - c_0001 - c_0010)*ks_00 + a_00*k*s_00,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":34 * cdef double s_11 = x[11] * return np.array( * [(-b_00 - c_0001 - c_0010)*ks_00 + a_00*k*s_00, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (-b_01 - c_0111)*ks_01 + a_01*k*s_01, * (-b_10 - c_1011)*ks_10 + a_10*k*s_10, */ __pyx_t_14 = PyFloat_FromDouble((((((-__pyx_v_b_00) - __pyx_v_c_0001) - __pyx_v_c_0010) * __pyx_v_ks_00) + ((__pyx_v_a_00 * __pyx_v_k) * __pyx_v_s_00))); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14);
35: (-b_01 - c_0111)*ks_01 + a_01*k*s_01,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":35 * return np.array( * [(-b_00 - c_0001 - c_0010)*ks_00 + a_00*k*s_00, * (-b_01 - c_0111)*ks_01 + a_01*k*s_01, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (-b_10 - c_1011)*ks_10 + a_10*k*s_10, * (-beta_01 - gamma_0100)*fs_01 + alpha_01*f*s_01, */ __pyx_t_16 = PyFloat_FromDouble(((((-__pyx_v_b_01) - __pyx_v_c_0111) * __pyx_v_ks_01) + ((__pyx_v_a_01 * __pyx_v_k) * __pyx_v_s_01))); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16);
36: (-b_10 - c_1011)*ks_10 + a_10*k*s_10,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":36 * [(-b_00 - c_0001 - c_0010)*ks_00 + a_00*k*s_00, * (-b_01 - c_0111)*ks_01 + a_01*k*s_01, * (-b_10 - c_1011)*ks_10 + a_10*k*s_10, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (-beta_01 - gamma_0100)*fs_01 + alpha_01*f*s_01, * (-beta_10 - gamma_1000)*fs_10 + alpha_10*f*s_10, */ __pyx_t_17 = PyFloat_FromDouble(((((-__pyx_v_b_10) - __pyx_v_c_1011) * __pyx_v_ks_10) + ((__pyx_v_a_10 * __pyx_v_k) * __pyx_v_s_10))); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17);
37: (-beta_01 - gamma_0100)*fs_01 + alpha_01*f*s_01,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":37 * (-b_01 - c_0111)*ks_01 + a_01*k*s_01, * (-b_10 - c_1011)*ks_10 + a_10*k*s_10, * (-beta_01 - gamma_0100)*fs_01 + alpha_01*f*s_01, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (-beta_10 - gamma_1000)*fs_10 + alpha_10*f*s_10, * (-beta_11 - gamma_1100 - gamma_1101 - gamma_1110)*fs_11 + alpha_11*f*s_11, */ __pyx_t_18 = PyFloat_FromDouble(((((-__pyx_v_beta_01) - __pyx_v_gamma_0100) * __pyx_v_fs_01) + ((__pyx_v_alpha_01 * __pyx_v_f) * __pyx_v_s_01))); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18);
38: (-beta_10 - gamma_1000)*fs_10 + alpha_10*f*s_10,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":38 * (-b_10 - c_1011)*ks_10 + a_10*k*s_10, * (-beta_01 - gamma_0100)*fs_01 + alpha_01*f*s_01, * (-beta_10 - gamma_1000)*fs_10 + alpha_10*f*s_10, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (-beta_11 - gamma_1100 - gamma_1101 - gamma_1110)*fs_11 + alpha_11*f*s_11, * (b_00 + c_0001 + c_0010)*ks_00 + (b_01 + c_0111)*ks_01 + (b_10 + c_1011)*ks_10 + (-a_00)*k*s_00 + (-a_01)*k*s_01 + (-a_10)*k*s_10, */ __pyx_t_19 = PyFloat_FromDouble(((((-__pyx_v_beta_10) - __pyx_v_gamma_1000) * __pyx_v_fs_10) + ((__pyx_v_alpha_10 * __pyx_v_f) * __pyx_v_s_10))); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19);
39: (-beta_11 - gamma_1100 - gamma_1101 - gamma_1110)*fs_11 + alpha_11*f*s_11,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":39 * (-beta_01 - gamma_0100)*fs_01 + alpha_01*f*s_01, * (-beta_10 - gamma_1000)*fs_10 + alpha_10*f*s_10, * (-beta_11 - gamma_1100 - gamma_1101 - gamma_1110)*fs_11 + alpha_11*f*s_11, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (b_00 + c_0001 + c_0010)*ks_00 + (b_01 + c_0111)*ks_01 + (b_10 + c_1011)*ks_10 + (-a_00)*k*s_00 + (-a_01)*k*s_01 + (-a_10)*k*s_10, * (beta_01 + gamma_0100)*fs_01 + (beta_10 + gamma_1000)*fs_10 + (beta_11 + gamma_1100 + gamma_1101 + gamma_1110)*fs_11 + (-alpha_01)*f*s_01 + (-alpha_10)*f*s_10 + (-alpha_11)*f*s_11, */ __pyx_t_20 = PyFloat_FromDouble(((((((-__pyx_v_beta_11) - __pyx_v_gamma_1100) - __pyx_v_gamma_1101) - __pyx_v_gamma_1110) * __pyx_v_fs_11) + ((__pyx_v_alpha_11 * __pyx_v_f) * __pyx_v_s_11))); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20);
40: (b_00 + c_0001 + c_0010)*ks_00 + (b_01 + c_0111)*ks_01 + (b_10 + c_1011)*ks_10 + (-a_00)*k*s_00 + (-a_01)*k*s_01 + (-a_10)*k*s_10,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":40 * (-beta_10 - gamma_1000)*fs_10 + alpha_10*f*s_10, * (-beta_11 - gamma_1100 - gamma_1101 - gamma_1110)*fs_11 + alpha_11*f*s_11, * (b_00 + c_0001 + c_0010)*ks_00 + (b_01 + c_0111)*ks_01 + (b_10 + c_1011)*ks_10 + (-a_00)*k*s_00 + (-a_01)*k*s_01 + (-a_10)*k*s_10, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* (beta_01 + gamma_0100)*fs_01 + (beta_10 + gamma_1000)*fs_10 + (beta_11 + gamma_1100 + gamma_1101 + gamma_1110)*fs_11 + (-alpha_01)*f*s_01 + (-alpha_10)*f*s_10 + (-alpha_11)*f*s_11, * b_00*ks_00 + gamma_0100*fs_01 + gamma_1000*fs_10 + gamma_1100*fs_11 + (-a_00)*k*s_00, */ __pyx_t_21 = PyFloat_FromDouble(((((((((__pyx_v_b_00 + __pyx_v_c_0001) + __pyx_v_c_0010) * __pyx_v_ks_00) + ((__pyx_v_b_01 + __pyx_v_c_0111) * __pyx_v_ks_01)) + ((__pyx_v_b_10 + __pyx_v_c_1011) * __pyx_v_ks_10)) + (((-__pyx_v_a_00) * __pyx_v_k) * __pyx_v_s_00)) + (((-__pyx_v_a_01) * __pyx_v_k) * __pyx_v_s_01)) + (((-__pyx_v_a_10) * __pyx_v_k) * __pyx_v_s_10))); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_21);
41: (beta_01 + gamma_0100)*fs_01 + (beta_10 + gamma_1000)*fs_10 + (beta_11 + gamma_1100 + gamma_1101 + gamma_1110)*fs_11 + (-alpha_01)*f*s_01 + (-alpha_10)*f*s_10 + (-alpha_11)*f*s_11,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":41 * (-beta_11 - gamma_1100 - gamma_1101 - gamma_1110)*fs_11 + alpha_11*f*s_11, * (b_00 + c_0001 + c_0010)*ks_00 + (b_01 + c_0111)*ks_01 + (b_10 + c_1011)*ks_10 + (-a_00)*k*s_00 + (-a_01)*k*s_01 + (-a_10)*k*s_10, * (beta_01 + gamma_0100)*fs_01 + (beta_10 + gamma_1000)*fs_10 + (beta_11 + gamma_1100 + gamma_1101 + gamma_1110)*fs_11 + (-alpha_01)*f*s_01 + (-alpha_10)*f*s_10 + (-alpha_11)*f*s_11, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* b_00*ks_00 + gamma_0100*fs_01 + gamma_1000*fs_10 + gamma_1100*fs_11 + (-a_00)*k*s_00, * c_0001*ks_00 + b_01*ks_01 + beta_01*fs_01 + gamma_1101*fs_11 + (-a_01)*k*s_01 + (-alpha_01)*f*s_01, */ __pyx_t_22 = PyFloat_FromDouble((((((((__pyx_v_beta_01 + __pyx_v_gamma_0100) * __pyx_v_fs_01) + ((__pyx_v_beta_10 + __pyx_v_gamma_1000) * __pyx_v_fs_10)) + ((((__pyx_v_beta_11 + __pyx_v_gamma_1100) + __pyx_v_gamma_1101) + __pyx_v_gamma_1110) * __pyx_v_fs_11)) + (((-__pyx_v_alpha_01) * __pyx_v_f) * __pyx_v_s_01)) + (((-__pyx_v_alpha_10) * __pyx_v_f) * __pyx_v_s_10)) + (((-__pyx_v_alpha_11) * __pyx_v_f) * __pyx_v_s_11))); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22);
42: b_00*ks_00 + gamma_0100*fs_01 + gamma_1000*fs_10 + gamma_1100*fs_11 + (-a_00)*k*s_00,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":42 * (b_00 + c_0001 + c_0010)*ks_00 + (b_01 + c_0111)*ks_01 + (b_10 + c_1011)*ks_10 + (-a_00)*k*s_00 + (-a_01)*k*s_01 + (-a_10)*k*s_10, * (beta_01 + gamma_0100)*fs_01 + (beta_10 + gamma_1000)*fs_10 + (beta_11 + gamma_1100 + gamma_1101 + gamma_1110)*fs_11 + (-alpha_01)*f*s_01 + (-alpha_10)*f*s_10 + (-alpha_11)*f*s_11, * b_00*ks_00 + gamma_0100*fs_01 + gamma_1000*fs_10 + gamma_1100*fs_11 + (-a_00)*k*s_00, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* c_0001*ks_00 + b_01*ks_01 + beta_01*fs_01 + gamma_1101*fs_11 + (-a_01)*k*s_01 + (-alpha_01)*f*s_01, * c_0010*ks_00 + b_10*ks_10 + beta_10*fs_10 + gamma_1110*fs_11 + (-a_10)*k*s_10 + (-alpha_10)*f*s_10, */ __pyx_t_23 = PyFloat_FromDouble((((((__pyx_v_b_00 * __pyx_v_ks_00) + (__pyx_v_gamma_0100 * __pyx_v_fs_01)) + (__pyx_v_gamma_1000 * __pyx_v_fs_10)) + (__pyx_v_gamma_1100 * __pyx_v_fs_11)) + (((-__pyx_v_a_00) * __pyx_v_k) * __pyx_v_s_00))); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 42; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23);
43: c_0001*ks_00 + b_01*ks_01 + beta_01*fs_01 + gamma_1101*fs_11 + (-a_01)*k*s_01 + (-alpha_01)*f*s_01,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":43 * (beta_01 + gamma_0100)*fs_01 + (beta_10 + gamma_1000)*fs_10 + (beta_11 + gamma_1100 + gamma_1101 + gamma_1110)*fs_11 + (-alpha_01)*f*s_01 + (-alpha_10)*f*s_10 + (-alpha_11)*f*s_11, * b_00*ks_00 + gamma_0100*fs_01 + gamma_1000*fs_10 + gamma_1100*fs_11 + (-a_00)*k*s_00, * c_0001*ks_00 + b_01*ks_01 + beta_01*fs_01 + gamma_1101*fs_11 + (-a_01)*k*s_01 + (-alpha_01)*f*s_01, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* c_0010*ks_00 + b_10*ks_10 + beta_10*fs_10 + gamma_1110*fs_11 + (-a_10)*k*s_10 + (-alpha_10)*f*s_10, * c_0111*ks_01 + c_1011*ks_10 + beta_11*fs_11 + (-alpha_11)*f*s_11]) */ __pyx_t_24 = PyFloat_FromDouble(((((((__pyx_v_c_0001 * __pyx_v_ks_00) + (__pyx_v_b_01 * __pyx_v_ks_01)) + (__pyx_v_beta_01 * __pyx_v_fs_01)) + (__pyx_v_gamma_1101 * __pyx_v_fs_11)) + (((-__pyx_v_a_01) * __pyx_v_k) * __pyx_v_s_01)) + (((-__pyx_v_alpha_01) * __pyx_v_f) * __pyx_v_s_01))); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 43; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24);
44: c_0010*ks_00 + b_10*ks_10 + beta_10*fs_10 + gamma_1110*fs_11 + (-a_10)*k*s_10 + (-alpha_10)*f*s_10,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":44 * b_00*ks_00 + gamma_0100*fs_01 + gamma_1000*fs_10 + gamma_1100*fs_11 + (-a_00)*k*s_00, * c_0001*ks_00 + b_01*ks_01 + beta_01*fs_01 + gamma_1101*fs_11 + (-a_01)*k*s_01 + (-alpha_01)*f*s_01, * c_0010*ks_00 + b_10*ks_10 + beta_10*fs_10 + gamma_1110*fs_11 + (-a_10)*k*s_10 + (-alpha_10)*f*s_10, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* c_0111*ks_01 + c_1011*ks_10 + beta_11*fs_11 + (-alpha_11)*f*s_11]) */ __pyx_t_25 = PyFloat_FromDouble(((((((__pyx_v_c_0010 * __pyx_v_ks_00) + (__pyx_v_b_10 * __pyx_v_ks_10)) + (__pyx_v_beta_10 * __pyx_v_fs_10)) + (__pyx_v_gamma_1110 * __pyx_v_fs_11)) + (((-__pyx_v_a_10) * __pyx_v_k) * __pyx_v_s_10)) + (((-__pyx_v_alpha_10) * __pyx_v_f) * __pyx_v_s_10))); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 44; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25);
45: c_0111*ks_01 + c_1011*ks_10 + beta_11*fs_11 + (-alpha_11)*f*s_11])
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.pyx":45 * c_0001*ks_00 + b_01*ks_01 + beta_01*fs_01 + gamma_1101*fs_11 + (-a_01)*k*s_01 + (-alpha_01)*f*s_01, * c_0010*ks_00 + b_10*ks_10 + beta_10*fs_10 + gamma_1110*fs_11 + (-a_10)*k*s_10 + (-alpha_10)*f*s_10, * c_0111*ks_01 + c_1011*ks_10 + beta_11*fs_11 + (-alpha_11)*f*s_11]) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
*/ __pyx_t_26 = PyFloat_FromDouble(((((__pyx_v_c_0111 * __pyx_v_ks_01) + (__pyx_v_c_1011 * __pyx_v_ks_10)) + (__pyx_v_beta_11 * __pyx_v_fs_11)) + (((-__pyx_v_alpha_11) * __pyx_v_f) * __pyx_v_s_11))); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 45; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_27 = PyList_New(12); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_27)); PyList_SET_ITEM(__pyx_t_27, 0, __pyx_t_14); __Pyx_GIVEREF(__pyx_t_14); PyList_SET_ITEM(__pyx_t_27, 1, __pyx_t_16); __Pyx_GIVEREF(__pyx_t_16); PyList_SET_ITEM(__pyx_t_27, 2, __pyx_t_17); __Pyx_GIVEREF(__pyx_t_17); PyList_SET_ITEM(__pyx_t_27, 3, __pyx_t_18); __Pyx_GIVEREF(__pyx_t_18); PyList_SET_ITEM(__pyx_t_27, 4, __pyx_t_19); __Pyx_GIVEREF(__pyx_t_19); PyList_SET_ITEM(__pyx_t_27, 5, __pyx_t_20); __Pyx_GIVEREF(__pyx_t_20); PyList_SET_ITEM(__pyx_t_27, 6, __pyx_t_21); __Pyx_GIVEREF(__pyx_t_21); PyList_SET_ITEM(__pyx_t_27, 7, __pyx_t_22); __Pyx_GIVEREF(__pyx_t_22); PyList_SET_ITEM(__pyx_t_27, 8, __pyx_t_23); __Pyx_GIVEREF(__pyx_t_23); PyList_SET_ITEM(__pyx_t_27, 9, __pyx_t_24); __Pyx_GIVEREF(__pyx_t_24); PyList_SET_ITEM(__pyx_t_27, 10, __pyx_t_25); __Pyx_GIVEREF(__pyx_t_25); PyList_SET_ITEM(__pyx_t_27, 11, __pyx_t_26); __Pyx_GIVEREF(__pyx_t_26); __pyx_t_14 = 0; __pyx_t_16 = 0; __pyx_t_17 = 0; __pyx_t_18 = 0; __pyx_t_19 = 0; __pyx_t_20 = 0; __pyx_t_21 = 0; __pyx_t_22 = 0; __pyx_t_23 = 0; __pyx_t_24 = 0; __pyx_t_25 = 0; __pyx_t_26 = 0; __pyx_t_26 = PyTuple_New(1); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_26)); PyTuple_SET_ITEM(__pyx_t_26, 0, ((PyObject *)__pyx_t_27)); __Pyx_GIVEREF(((PyObject *)__pyx_t_27)); __pyx_t_27 = 0; __pyx_t_27 = PyObject_Call(__pyx_t_15, ((PyObject *)__pyx_t_26), NULL); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_26)); __pyx_t_26 = 0; if (!(likely(((__pyx_t_27) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_27, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_r = ((PyArrayObject *)__pyx_t_27); __pyx_t_27 = 0; goto __pyx_L0; __pyx_r = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_14); __Pyx_XDECREF(__pyx_t_15); __Pyx_XDECREF(__pyx_t_16); __Pyx_XDECREF(__pyx_t_17); __Pyx_XDECREF(__pyx_t_18); __Pyx_XDECREF(__pyx_t_19); __Pyx_XDECREF(__pyx_t_20); __Pyx_XDECREF(__pyx_t_21); __Pyx_XDECREF(__pyx_t_22); __Pyx_XDECREF(__pyx_t_23); __Pyx_XDECREF(__pyx_t_24); __Pyx_XDECREF(__pyx_t_25); __Pyx_XDECREF(__pyx_t_26); __Pyx_XDECREF(__pyx_t_27); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);} __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage21_spyx_0.ode_DP", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = 0; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XGIVEREF((PyObject *)__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }