Generated by Cython 0.15.1 on Thu May 24 21:36:21 2012
Raw output: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.c
1:
2: include "interrupt.pxi" # ctrl-c interrupt block support
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":2 * * include "interrupt.pxi" # ctrl-c interrupt block support #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* include "stdsage.pxi" # ctrl-c interrupt block support * */ __pyx_t_1 = PyDict_New(); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 2; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); if (PyObject_SetAttr(__pyx_m, __pyx_n_s____test__, ((PyObject *)__pyx_t_1))<
0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 2; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 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_sage79_spyx_0.pyx":6 * * include "cdefs.pxi" * import numpy as np #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cimport numpy as np * from scipy.special import gamma */ __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: from scipy.special import gamma
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":8 * import numpy as np * cimport numpy as np * from scipy.special import gamma #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* * cdouble = np.double */ __pyx_t_1 = PyList_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 8; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); __Pyx_INCREF(((PyObject *)__pyx_n_s__gamma)); PyList_SET_ITEM(__pyx_t_1, 0, ((PyObject *)__pyx_n_s__gamma)); __Pyx_GIVEREF(((PyObject *)__pyx_n_s__gamma)); __pyx_t_2 = __Pyx_Import(((PyObject *)__pyx_n_s_15), ((PyObject *)__pyx_t_1), -1); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 8; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __pyx_t_1 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__gamma); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 8; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (PyObject_SetAttr(__pyx_m, __pyx_n_s__gamma, __pyx_t_1)<
0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 8; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0;
9:
10: cdouble = np.double
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":10 * from scipy.special import gamma * * cdouble = np.double #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* ctypedef np.double_t cdouble_t * */ __pyx_t_2 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_1 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__double); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; if (PyObject_SetAttr(__pyx_m, __pyx_n_s__cdouble, __pyx_t_1)<
0) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 10; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0;
11: ctypedef np.double_t cdouble_t
12:
13: cpdef coeffs_clust(np.ndarray[cdouble_t] a):
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":13 * ctypedef np.double_t cdouble_t * * cpdef coeffs_clust(np.ndarray[cdouble_t] a): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_o = a[0] * cdef double k_c = a[1] */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(PyObject *__pyx_self, PyObject *__pyx_v_a); /*proto*/ static PyObject *__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(PyArrayObject *__pyx_v_a, int __pyx_skip_dispatch) { double __pyx_v_k_o; double __pyx_v_k_c; double __pyx_v_k_u; double __pyx_v_k_s1; double __pyx_v_k_s2; double __pyx_v_k_s3; double __pyx_v_k_l1; double __pyx_v_k_l2; double __pyx_v_k_l3; Py_buffer __pyx_bstruct_a; Py_ssize_t __pyx_bstride_0_a = 0; Py_ssize_t __pyx_bshape_0_a = 0; PyObject *__pyx_r = NULL; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("coeffs_clust"); __pyx_bstruct_a.buf = NULL; { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_a, (PyObject*)__pyx_v_a, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_a = __pyx_bstruct_a.strides[0]; __pyx_bshape_0_a = __pyx_bstruct_a.shape[0]; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":13 * ctypedef np.double_t cdouble_t * * cpdef coeffs_clust(np.ndarray[cdouble_t] a): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_o = a[0] * cdef double k_c = a[1] */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(PyObject *__pyx_self, PyObject *__pyx_v_a); /*proto*/ static char __pyx_doc_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust[] = "File: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx (starting at line 13)"; static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(PyObject *__pyx_self, PyObject *__pyx_v_a) { Py_buffer __pyx_bstruct_a; Py_ssize_t __pyx_bstride_0_a = 0; Py_ssize_t __pyx_bshape_0_a = 0; PyObject *__pyx_r = NULL; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("coeffs_clust"); __pyx_self = __pyx_self; __pyx_bstruct_a.buf = NULL; if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_a), __pyx_ptype_5numpy_ndarray, 1, "a", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L1_error;} { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_a, (PyObject*)__pyx_v_a, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_a = __pyx_bstruct_a.strides[0]; __pyx_bshape_0_a = __pyx_bstruct_a.shape[0]; __Pyx_XDECREF(__pyx_r); __pyx_t_1 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(((PyArrayObject *)__pyx_v_a), 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 13; __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_a); __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);} __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.coeffs_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = NULL; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
14: cdef double k_o = a[0]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":14 * * cpdef coeffs_clust(np.ndarray[cdouble_t] a): * cdef double k_o = a[0] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_c = a[1] * cdef double k_u = a[2] */ __pyx_t_1 = 0; __pyx_t_2 = -1; if (__pyx_t_1<
0) { __pyx_t_1 += __pyx_bshape_0_a; if (unlikely(__pyx_t_1<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_1 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 14; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_o = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_1, __pyx_bstride_0_a));
15: cdef double k_c = a[1]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":15 * cpdef coeffs_clust(np.ndarray[cdouble_t] a): * cdef double k_o = a[0] * cdef double k_c = a[1] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_u = a[2] * cdef double k_s1 = a[3] */ __pyx_t_3 = 1; __pyx_t_2 = -1; if (__pyx_t_3<
0) { __pyx_t_3 += __pyx_bshape_0_a; if (unlikely(__pyx_t_3<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_3 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 15; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_c = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_3, __pyx_bstride_0_a));
16: cdef double k_u = a[2]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":16 * cdef double k_o = a[0] * cdef double k_c = a[1] * cdef double k_u = a[2] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_s1 = a[3] * cdef double k_s2 = a[4] */ __pyx_t_4 = 2; __pyx_t_2 = -1; if (__pyx_t_4<
0) { __pyx_t_4 += __pyx_bshape_0_a; if (unlikely(__pyx_t_4<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_4 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 16; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_u = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_4, __pyx_bstride_0_a));
17: cdef double k_s1 = a[3]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":17 * cdef double k_c = a[1] * cdef double k_u = a[2] * cdef double k_s1 = a[3] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_s2 = a[4] * cdef double k_s3 = a[5] */ __pyx_t_5 = 3; __pyx_t_2 = -1; if (__pyx_t_5<
0) { __pyx_t_5 += __pyx_bshape_0_a; if (unlikely(__pyx_t_5<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_5 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 17; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_s1 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_5, __pyx_bstride_0_a));
18: cdef double k_s2 = a[4]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":18 * cdef double k_u = a[2] * cdef double k_s1 = a[3] * cdef double k_s2 = a[4] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_s3 = a[5] * cdef double k_l1 = a[6] */ __pyx_t_6 = 4; __pyx_t_2 = -1; if (__pyx_t_6<
0) { __pyx_t_6 += __pyx_bshape_0_a; if (unlikely(__pyx_t_6<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_6 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 18; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_s2 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_6, __pyx_bstride_0_a));
19: cdef double k_s3 = a[5]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":19 * cdef double k_s1 = a[3] * cdef double k_s2 = a[4] * cdef double k_s3 = a[5] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_l1 = a[6] * cdef double k_l2 = a[7] */ __pyx_t_7 = 5; __pyx_t_2 = -1; if (__pyx_t_7<
0) { __pyx_t_7 += __pyx_bshape_0_a; if (unlikely(__pyx_t_7<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_7 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 19; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_s3 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_7, __pyx_bstride_0_a));
20: cdef double k_l1 = a[6]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":20 * cdef double k_s2 = a[4] * cdef double k_s3 = a[5] * cdef double k_l1 = a[6] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_l2 = a[7] * cdef double k_l3 = a[8] */ __pyx_t_8 = 6; __pyx_t_2 = -1; if (__pyx_t_8<
0) { __pyx_t_8 += __pyx_bshape_0_a; if (unlikely(__pyx_t_8<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_8 >= __pyx_bshape_0_a)) __pyx_t_2 = 0; if (unlikely(__pyx_t_2 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 20; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_k_l1 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_8, __pyx_bstride_0_a));
21: cdef double k_l2 = a[7]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":21 * cdef double k_s3 = a[5] * cdef double k_l1 = a[6] * cdef double k_l2 = a[7] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double k_l3 = a[8] * return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 + */ __pyx_t_9 = 7; __pyx_t_2 = -1; if (__pyx_t_9<
0) { __pyx_t_9 += __pyx_bshape_0_a; if (unlikely(__pyx_t_9<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_9 >= __pyx_bshape_0_a)) __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_k_l2 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_9, __pyx_bstride_0_a));
22: cdef double k_l3 = a[8]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":22 * cdef double k_l1 = a[6] * cdef double k_l2 = a[7] * cdef double k_l3 = a[8] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 + * 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + */ __pyx_t_10 = 8; __pyx_t_2 = -1; if (__pyx_t_10<
0) { __pyx_t_10 += __pyx_bshape_0_a; if (unlikely(__pyx_t_10<
0)) __pyx_t_2 = 0; } else if (unlikely(__pyx_t_10 >= __pyx_bshape_0_a)) __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_k_l3 = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_a.buf, __pyx_t_10, __pyx_bstride_0_a));
23: return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 +
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":23 * cdef double k_l2 = a[7] * cdef double k_l3 = a[8] * return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + * 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - */ __Pyx_XDECREF(__pyx_r); __pyx_t_11 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_12 = PyObject_GetAttr(__pyx_t_11, __pyx_n_s__array); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_12); __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_13 = PyNumber_Multiply(__pyx_int_78364164096, __pyx_t_11); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_14 = PyNumber_Multiply(__pyx_t_13, __pyx_t_11); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_13 = PyNumber_Multiply(__pyx_int_neg_195910410240, __pyx_t_11); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_15 = PyNumber_Multiply(__pyx_t_13, __pyx_t_11); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0;
24: 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 +
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":24 * cdef double k_l3 = a[8] * return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 + * 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - * 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, */ __pyx_t_11 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_13 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_11); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_16 = PyNumber_Multiply(__pyx_t_13, __pyx_t_11); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_13 = PyNumber_Multiply(__pyx_t_16, __pyx_t_11); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __Pyx_DECREF(__pyx_t_11); __pyx_t_11 = 0; __pyx_t_11 = PyNumber_Add(__pyx_t_15, __pyx_t_13); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __pyx_t_13 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_15 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_13); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __pyx_t_13 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_16 = PyNumber_Multiply(__pyx_t_15, __pyx_t_13); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0;
25: 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 -
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":25 * return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 + * 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + * 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + */ __pyx_t_13 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_15 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_13); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __pyx_t_13 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_17 = PyNumber_Multiply(__pyx_t_15, __pyx_t_13); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __pyx_t_13 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_15 = PyNumber_Multiply(__pyx_t_17, __pyx_t_13); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __Pyx_DECREF(__pyx_t_13); __pyx_t_13 = 0; __pyx_t_13 = PyNumber_Add(__pyx_t_16, __pyx_t_15); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 24; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":26 * 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + * 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - * 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + */ __pyx_t_15 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_16 = PyNumber_Multiply(__pyx_int_169789022208, __pyx_t_15); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":25 * return np.array([78364164096*k_o**3*k_l3, -195910410240*k_o**3*k_l3 + * 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + * 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + */ __pyx_t_15 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_17 = PyNumber_Multiply(__pyx_t_16, __pyx_t_15); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0;
26: 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":26 * 39182082048*k_o**2*k_c*k_l3, 39182082048*k_o**3*k_l2 + * 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - * 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + */ __pyx_t_15 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_16 = PyNumber_Multiply(__pyx_int_52242776064, __pyx_t_15); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __pyx_t_15 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_18 = PyNumber_Multiply(__pyx_t_16, __pyx_t_15); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __pyx_t_15 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_16 = PyNumber_Multiply(__pyx_t_18, __pyx_t_15); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __pyx_t_15 = PyNumber_Subtract(__pyx_t_17, __pyx_t_16); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 25; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __pyx_t_16 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __pyx_t_17 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_16); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __pyx_t_16 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __pyx_t_18 = PyNumber_Multiply(__pyx_t_17, __pyx_t_16); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __pyx_t_16 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __pyx_t_17 = PyNumber_Multiply(__pyx_t_18, __pyx_t_16); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_16); __pyx_t_16 = 0; __pyx_t_16 = PyNumber_Add(__pyx_t_15, __pyx_t_17); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 26; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0;
27: -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 +
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":27 * 39182082048*k_o**2*k_c*k_l2, 169789022208*k_o**3*k_l3 - * 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + * 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, */ __pyx_t_17 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_15 = PyNumber_Multiply(__pyx_int_neg_65303470080, __pyx_t_17); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_18 = PyNumber_Multiply(__pyx_t_15, __pyx_t_17); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_15 = PyNumber_Multiply(__pyx_int_52242776064, __pyx_t_17); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_19 = PyNumber_Multiply(__pyx_t_15, __pyx_t_17); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_15 = PyNumber_Multiply(__pyx_t_19, __pyx_t_17); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyNumber_Subtract(__pyx_t_18, __pyx_t_15); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0;
28: 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 +
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":28 * 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, * 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - */ __pyx_t_15 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_18 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_15); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __pyx_t_15 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_19 = PyNumber_Multiply(__pyx_t_18, __pyx_t_15); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __pyx_t_15 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_18 = PyNumber_Multiply(__pyx_t_19, __pyx_t_15); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_15); __pyx_t_15 = 0; __pyx_t_15 = PyNumber_Add(__pyx_t_17, __pyx_t_18); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 27; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":29 * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + * 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - * 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3, */ __pyx_t_18 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_17 = PyNumber_Multiply(__pyx_int_neg_52242776064, __pyx_t_18); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":28 * 52242776064*k_o**2*k_c*k_l3 + 13060694016*k_o*k_c**2*k_l3, * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, * 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - */ __pyx_t_18 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_19 = PyNumber_Multiply(__pyx_t_17, __pyx_t_18); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0;
29: 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":29 * -65303470080*k_o**3*k_l2 - 52242776064*k_o**2*k_c*k_l2 + * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + * 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - * 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3, */ __pyx_t_18 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_17 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_18); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __pyx_t_18 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_20 = PyNumber_Multiply(__pyx_t_17, __pyx_t_18); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __pyx_t_18 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_17 = PyNumber_Multiply(__pyx_t_20, __pyx_t_18); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __pyx_t_18 = PyNumber_Add(__pyx_t_19, __pyx_t_17); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 28; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_19 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_17); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_20 = PyNumber_Multiply(__pyx_t_19, __pyx_t_17); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyFloat_FromDouble(__pyx_v_k_l3); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_19 = PyNumber_Multiply(__pyx_t_20, __pyx_t_17); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __Pyx_DECREF(__pyx_t_17); __pyx_t_17 = 0; __pyx_t_17 = PyNumber_Subtract(__pyx_t_18, __pyx_t_19); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 29; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0;
30: 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 -
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":30 * 13060694016*k_o*k_c**2*k_l2, -52242776064*k_o**3*k_l3 + * 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, * 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3, * -195910410240*k_o**3*k_s3 + 39182082048*k_o**2*k_c*k_s3, */ __pyx_t_19 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_18 = PyNumber_Multiply(__pyx_int_26121388032, __pyx_t_19); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_20 = PyNumber_Multiply(__pyx_t_18, __pyx_t_19); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_18 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_19); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_21 = PyNumber_Multiply(__pyx_t_18, __pyx_t_19); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_21); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_18 = PyNumber_Multiply(__pyx_t_21, __pyx_t_19); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_21); __pyx_t_21 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyNumber_Add(__pyx_t_20, __pyx_t_18); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0;
31: 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":31 * 13060694016*k_o**2*k_c*k_l3 - 13060694016*k_o*k_c**2*k_l3, * 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - * 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -195910410240*k_o**3*k_s3 + 39182082048*k_o**2*k_c*k_s3, * 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2, */ __pyx_t_18 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_20 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_18); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __pyx_t_18 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_21 = PyNumber_Multiply(__pyx_t_20, __pyx_t_18); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_21); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __pyx_t_18 = PyFloat_FromDouble(__pyx_v_k_l2); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_20 = PyNumber_Multiply(__pyx_t_21, __pyx_t_18); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __Pyx_DECREF(__pyx_t_21); __pyx_t_21 = 0; __Pyx_DECREF(__pyx_t_18); __pyx_t_18 = 0; __pyx_t_18 = PyNumber_Subtract(__pyx_t_19, __pyx_t_20); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 30; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_19 = PyNumber_Multiply(__pyx_int_78364164096, __pyx_t_20); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_21 = PyNumber_Multiply(__pyx_t_19, __pyx_t_20); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 31; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_21); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0;
32: -195910410240*k_o**3*k_s3 + 39182082048*k_o**2*k_c*k_s3,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":32 * 26121388032*k_o**3*k_l2 + 13060694016*k_o**2*k_c*k_l2 - * 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3, * -195910410240*k_o**3*k_s3 + 39182082048*k_o**2*k_c*k_s3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2, * 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + */ __pyx_t_20 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_19 = PyNumber_Multiply(__pyx_int_neg_195910410240, __pyx_t_20); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_22 = PyNumber_Multiply(__pyx_t_19, __pyx_t_20); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_19 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_20); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_23 = PyNumber_Multiply(__pyx_t_19, __pyx_t_20); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_19 = PyNumber_Multiply(__pyx_t_23, __pyx_t_20); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_t_20 = PyNumber_Add(__pyx_t_22, __pyx_t_19); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 32; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0;
33: 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":33 * 13060694016*k_o*k_c**2*k_l2, 78364164096*k_o**3*k_s3, * -195910410240*k_o**3*k_s3 + 39182082048*k_o**2*k_c*k_s3, * 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + * 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - */ __pyx_t_19 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_22 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_19); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_23 = PyNumber_Multiply(__pyx_t_22, __pyx_t_19); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_22 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_19); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_24 = PyNumber_Multiply(__pyx_t_22, __pyx_t_19); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_22 = PyNumber_Multiply(__pyx_t_24, __pyx_t_19); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __Pyx_DECREF(__pyx_t_19); __pyx_t_19 = 0; __pyx_t_19 = PyNumber_Add(__pyx_t_23, __pyx_t_22); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 33; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0;
34: 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 +
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":34 * -195910410240*k_o**3*k_s3 + 39182082048*k_o**2*k_c*k_s3, * 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2, * 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - * 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, */ __pyx_t_22 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_23 = PyNumber_Multiply(__pyx_int_169789022208, __pyx_t_22); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_24 = PyNumber_Multiply(__pyx_t_23, __pyx_t_22); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_23 = PyNumber_Multiply(__pyx_int_52242776064, __pyx_t_22); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_25 = PyNumber_Multiply(__pyx_t_23, __pyx_t_22); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_23 = PyNumber_Multiply(__pyx_t_25, __pyx_t_22); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyNumber_Subtract(__pyx_t_24, __pyx_t_23); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0;
35: 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 -
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":35 * 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2, * 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + * 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - */ __pyx_t_23 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __pyx_t_24 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_23); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __pyx_t_23 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __pyx_t_25 = PyNumber_Multiply(__pyx_t_24, __pyx_t_23); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __pyx_t_23 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __pyx_t_24 = PyNumber_Multiply(__pyx_t_25, __pyx_t_23); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_23); __pyx_t_23 = 0; __pyx_t_23 = PyNumber_Add(__pyx_t_22, __pyx_t_24); if (unlikely(!__pyx_t_23)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 34; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_23); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":36 * 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + * 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - * 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + */ __pyx_t_24 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_22 = PyNumber_Multiply(__pyx_int_neg_65303470080, __pyx_t_24); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":35 * 39182082048*k_o**3*k_s2 + 39182082048*k_o**2*k_c*k_s2, * 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + * 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - */ __pyx_t_24 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_25 = PyNumber_Multiply(__pyx_t_22, __pyx_t_24); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0;
36: 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":36 * 169789022208*k_o**3*k_s3 - 52242776064*k_o**2*k_c*k_s3 + * 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - * 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + */ __pyx_t_24 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_22 = PyNumber_Multiply(__pyx_int_52242776064, __pyx_t_24); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __pyx_t_24 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_26 = PyNumber_Multiply(__pyx_t_22, __pyx_t_24); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __pyx_t_24 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_22 = PyNumber_Multiply(__pyx_t_26, __pyx_t_24); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __pyx_t_24 = PyNumber_Subtract(__pyx_t_25, __pyx_t_22); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 35; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_25 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_22); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_26 = PyNumber_Multiply(__pyx_t_25, __pyx_t_22); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __pyx_t_25 = PyNumber_Multiply(__pyx_t_26, __pyx_t_22); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_22); __pyx_t_22 = 0; __pyx_t_22 = PyNumber_Add(__pyx_t_24, __pyx_t_25); if (unlikely(!__pyx_t_22)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 36; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_22); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0;
37: -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 -
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":37 * 13060694016*k_o*k_c**2*k_s3, -65303470080*k_o**3*k_s2 - * 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + * 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, */ __pyx_t_25 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_24 = PyNumber_Multiply(__pyx_int_neg_52242776064, __pyx_t_25); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_26 = PyNumber_Multiply(__pyx_t_24, __pyx_t_25); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_24 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_25); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_27 = PyNumber_Multiply(__pyx_t_24, __pyx_t_25); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_24 = PyNumber_Multiply(__pyx_t_27, __pyx_t_25); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyNumber_Add(__pyx_t_26, __pyx_t_24); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0;
38: 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 +
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":38 * 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, * -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u - */ __pyx_t_24 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_26 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_24); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __pyx_t_24 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_27 = PyNumber_Multiply(__pyx_t_26, __pyx_t_24); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __pyx_t_24 = PyFloat_FromDouble(__pyx_v_k_s3); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __pyx_t_26 = PyNumber_Multiply(__pyx_t_27, __pyx_t_24); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __Pyx_DECREF(__pyx_t_24); __pyx_t_24 = 0; __pyx_t_24 = PyNumber_Subtract(__pyx_t_25, __pyx_t_26); if (unlikely(!__pyx_t_24)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 37; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_24); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":39 * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + * 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u - * 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u]) */ __pyx_t_26 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_25 = PyNumber_Multiply(__pyx_int_26121388032, __pyx_t_26); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":38 * 52242776064*k_o**2*k_c*k_s2 + 13060694016*k_o*k_c**2*k_s2, * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, * -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u - */ __pyx_t_26 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_27 = PyNumber_Multiply(__pyx_t_25, __pyx_t_26); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0;
39: 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":39 * -52242776064*k_o**3*k_s3 + 13060694016*k_o**2*k_c*k_s3 - * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + * 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u - * 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u]) */ __pyx_t_26 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_25 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_26); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __pyx_t_26 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_28 = PyNumber_Multiply(__pyx_t_25, __pyx_t_26); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __pyx_t_26 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_25 = PyNumber_Multiply(__pyx_t_28, __pyx_t_26); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __pyx_t_26 = PyNumber_Add(__pyx_t_27, __pyx_t_25); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 38; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_27 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_25); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_28 = PyNumber_Multiply(__pyx_t_27, __pyx_t_25); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyFloat_FromDouble(__pyx_v_k_s2); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __pyx_t_27 = PyNumber_Multiply(__pyx_t_28, __pyx_t_25); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __Pyx_DECREF(__pyx_t_25); __pyx_t_25 = 0; __pyx_t_25 = PyNumber_Subtract(__pyx_t_26, __pyx_t_27); if (unlikely(!__pyx_t_25)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 39; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_25); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0;
40: -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u -
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":40 * 13060694016*k_o*k_c**2*k_s3, 26121388032*k_o**3*k_s2 + * 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, * -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u - #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u]) * */ __pyx_t_27 = PyFloat_FromDouble(pow(__pyx_v_k_o, 3.0)); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __pyx_t_26 = PyNumber_Multiply(__pyx_int_neg_13060694016, __pyx_t_27); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __pyx_t_27 = PyFloat_FromDouble(__pyx_v_k_u); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __pyx_t_28 = PyNumber_Multiply(__pyx_t_26, __pyx_t_27); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __pyx_t_27 = PyFloat_FromDouble(pow(__pyx_v_k_o, 2.0)); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __pyx_t_26 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_27); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __pyx_t_27 = PyFloat_FromDouble(__pyx_v_k_c); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __pyx_t_29 = PyNumber_Multiply(__pyx_t_26, __pyx_t_27); if (unlikely(!__pyx_t_29)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_29); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __pyx_t_27 = PyFloat_FromDouble(__pyx_v_k_u); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __pyx_t_26 = PyNumber_Multiply(__pyx_t_29, __pyx_t_27); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_29); __pyx_t_29 = 0; __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __pyx_t_27 = PyNumber_Subtract(__pyx_t_28, __pyx_t_26); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0;
41: 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u])
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":41 * 13060694016*k_o**2*k_c*k_s2 - 13060694016*k_o*k_c**2*k_s2, * -13060694016*k_o**3*k_u - 39182082048*k_o**2*k_c*k_u - * 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u]) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* * cpdef monoms_clust(np.ndarray[cdouble_t] x): */ __pyx_t_26 = PyFloat_FromDouble(__pyx_v_k_o); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_28 = PyNumber_Multiply(__pyx_int_39182082048, __pyx_t_26); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __pyx_t_26 = PyFloat_FromDouble(pow(__pyx_v_k_c, 2.0)); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_29 = PyNumber_Multiply(__pyx_t_28, __pyx_t_26); if (unlikely(!__pyx_t_29)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_29); __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __pyx_t_26 = PyFloat_FromDouble(__pyx_v_k_u); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __pyx_t_28 = PyNumber_Multiply(__pyx_t_29, __pyx_t_26); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __Pyx_DECREF(__pyx_t_29); __pyx_t_29 = 0; __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __pyx_t_26 = PyNumber_Subtract(__pyx_t_27, __pyx_t_28); if (unlikely(!__pyx_t_26)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 40; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_26); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __pyx_t_28 = PyFloat_FromDouble(pow(__pyx_v_k_c, 3.0)); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __pyx_t_27 = PyNumber_Multiply(__pyx_int_13060694016, __pyx_t_28); if (unlikely(!__pyx_t_27)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_27); __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __pyx_t_28 = PyFloat_FromDouble(__pyx_v_k_u); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __pyx_t_29 = PyNumber_Multiply(__pyx_t_27, __pyx_t_28); if (unlikely(!__pyx_t_29)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_29); __Pyx_DECREF(__pyx_t_27); __pyx_t_27 = 0; __Pyx_DECREF(__pyx_t_28); __pyx_t_28 = 0; __pyx_t_28 = PyNumber_Subtract(__pyx_t_26, __pyx_t_29); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 41; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_28); __Pyx_DECREF(__pyx_t_26); __pyx_t_26 = 0; __Pyx_DECREF(__pyx_t_29); __pyx_t_29 = 0; __pyx_t_29 = PyList_New(15); if (unlikely(!__pyx_t_29)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_29)); PyList_SET_ITEM(__pyx_t_29, 0, __pyx_t_14); __Pyx_GIVEREF(__pyx_t_14); PyList_SET_ITEM(__pyx_t_29, 1, __pyx_t_11); __Pyx_GIVEREF(__pyx_t_11); PyList_SET_ITEM(__pyx_t_29, 2, __pyx_t_13); __Pyx_GIVEREF(__pyx_t_13); PyList_SET_ITEM(__pyx_t_29, 3, __pyx_t_16); __Pyx_GIVEREF(__pyx_t_16); PyList_SET_ITEM(__pyx_t_29, 4, __pyx_t_15); __Pyx_GIVEREF(__pyx_t_15); PyList_SET_ITEM(__pyx_t_29, 5, __pyx_t_17); __Pyx_GIVEREF(__pyx_t_17); PyList_SET_ITEM(__pyx_t_29, 6, __pyx_t_18); __Pyx_GIVEREF(__pyx_t_18); PyList_SET_ITEM(__pyx_t_29, 7, __pyx_t_21); __Pyx_GIVEREF(__pyx_t_21); PyList_SET_ITEM(__pyx_t_29, 8, __pyx_t_20); __Pyx_GIVEREF(__pyx_t_20); PyList_SET_ITEM(__pyx_t_29, 9, __pyx_t_19); __Pyx_GIVEREF(__pyx_t_19); PyList_SET_ITEM(__pyx_t_29, 10, __pyx_t_23); __Pyx_GIVEREF(__pyx_t_23); PyList_SET_ITEM(__pyx_t_29, 11, __pyx_t_22); __Pyx_GIVEREF(__pyx_t_22); PyList_SET_ITEM(__pyx_t_29, 12, __pyx_t_24); __Pyx_GIVEREF(__pyx_t_24); PyList_SET_ITEM(__pyx_t_29, 13, __pyx_t_25); __Pyx_GIVEREF(__pyx_t_25); PyList_SET_ITEM(__pyx_t_29, 14, __pyx_t_28); __Pyx_GIVEREF(__pyx_t_28); __pyx_t_14 = 0; __pyx_t_11 = 0; __pyx_t_13 = 0; __pyx_t_16 = 0; __pyx_t_15 = 0; __pyx_t_17 = 0; __pyx_t_18 = 0; __pyx_t_21 = 0; __pyx_t_20 = 0; __pyx_t_19 = 0; __pyx_t_23 = 0; __pyx_t_22 = 0; __pyx_t_24 = 0; __pyx_t_25 = 0; __pyx_t_28 = 0; __pyx_t_28 = PyTuple_New(1); if (unlikely(!__pyx_t_28)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_28)); PyTuple_SET_ITEM(__pyx_t_28, 0, ((PyObject *)__pyx_t_29)); __Pyx_GIVEREF(((PyObject *)__pyx_t_29)); __pyx_t_29 = 0; __pyx_t_29 = PyObject_Call(__pyx_t_12, ((PyObject *)__pyx_t_28), NULL); if (unlikely(!__pyx_t_29)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 23; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_29); __Pyx_DECREF(__pyx_t_12); __pyx_t_12 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_28)); __pyx_t_28 = 0; __pyx_r = __pyx_t_29; __pyx_t_29 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_11); __Pyx_XDECREF(__pyx_t_12); __Pyx_XDECREF(__pyx_t_13); __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); __Pyx_XDECREF(__pyx_t_28); __Pyx_XDECREF(__pyx_t_29); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_ErrRestore(__pyx_type, __pyx_value, __pyx_tb);} __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.coeffs_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = 0; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
42:
43: cpdef monoms_clust(np.ndarray[cdouble_t] x):
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":43 * 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u]) * * cpdef monoms_clust(np.ndarray[cdouble_t] x): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double Lambda = x[0] * cdef double rho = x[1] */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_1monoms_clust(PyObject *__pyx_self, PyObject *__pyx_v_x); /*proto*/ static PyObject *__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_monoms_clust(PyArrayObject *__pyx_v_x, int __pyx_skip_dispatch) { double __pyx_v_Lambda; double __pyx_v_rho; double __pyx_v_zeta; 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; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("monoms_clust"); __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 43; __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_sage79_spyx_0.pyx":43 * 39182082048*k_o*k_c**2*k_u - 13060694016*k_c**3*k_u]) * * cpdef monoms_clust(np.ndarray[cdouble_t] x): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double Lambda = x[0] * cdef double rho = x[1] */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_1monoms_clust(PyObject *__pyx_self, PyObject *__pyx_v_x); /*proto*/ static char __pyx_doc_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_1monoms_clust[] = "File: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx (starting at line 43)"; static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_1monoms_clust(PyObject *__pyx_self, PyObject *__pyx_v_x) { 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; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("monoms_clust"); __pyx_self = __pyx_self; __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 = 43; __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 43; __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 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_monoms_clust(((PyArrayObject *)__pyx_v_x), 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 43; __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_sage79_spyx_0.monoms_clust", __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; }
44: cdef double Lambda = x[0]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":44 * * cpdef monoms_clust(np.ndarray[cdouble_t] x): * cdef double Lambda = x[0] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double rho = x[1] * cdef double zeta = 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 = 44; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_Lambda = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_1, __pyx_bstride_0_x));
45: cdef double rho = x[1]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":45 * cpdef monoms_clust(np.ndarray[cdouble_t] x): * cdef double Lambda = x[0] * cdef double rho = x[1] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double zeta = x[2] * return np.array([Lambda*rho**3, Lambda*rho**2*zeta, Lambda*rho**2, Lambda*rho*zeta**2, */ __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 = 45; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_rho = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_3, __pyx_bstride_0_x));
46: cdef double zeta = x[2]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":46 * cdef double Lambda = x[0] * cdef double rho = x[1] * cdef double zeta = x[2] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return np.array([Lambda*rho**3, Lambda*rho**2*zeta, Lambda*rho**2, Lambda*rho*zeta**2, * Lambda*rho*zeta, Lambda*zeta**3, Lambda*zeta**2, rho**3, rho**2*zeta, rho**2, */ __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 = 46; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_zeta = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_4, __pyx_bstride_0_x));
47: return np.array([Lambda*rho**3, Lambda*rho**2*zeta, Lambda*rho**2, Lambda*rho*zeta**2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":47 * cdef double rho = x[1] * cdef double zeta = x[2] * return np.array([Lambda*rho**3, Lambda*rho**2*zeta, Lambda*rho**2, Lambda*rho*zeta**2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* Lambda*rho*zeta, Lambda*zeta**3, Lambda*zeta**2, rho**3, rho**2*zeta, rho**2, * rho*zeta**2, rho*zeta, zeta**3, zeta**2, zeta]) */ __Pyx_XDECREF(__pyx_r); __pyx_t_5 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_5)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_5); __pyx_t_6 = PyObject_GetAttr(__pyx_t_5, __pyx_n_s__array); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_5); __pyx_t_5 = 0; __pyx_t_5 = PyFloat_FromDouble((__pyx_v_Lambda * pow(__pyx_v_rho, 3.0))); if (unlikely(!__pyx_t_5)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_5); __pyx_t_7 = PyFloat_FromDouble(((__pyx_v_Lambda * pow(__pyx_v_rho, 2.0)) * __pyx_v_zeta)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __pyx_t_8 = PyFloat_FromDouble((__pyx_v_Lambda * pow(__pyx_v_rho, 2.0))); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_9 = PyFloat_FromDouble(((__pyx_v_Lambda * __pyx_v_rho) * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_9);
48: Lambda*rho*zeta, Lambda*zeta**3, Lambda*zeta**2, rho**3, rho**2*zeta, rho**2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":48 * cdef double zeta = x[2] * return np.array([Lambda*rho**3, Lambda*rho**2*zeta, Lambda*rho**2, Lambda*rho*zeta**2, * Lambda*rho*zeta, Lambda*zeta**3, Lambda*zeta**2, rho**3, rho**2*zeta, rho**2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* rho*zeta**2, rho*zeta, zeta**3, zeta**2, zeta]) * */ __pyx_t_10 = PyFloat_FromDouble(((__pyx_v_Lambda * __pyx_v_rho) * __pyx_v_zeta)); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_10); __pyx_t_11 = PyFloat_FromDouble((__pyx_v_Lambda * pow(__pyx_v_zeta, 3.0))); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11); __pyx_t_12 = PyFloat_FromDouble((__pyx_v_Lambda * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_12); __pyx_t_13 = PyFloat_FromDouble(pow(__pyx_v_rho, 3.0)); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_14 = PyFloat_FromDouble((pow(__pyx_v_rho, 2.0) * __pyx_v_zeta)); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __pyx_t_15 = PyFloat_FromDouble(pow(__pyx_v_rho, 2.0)); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 48; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15);
49: rho*zeta**2, rho*zeta, zeta**3, zeta**2, zeta])
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":49 * return np.array([Lambda*rho**3, Lambda*rho**2*zeta, Lambda*rho**2, Lambda*rho*zeta**2, * Lambda*rho*zeta, Lambda*zeta**3, Lambda*zeta**2, rho**3, rho**2*zeta, rho**2, * rho*zeta**2, rho*zeta, zeta**3, zeta**2, zeta]) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* * cpdef dxi_clust(np.ndarray[cdouble_t] x, double eps): */ __pyx_t_16 = PyFloat_FromDouble((__pyx_v_rho * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 49; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16); __pyx_t_17 = PyFloat_FromDouble((__pyx_v_rho * __pyx_v_zeta)); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 49; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_18 = PyFloat_FromDouble(pow(__pyx_v_zeta, 3.0)); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 49; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_19 = PyFloat_FromDouble(pow(__pyx_v_zeta, 2.0)); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 49; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_20 = PyFloat_FromDouble(__pyx_v_zeta); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 49; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_21 = PyList_New(15); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_21)); PyList_SET_ITEM(__pyx_t_21, 0, __pyx_t_5); __Pyx_GIVEREF(__pyx_t_5); PyList_SET_ITEM(__pyx_t_21, 1, __pyx_t_7); __Pyx_GIVEREF(__pyx_t_7); PyList_SET_ITEM(__pyx_t_21, 2, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); PyList_SET_ITEM(__pyx_t_21, 3, __pyx_t_9); __Pyx_GIVEREF(__pyx_t_9); PyList_SET_ITEM(__pyx_t_21, 4, __pyx_t_10); __Pyx_GIVEREF(__pyx_t_10); PyList_SET_ITEM(__pyx_t_21, 5, __pyx_t_11); __Pyx_GIVEREF(__pyx_t_11); PyList_SET_ITEM(__pyx_t_21, 6, __pyx_t_12); __Pyx_GIVEREF(__pyx_t_12); PyList_SET_ITEM(__pyx_t_21, 7, __pyx_t_13); __Pyx_GIVEREF(__pyx_t_13); PyList_SET_ITEM(__pyx_t_21, 8, __pyx_t_14); __Pyx_GIVEREF(__pyx_t_14); PyList_SET_ITEM(__pyx_t_21, 9, __pyx_t_15); __Pyx_GIVEREF(__pyx_t_15); PyList_SET_ITEM(__pyx_t_21, 10, __pyx_t_16); __Pyx_GIVEREF(__pyx_t_16); PyList_SET_ITEM(__pyx_t_21, 11, __pyx_t_17); __Pyx_GIVEREF(__pyx_t_17); PyList_SET_ITEM(__pyx_t_21, 12, __pyx_t_18); __Pyx_GIVEREF(__pyx_t_18); PyList_SET_ITEM(__pyx_t_21, 13, __pyx_t_19); __Pyx_GIVEREF(__pyx_t_19); PyList_SET_ITEM(__pyx_t_21, 14, __pyx_t_20); __Pyx_GIVEREF(__pyx_t_20); __pyx_t_5 = 0; __pyx_t_7 = 0; __pyx_t_8 = 0; __pyx_t_9 = 0; __pyx_t_10 = 0; __pyx_t_11 = 0; __pyx_t_12 = 0; __pyx_t_13 = 0; __pyx_t_14 = 0; __pyx_t_15 = 0; __pyx_t_16 = 0; __pyx_t_17 = 0; __pyx_t_18 = 0; __pyx_t_19 = 0; __pyx_t_20 = 0; __pyx_t_20 = PyTuple_New(1); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_20)); PyTuple_SET_ITEM(__pyx_t_20, 0, ((PyObject *)__pyx_t_21)); __Pyx_GIVEREF(((PyObject *)__pyx_t_21)); __pyx_t_21 = 0; __pyx_t_21 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_20), NULL); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 47; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_21); __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_20)); __pyx_t_20 = 0; __pyx_r = __pyx_t_21; __pyx_t_21 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_5); __Pyx_XDECREF(__pyx_t_6); __Pyx_XDECREF(__pyx_t_7); __Pyx_XDECREF(__pyx_t_8); __Pyx_XDECREF(__pyx_t_9); __Pyx_XDECREF(__pyx_t_10); __Pyx_XDECREF(__pyx_t_11); __Pyx_XDECREF(__pyx_t_12); __Pyx_XDECREF(__pyx_t_13); __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); { 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_sage79_spyx_0.monoms_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = 0; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
50:
51: cpdef dxi_clust(np.ndarray[cdouble_t] x, double eps):
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":51 * rho*zeta**2, rho*zeta, zeta**3, zeta**2, zeta]) * * cpdef dxi_clust(np.ndarray[cdouble_t] x, double eps): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double Lambda = x[0] * cdef double rho = x[1] */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_2dxi_clust(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static PyObject *__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_dxi_clust(PyArrayObject *__pyx_v_x, double __pyx_v_eps, int __pyx_skip_dispatch) { double __pyx_v_Lambda; double __pyx_v_rho; double __pyx_v_zeta; 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; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("dxi_clust"); __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __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_sage79_spyx_0.pyx":51 * rho*zeta**2, rho*zeta, zeta**3, zeta**2, zeta]) * * cpdef dxi_clust(np.ndarray[cdouble_t] x, double eps): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double Lambda = x[0] * cdef double rho = x[1] */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_2dxi_clust(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_sage79_spyx_0_2dxi_clust[] = "File: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx (starting at line 51)"; static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_2dxi_clust(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) { PyArrayObject *__pyx_v_x = 0; double __pyx_v_eps; 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__eps,0}; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("dxi_clust"); __pyx_self = __pyx_self; { PyObject* values[2] = {0,0}; if (unlikely(__pyx_kwds)) { Py_ssize_t kw_args; switch (PyTuple_GET_SIZE(__pyx_args)) { 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__eps); if (likely(values[1])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("dxi_clust", 1, 2, 2, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __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), "dxi_clust")<
0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } } else if (PyTuple_GET_SIZE(__pyx_args) != 2) { goto __pyx_L5_argtuple_error; } else { values[0] = PyTuple_GET_ITEM(__pyx_args, 0); values[1] = PyTuple_GET_ITEM(__pyx_args, 1); } __pyx_v_x = ((PyArrayObject *)values[0]); __pyx_v_eps = __pyx_PyFloat_AsDouble(values[1]); if (unlikely((__pyx_v_eps == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } goto __pyx_L4_argument_unpacking_done; __pyx_L5_argtuple_error:; __Pyx_RaiseArgtupleInvalid("dxi_clust", 1, 2, 2, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_L3_error:; __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.dxi_clust", __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 = 51; __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __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 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_dxi_clust(__pyx_v_x, __pyx_v_eps, 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 51; __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_sage79_spyx_0.dxi_clust", __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; }
52: cdef double Lambda = x[0]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":52 * * cpdef dxi_clust(np.ndarray[cdouble_t] x, double eps): * cdef double Lambda = x[0] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double rho = x[1] * cdef double zeta = 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 = 52; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_Lambda = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_1, __pyx_bstride_0_x));
53: cdef double rho = x[1]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":53 * cpdef dxi_clust(np.ndarray[cdouble_t] x, double eps): * cdef double Lambda = x[0] * cdef double rho = x[1] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double zeta = x[2] * return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2, */ __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 = 53; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_rho = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_3, __pyx_bstride_0_x));
54: cdef double zeta = x[2]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":54 * cdef double Lambda = x[0] * cdef double rho = x[1] * cdef double zeta = x[2] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2, * 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, */ __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 = 54; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_v_zeta = (*__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_x.buf, __pyx_t_4, __pyx_bstride_0_x));
55: return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":55 * cdef double rho = x[1] * cdef double zeta = x[2] * return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, * 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2, */ __Pyx_XDECREF(__pyx_r); /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":58 * 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, * 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2, * 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* * cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x): */ __pyx_t_5 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_5)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_5); __pyx_t_6 = PyObject_GetAttr(__pyx_t_5, __pyx_n_s__array); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_5); __pyx_t_5 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":55 * cdef double rho = x[1] * cdef double zeta = x[2] * return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, * 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2, */ __pyx_t_5 = PyFloat_FromDouble(((4.0 * __pyx_v_Lambda) * pow(__pyx_v_rho, 3.0))); if (unlikely(!__pyx_t_5)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_5); __pyx_t_7 = PyFloat_FromDouble((((4.0 * __pyx_v_Lambda) * pow(__pyx_v_rho, 2.0)) * __pyx_v_zeta)); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __pyx_t_8 = PyFloat_FromDouble(((3.0 * __pyx_v_Lambda) * pow(__pyx_v_rho, 2.0))); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8);
56: 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":56 * cdef double zeta = x[2] * return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2, * 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2, * 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps */ __pyx_t_9 = PyFloat_FromDouble((((4.0 * __pyx_v_Lambda) * __pyx_v_rho) * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_9)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 56; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_9); __pyx_t_10 = PyFloat_FromDouble((((3.0 * __pyx_v_Lambda) * __pyx_v_rho) * __pyx_v_zeta)); if (unlikely(!__pyx_t_10)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 56; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_10); __pyx_t_11 = PyFloat_FromDouble(((4.0 * __pyx_v_Lambda) * pow(__pyx_v_zeta, 3.0))); if (unlikely(!__pyx_t_11)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 56; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_11);
57: 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2,
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":57 * return np.array([4*Lambda*rho**3, 4*Lambda*rho**2*zeta, 3*Lambda*rho**2, * 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, * 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2, #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps * */ __pyx_t_12 = PyFloat_FromDouble(((3.0 * __pyx_v_Lambda) * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_12)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 57; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_12); __pyx_t_13 = PyFloat_FromDouble((3.0 * pow(__pyx_v_rho, 3.0))); if (unlikely(!__pyx_t_13)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 57; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_13); __pyx_t_14 = PyFloat_FromDouble(((3.0 * pow(__pyx_v_rho, 2.0)) * __pyx_v_zeta)); if (unlikely(!__pyx_t_14)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 57; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_14); __pyx_t_15 = PyFloat_FromDouble((2.0 * pow(__pyx_v_rho, 2.0))); if (unlikely(!__pyx_t_15)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 57; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_15); __pyx_t_16 = PyFloat_FromDouble(((3.0 * __pyx_v_rho) * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_16)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 57; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_16);
58: 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":58 * 4*Lambda*rho*zeta**2, 3*Lambda*rho*zeta, 4*Lambda*zeta**3, * 3*Lambda*zeta**2, 3*rho**3, 3*rho**2*zeta, 2*rho**2, 3*rho*zeta**2, * 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* * cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x): */ __pyx_t_17 = PyFloat_FromDouble(((2.0 * __pyx_v_rho) * __pyx_v_zeta)); if (unlikely(!__pyx_t_17)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_17); __pyx_t_18 = PyFloat_FromDouble((3.0 * pow(__pyx_v_zeta, 3.0))); if (unlikely(!__pyx_t_18)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_18); __pyx_t_19 = PyFloat_FromDouble((2.0 * pow(__pyx_v_zeta, 2.0))); if (unlikely(!__pyx_t_19)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_19); __pyx_t_20 = PyFloat_FromDouble(__pyx_v_zeta); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_21 = PyList_New(15); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_21)); PyList_SET_ITEM(__pyx_t_21, 0, __pyx_t_5); __Pyx_GIVEREF(__pyx_t_5); PyList_SET_ITEM(__pyx_t_21, 1, __pyx_t_7); __Pyx_GIVEREF(__pyx_t_7); PyList_SET_ITEM(__pyx_t_21, 2, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); PyList_SET_ITEM(__pyx_t_21, 3, __pyx_t_9); __Pyx_GIVEREF(__pyx_t_9); PyList_SET_ITEM(__pyx_t_21, 4, __pyx_t_10); __Pyx_GIVEREF(__pyx_t_10); PyList_SET_ITEM(__pyx_t_21, 5, __pyx_t_11); __Pyx_GIVEREF(__pyx_t_11); PyList_SET_ITEM(__pyx_t_21, 6, __pyx_t_12); __Pyx_GIVEREF(__pyx_t_12); PyList_SET_ITEM(__pyx_t_21, 7, __pyx_t_13); __Pyx_GIVEREF(__pyx_t_13); PyList_SET_ITEM(__pyx_t_21, 8, __pyx_t_14); __Pyx_GIVEREF(__pyx_t_14); PyList_SET_ITEM(__pyx_t_21, 9, __pyx_t_15); __Pyx_GIVEREF(__pyx_t_15); PyList_SET_ITEM(__pyx_t_21, 10, __pyx_t_16); __Pyx_GIVEREF(__pyx_t_16); PyList_SET_ITEM(__pyx_t_21, 11, __pyx_t_17); __Pyx_GIVEREF(__pyx_t_17); PyList_SET_ITEM(__pyx_t_21, 12, __pyx_t_18); __Pyx_GIVEREF(__pyx_t_18); PyList_SET_ITEM(__pyx_t_21, 13, __pyx_t_19); __Pyx_GIVEREF(__pyx_t_19); PyList_SET_ITEM(__pyx_t_21, 14, __pyx_t_20); __Pyx_GIVEREF(__pyx_t_20); __pyx_t_5 = 0; __pyx_t_7 = 0; __pyx_t_8 = 0; __pyx_t_9 = 0; __pyx_t_10 = 0; __pyx_t_11 = 0; __pyx_t_12 = 0; __pyx_t_13 = 0; __pyx_t_14 = 0; __pyx_t_15 = 0; __pyx_t_16 = 0; __pyx_t_17 = 0; __pyx_t_18 = 0; __pyx_t_19 = 0; __pyx_t_20 = 0; __pyx_t_20 = PyTuple_New(1); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_20)); PyTuple_SET_ITEM(__pyx_t_20, 0, ((PyObject *)__pyx_t_21)); __Pyx_GIVEREF(((PyObject *)__pyx_t_21)); __pyx_t_21 = 0; __pyx_t_21 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_20), NULL); if (unlikely(!__pyx_t_21)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 55; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_21); __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_20)); __pyx_t_20 = 0; __pyx_t_20 = PyFloat_FromDouble(__pyx_v_eps); if (unlikely(!__pyx_t_20)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_20); __pyx_t_6 = PyNumber_Multiply(__pyx_t_21, __pyx_t_20); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 58; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_21); __pyx_t_21 = 0; __Pyx_DECREF(__pyx_t_20); __pyx_t_20 = 0; __pyx_r = __pyx_t_6; __pyx_t_6 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_5); __Pyx_XDECREF(__pyx_t_6); __Pyx_XDECREF(__pyx_t_7); __Pyx_XDECREF(__pyx_t_8); __Pyx_XDECREF(__pyx_t_9); __Pyx_XDECREF(__pyx_t_10); __Pyx_XDECREF(__pyx_t_11); __Pyx_XDECREF(__pyx_t_12); __Pyx_XDECREF(__pyx_t_13); __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); { 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_sage79_spyx_0.dxi_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = 0; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
59:
60: cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x):
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":60 * 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps * * cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef int m = x.shape[0] * cdef np.ndarray[cdouble_t] I = np.empty(m) */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_3inv_clust(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static PyObject *__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_inv_clust(PyArrayObject *__pyx_v_a, PyArrayObject *__pyx_v_x, int __pyx_skip_dispatch) { int __pyx_v_m; PyArrayObject *__pyx_v_I = 0; int __pyx_v_i; PyObject *__pyx_v_coeffs = NULL; PyObject *__pyx_v_monoms = NULL; Py_buffer __pyx_bstruct_a; Py_ssize_t __pyx_bstride_0_a = 0; Py_ssize_t __pyx_bshape_0_a = 0; Py_buffer __pyx_bstruct_I; Py_ssize_t __pyx_bstride_0_I = 0; Py_ssize_t __pyx_bshape_0_I = 0; Py_buffer __pyx_bstruct_x; Py_ssize_t __pyx_bstride_0_x = 0; Py_ssize_t __pyx_bstride_1_x = 0; Py_ssize_t __pyx_bshape_0_x = 0; Py_ssize_t __pyx_bshape_1_x = 0; PyObject *__pyx_r = NULL; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("inv_clust"); __pyx_bstruct_I.buf = NULL; __pyx_bstruct_a.buf = NULL; __pyx_bstruct_x.buf = NULL; { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_a, (PyObject*)__pyx_v_a, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_a = __pyx_bstruct_a.strides[0]; __pyx_bshape_0_a = __pyx_bstruct_a.shape[0]; { __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 2, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_x = __pyx_bstruct_x.strides[0]; __pyx_bstride_1_x = __pyx_bstruct_x.strides[1]; __pyx_bshape_0_x = __pyx_bstruct_x.shape[0]; __pyx_bshape_1_x = __pyx_bstruct_x.shape[1]; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":60 * 2*rho*zeta, 3*zeta**3, 2*zeta**2, zeta]) * eps * * cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef int m = x.shape[0] * cdef np.ndarray[cdouble_t] I = np.empty(m) */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_3inv_clust(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_sage79_spyx_0_3inv_clust[] = "File: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx (starting at line 60)"; static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_3inv_clust(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) { PyArrayObject *__pyx_v_a = 0; PyArrayObject *__pyx_v_x = 0; Py_buffer __pyx_bstruct_a; Py_ssize_t __pyx_bstride_0_a = 0; Py_ssize_t __pyx_bshape_0_a = 0; Py_buffer __pyx_bstruct_x; Py_ssize_t __pyx_bstride_0_x = 0; Py_ssize_t __pyx_bstride_1_x = 0; Py_ssize_t __pyx_bshape_0_x = 0; Py_ssize_t __pyx_bshape_1_x = 0; PyObject *__pyx_r = NULL; static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__a,&__pyx_n_s__x,0}; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("inv_clust"); __pyx_self = __pyx_self; { PyObject* values[2] = {0,0}; if (unlikely(__pyx_kwds)) { Py_ssize_t kw_args; switch (PyTuple_GET_SIZE(__pyx_args)) { 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__a); if (likely(values[0])) kw_args--; else goto __pyx_L5_argtuple_error; case 1: values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__x); if (likely(values[1])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("inv_clust", 1, 2, 2, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __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), "inv_clust")<
0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } } else if (PyTuple_GET_SIZE(__pyx_args) != 2) { goto __pyx_L5_argtuple_error; } else { values[0] = PyTuple_GET_ITEM(__pyx_args, 0); values[1] = PyTuple_GET_ITEM(__pyx_args, 1); } __pyx_v_a = ((PyArrayObject *)values[0]); __pyx_v_x = ((PyArrayObject *)values[1]); } goto __pyx_L4_argument_unpacking_done; __pyx_L5_argtuple_error:; __Pyx_RaiseArgtupleInvalid("inv_clust", 1, 2, 2, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_L3_error:; __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.inv_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __Pyx_RefNannyFinishContext(); return NULL; __pyx_L4_argument_unpacking_done:; __pyx_bstruct_a.buf = NULL; __pyx_bstruct_x.buf = NULL; if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_a), __pyx_ptype_5numpy_ndarray, 1, "a", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;} if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x), __pyx_ptype_5numpy_ndarray, 1, "x", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;} { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_a, (PyObject*)__pyx_v_a, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_a = __pyx_bstruct_a.strides[0]; __pyx_bshape_0_a = __pyx_bstruct_a.shape[0]; { __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 2, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_x = __pyx_bstruct_x.strides[0]; __pyx_bstride_1_x = __pyx_bstruct_x.strides[1]; __pyx_bshape_0_x = __pyx_bstruct_x.shape[0]; __pyx_bshape_1_x = __pyx_bstruct_x.shape[1]; __Pyx_XDECREF(__pyx_r); __pyx_t_1 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_inv_clust(__pyx_v_a, __pyx_v_x, 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 60; __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_a); __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_sage79_spyx_0.inv_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = NULL; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
61: cdef int m = x.shape[0]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":61 * * cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x): * cdef int m = x.shape[0] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef np.ndarray[cdouble_t] I = np.empty(m) * cdef int i */ __pyx_v_m = (__pyx_v_x->dimensions[0]);
62: cdef np.ndarray[cdouble_t] I = np.empty(m)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":62 * cpdef inv_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x): * cdef int m = x.shape[0] * cdef np.ndarray[cdouble_t] I = np.empty(m) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef int i * coeffs = coeffs_clust(a) */ __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_2 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__empty); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyInt_FromLong(__pyx_v_m); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_1); __Pyx_GIVEREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; if (!(likely(((__pyx_t_1) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_1, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_4 = ((PyArrayObject *)__pyx_t_1); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_I, (PyObject*)__pyx_t_4, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) { __pyx_v_I = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_I.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 62; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_I = __pyx_bstruct_I.strides[0]; __pyx_bshape_0_I = __pyx_bstruct_I.shape[0]; } } __pyx_t_4 = 0; __pyx_v_I = ((PyArrayObject *)__pyx_t_1); __pyx_t_1 = 0;
63: cdef int i
64: coeffs = coeffs_clust(a)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":64 * cdef np.ndarray[cdouble_t] I = np.empty(m) * cdef int i * coeffs = coeffs_clust(a) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* for i from 0<
= i<
m: * monoms = monoms_clust(x[i]) */ __pyx_t_1 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(((PyArrayObject *)__pyx_v_a), 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 64; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_v_coeffs = __pyx_t_1; __pyx_t_1 = 0;
65: for i from 0 <= i < m:
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":65 * cdef int i * coeffs = coeffs_clust(a) * for i from 0<
= i<
m: #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* monoms = monoms_clust(x[i]) * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(monoms)) */ __pyx_t_5 = __pyx_v_m; for (__pyx_v_i = 0; __pyx_v_i<
__pyx_t_5; __pyx_v_i++) {
66: monoms = monoms_clust(x[i])
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":66 * coeffs = coeffs_clust(a) * for i from 0<
= i<
m: * monoms = monoms_clust(x[i]) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(monoms)) * return np.linalg.norm(I) */ __pyx_t_1 = __Pyx_GetItemInt(((PyObject *)__pyx_v_x), __pyx_v_i, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 66; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (!(likely(((__pyx_t_1) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_1, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 66; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_3 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_monoms_clust(((PyArrayObject *)__pyx_t_1), 0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 66; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_XDECREF(__pyx_v_monoms); __pyx_v_monoms = __pyx_t_3; __pyx_t_3 = 0;
67: I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(monoms))
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":67 * for i from 0<
= i<
m: * monoms = monoms_clust(x[i]) * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(monoms)) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return np.linalg.norm(I) * */ __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__dot); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyTuple_New(2); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); __Pyx_INCREF(__pyx_v_coeffs); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_v_coeffs); __Pyx_GIVEREF(__pyx_v_coeffs); __Pyx_INCREF(__pyx_v_monoms); PyTuple_SET_ITEM(__pyx_t_3, 1, __pyx_v_monoms); __Pyx_GIVEREF(__pyx_v_monoms); __pyx_t_2 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_1 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__dot); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_6 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__abs); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); __Pyx_INCREF(__pyx_v_coeffs); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_v_coeffs); __Pyx_GIVEREF(__pyx_v_coeffs); __pyx_t_7 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_6 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__abs); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); __Pyx_INCREF(__pyx_v_monoms); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_v_monoms); __Pyx_GIVEREF(__pyx_v_monoms); __pyx_t_8 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; __pyx_t_3 = PyTuple_New(2); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_7); __Pyx_GIVEREF(__pyx_t_7); PyTuple_SET_ITEM(__pyx_t_3, 1, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); __pyx_t_7 = 0; __pyx_t_8 = 0; __pyx_t_8 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; __pyx_t_3 = __Pyx_PyNumber_Divide(__pyx_t_2, __pyx_t_8); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_9 = __pyx_PyFloat_AsDouble(__pyx_t_3); if (unlikely((__pyx_t_9 == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_10 = __pyx_v_i; __pyx_t_11 = -1; if (__pyx_t_10<
0) { __pyx_t_10 += __pyx_bshape_0_I; if (unlikely(__pyx_t_10<
0)) __pyx_t_11 = 0; } else if (unlikely(__pyx_t_10 >= __pyx_bshape_0_I)) __pyx_t_11 = 0; if (unlikely(__pyx_t_11 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_11); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 67; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } *__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_I.buf, __pyx_t_10, __pyx_bstride_0_I) = __pyx_t_9; }
68: return np.linalg.norm(I)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":68 * monoms = monoms_clust(x[i]) * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(monoms)) * return np.linalg.norm(I) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* * cpdef logL_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x, double eps): */ __Pyx_XDECREF(__pyx_r); __pyx_t_3 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __pyx_t_8 = PyObject_GetAttr(__pyx_t_3, __pyx_n_s__linalg); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __pyx_t_3 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__norm); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyTuple_New(1); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_8)); __Pyx_INCREF(((PyObject *)__pyx_v_I)); PyTuple_SET_ITEM(__pyx_t_8, 0, ((PyObject *)__pyx_v_I)); __Pyx_GIVEREF(((PyObject *)__pyx_v_I)); __pyx_t_2 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_8), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 68; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0; __pyx_r = __pyx_t_2; __pyx_t_2 = 0; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_1); __Pyx_XDECREF(__pyx_t_2); __Pyx_XDECREF(__pyx_t_3); __Pyx_XDECREF(__pyx_t_6); __Pyx_XDECREF(__pyx_t_7); __Pyx_XDECREF(__pyx_t_8); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_I); __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_sage79_spyx_0.inv_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = 0; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_I); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XDECREF((PyObject *)__pyx_v_I); __Pyx_XDECREF(__pyx_v_coeffs); __Pyx_XDECREF(__pyx_v_monoms); __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
69:
70: cpdef logL_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x, double eps):
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":70 * return np.linalg.norm(I) * * cpdef logL_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x, double eps): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef int m = x.shape[0] * cdef np.ndarray[cdouble_t] I = np.empty(m) */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_4logL_clust(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds); /*proto*/ static PyObject *__pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_logL_clust(PyArrayObject *__pyx_v_a, PyArrayObject *__pyx_v_x, double __pyx_v_eps, int __pyx_skip_dispatch) { int __pyx_v_m; PyArrayObject *__pyx_v_I = 0; int __pyx_v_i; PyObject *__pyx_v_coeffs = NULL; PyObject *__pyx_v_monoms = NULL; PyObject *__pyx_v_dxi = NULL; double __pyx_v_nrm; PyObject *__pyx_v_logL = NULL; Py_buffer __pyx_bstruct_a; Py_ssize_t __pyx_bstride_0_a = 0; Py_ssize_t __pyx_bshape_0_a = 0; Py_buffer __pyx_bstruct_I; Py_ssize_t __pyx_bstride_0_I = 0; Py_ssize_t __pyx_bshape_0_I = 0; Py_buffer __pyx_bstruct_x; Py_ssize_t __pyx_bstride_0_x = 0; Py_ssize_t __pyx_bstride_1_x = 0; Py_ssize_t __pyx_bshape_0_x = 0; Py_ssize_t __pyx_bshape_1_x = 0; PyObject *__pyx_r = NULL; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("logL_clust"); __pyx_bstruct_I.buf = NULL; __pyx_bstruct_a.buf = NULL; __pyx_bstruct_x.buf = NULL; { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_a, (PyObject*)__pyx_v_a, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_a = __pyx_bstruct_a.strides[0]; __pyx_bshape_0_a = __pyx_bstruct_a.shape[0]; { __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 2, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_x = __pyx_bstruct_x.strides[0]; __pyx_bstride_1_x = __pyx_bstruct_x.strides[1]; __pyx_bshape_0_x = __pyx_bstruct_x.shape[0]; __pyx_bshape_1_x = __pyx_bstruct_x.shape[1]; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":70 * return np.linalg.norm(I) * * cpdef logL_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x, double eps): #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef int m = x.shape[0] * cdef np.ndarray[cdouble_t] I = np.empty(m) */ static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_4logL_clust(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_sage79_spyx_0_4logL_clust[] = "File: _sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx (starting at line 70)"; static PyObject *__pyx_pf_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_4logL_clust(PyObject *__pyx_self, PyObject *__pyx_args, PyObject *__pyx_kwds) { PyArrayObject *__pyx_v_a = 0; PyArrayObject *__pyx_v_x = 0; double __pyx_v_eps; Py_buffer __pyx_bstruct_a; Py_ssize_t __pyx_bstride_0_a = 0; Py_ssize_t __pyx_bshape_0_a = 0; Py_buffer __pyx_bstruct_x; Py_ssize_t __pyx_bstride_0_x = 0; Py_ssize_t __pyx_bstride_1_x = 0; Py_ssize_t __pyx_bshape_0_x = 0; Py_ssize_t __pyx_bshape_1_x = 0; PyObject *__pyx_r = NULL; static PyObject **__pyx_pyargnames[] = {&__pyx_n_s__a,&__pyx_n_s__x,&__pyx_n_s__eps,0}; __Pyx_RefNannyDeclarations __Pyx_RefNannySetupContext("logL_clust"); __pyx_self = __pyx_self; { PyObject* values[3] = {0,0,0}; if (unlikely(__pyx_kwds)) { Py_ssize_t kw_args; switch (PyTuple_GET_SIZE(__pyx_args)) { 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__a); if (likely(values[0])) kw_args--; else goto __pyx_L5_argtuple_error; case 1: values[1] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__x); if (likely(values[1])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("logL_clust", 1, 3, 3, 1); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } case 2: values[2] = PyDict_GetItem(__pyx_kwds, __pyx_n_s__eps); if (likely(values[2])) kw_args--; else { __Pyx_RaiseArgtupleInvalid("logL_clust", 1, 3, 3, 2); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __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), "logL_clust")<
0)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } } else if (PyTuple_GET_SIZE(__pyx_args) != 3) { 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); } __pyx_v_a = ((PyArrayObject *)values[0]); __pyx_v_x = ((PyArrayObject *)values[1]); __pyx_v_eps = __pyx_PyFloat_AsDouble(values[2]); if (unlikely((__pyx_v_eps == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L3_error;} } goto __pyx_L4_argument_unpacking_done; __pyx_L5_argtuple_error:; __Pyx_RaiseArgtupleInvalid("logL_clust", 1, 3, 3, PyTuple_GET_SIZE(__pyx_args)); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L3_error;} __pyx_L3_error:; __Pyx_AddTraceback("_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.logL_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __Pyx_RefNannyFinishContext(); return NULL; __pyx_L4_argument_unpacking_done:; __pyx_bstruct_a.buf = NULL; __pyx_bstruct_x.buf = NULL; if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_a), __pyx_ptype_5numpy_ndarray, 1, "a", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;} if (unlikely(!__Pyx_ArgTypeTest(((PyObject *)__pyx_v_x), __pyx_ptype_5numpy_ndarray, 1, "x", 0))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;} { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_a, (PyObject*)__pyx_v_a, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 1, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_a = __pyx_bstruct_a.strides[0]; __pyx_bshape_0_a = __pyx_bstruct_a.shape[0]; { __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_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES, 2, 0, __pyx_stack) == -1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } __pyx_bstride_0_x = __pyx_bstruct_x.strides[0]; __pyx_bstride_1_x = __pyx_bstruct_x.strides[1]; __pyx_bshape_0_x = __pyx_bstruct_x.shape[0]; __pyx_bshape_1_x = __pyx_bstruct_x.shape[1]; __Pyx_XDECREF(__pyx_r); __pyx_t_1 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_logL_clust(__pyx_v_a, __pyx_v_x, __pyx_v_eps, 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 70; __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_a); __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_sage79_spyx_0.logL_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = NULL; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }
71: cdef int m = x.shape[0]
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":71 * * cpdef logL_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x, double eps): * cdef int m = x.shape[0] #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef np.ndarray[cdouble_t] I = np.empty(m) * cdef int i */ __pyx_v_m = (__pyx_v_x->dimensions[0]);
72: cdef np.ndarray[cdouble_t] I = np.empty(m)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":72 * cpdef logL_clust(np.ndarray[cdouble_t] a, np.ndarray[cdouble_t, ndim=2] x, double eps): * cdef int m = x.shape[0] * cdef np.ndarray[cdouble_t] I = np.empty(m) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef int i * coeffs = coeffs_clust(a) */ __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_2 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__empty); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyInt_FromLong(__pyx_v_m); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_1); __Pyx_GIVEREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; if (!(likely(((__pyx_t_1) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_1, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_4 = ((PyArrayObject *)__pyx_t_1); { __Pyx_BufFmt_StackElem __pyx_stack[1]; if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_bstruct_I, (PyObject*)__pyx_t_4, &__Pyx_TypeInfo_nn___pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 1, 0, __pyx_stack) == -1)) { __pyx_v_I = ((PyArrayObject *)Py_None); __Pyx_INCREF(Py_None); __pyx_bstruct_I.buf = NULL; {__pyx_filename = __pyx_f[0]; __pyx_lineno = 72; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } else {__pyx_bstride_0_I = __pyx_bstruct_I.strides[0]; __pyx_bshape_0_I = __pyx_bstruct_I.shape[0]; } } __pyx_t_4 = 0; __pyx_v_I = ((PyArrayObject *)__pyx_t_1); __pyx_t_1 = 0;
73: cdef int i
74: coeffs = coeffs_clust(a)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":74 * cdef np.ndarray[cdouble_t] I = np.empty(m) * cdef int i * coeffs = coeffs_clust(a) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* for i from 0<
= i<
m: * monoms = monoms_clust(x[i]) */ __pyx_t_1 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_coeffs_clust(((PyArrayObject *)__pyx_v_a), 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 74; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_v_coeffs = __pyx_t_1; __pyx_t_1 = 0;
75: for i from 0 <= i < m:
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":75 * cdef int i * coeffs = coeffs_clust(a) * for i from 0<
= i<
m: #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* monoms = monoms_clust(x[i]) * dxi = dxi_clust(x[i], eps) */ __pyx_t_5 = __pyx_v_m; for (__pyx_v_i = 0; __pyx_v_i<
__pyx_t_5; __pyx_v_i++) {
76: monoms = monoms_clust(x[i])
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":76 * coeffs = coeffs_clust(a) * for i from 0<
= i<
m: * monoms = monoms_clust(x[i]) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* dxi = dxi_clust(x[i], eps) * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi)) */ __pyx_t_1 = __Pyx_GetItemInt(((PyObject *)__pyx_v_x), __pyx_v_i, sizeof(int), PyInt_FromLong); if (!__pyx_t_1) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); if (!(likely(((__pyx_t_1) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_1, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_3 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_monoms_clust(((PyArrayObject *)__pyx_t_1), 0); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 76; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_XDECREF(__pyx_v_monoms); __pyx_v_monoms = __pyx_t_3; __pyx_t_3 = 0;
77: dxi = dxi_clust(x[i], eps)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":77 * for i from 0<
= i<
m: * monoms = monoms_clust(x[i]) * dxi = dxi_clust(x[i], eps) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi)) * cdef double nrm = np.linalg.norm(I) */ __pyx_t_3 = __Pyx_GetItemInt(((PyObject *)__pyx_v_x), __pyx_v_i, sizeof(int), PyInt_FromLong); if (!__pyx_t_3) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 77; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); if (!(likely(((__pyx_t_3) == Py_None) || likely(__Pyx_TypeTest(__pyx_t_3, __pyx_ptype_5numpy_ndarray))))) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 77; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __pyx_t_1 = __pyx_f_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_dxi_clust(((PyArrayObject *)__pyx_t_3), __pyx_v_eps, 0); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 77; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_XDECREF(__pyx_v_dxi); __pyx_v_dxi = __pyx_t_1; __pyx_t_1 = 0;
78: I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi))
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":78 * monoms = monoms_clust(x[i]) * dxi = dxi_clust(x[i], eps) * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi)) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* cdef double nrm = np.linalg.norm(I) * logL = (1 - 0.5*m)*np.log(2) + (m - 1)*np.log(nrm) - 0.5*nrm**2 - np.log(gamma(0.5*m)) */ __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_3 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__dot); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyTuple_New(2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); __Pyx_INCREF(__pyx_v_coeffs); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_v_coeffs); __Pyx_GIVEREF(__pyx_v_coeffs); __Pyx_INCREF(__pyx_v_monoms); PyTuple_SET_ITEM(__pyx_t_1, 1, __pyx_v_monoms); __Pyx_GIVEREF(__pyx_v_monoms); __pyx_t_2 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_3 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__dot); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_6 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__abs); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); __Pyx_INCREF(__pyx_v_coeffs); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_v_coeffs); __Pyx_GIVEREF(__pyx_v_coeffs); __pyx_t_7 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_6 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__abs); if (unlikely(!__pyx_t_6)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_6); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyTuple_New(1); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); __Pyx_INCREF(__pyx_v_dxi); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_v_dxi); __Pyx_GIVEREF(__pyx_v_dxi); __pyx_t_8 = PyObject_Call(__pyx_t_6, ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_6); __pyx_t_6 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __pyx_t_1 = PyTuple_New(2); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_1)); PyTuple_SET_ITEM(__pyx_t_1, 0, __pyx_t_7); __Pyx_GIVEREF(__pyx_t_7); PyTuple_SET_ITEM(__pyx_t_1, 1, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); __pyx_t_7 = 0; __pyx_t_8 = 0; __pyx_t_8 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_1), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_1)); __pyx_t_1 = 0; __pyx_t_1 = __Pyx_PyNumber_Divide(__pyx_t_2, __pyx_t_8); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_9 = __pyx_PyFloat_AsDouble(__pyx_t_1); if (unlikely((__pyx_t_9 == (npy_double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_10 = __pyx_v_i; __pyx_t_11 = -1; if (__pyx_t_10<
0) { __pyx_t_10 += __pyx_bshape_0_I; if (unlikely(__pyx_t_10<
0)) __pyx_t_11 = 0; } else if (unlikely(__pyx_t_10 >= __pyx_bshape_0_I)) __pyx_t_11 = 0; if (unlikely(__pyx_t_11 != -1)) { __Pyx_RaiseBufferIndexError(__pyx_t_11); {__pyx_filename = __pyx_f[0]; __pyx_lineno = 78; __pyx_clineno = __LINE__; goto __pyx_L1_error;} } *__Pyx_BufPtrStrided1d(__pyx_t_75_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0_cdouble_t *, __pyx_bstruct_I.buf, __pyx_t_10, __pyx_bstride_0_I) = __pyx_t_9; }
79: cdef double nrm = np.linalg.norm(I)
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":79 * dxi = dxi_clust(x[i], eps) * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi)) * cdef double nrm = np.linalg.norm(I) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* logL = (1 - 0.5*m)*np.log(2) + (m - 1)*np.log(nrm) - 0.5*nrm**2 - np.log(gamma(0.5*m)) * return logL */ __pyx_t_1 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __pyx_t_8 = PyObject_GetAttr(__pyx_t_1, __pyx_n_s__linalg); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__norm); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyTuple_New(1); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_8)); __Pyx_INCREF(((PyObject *)__pyx_v_I)); PyTuple_SET_ITEM(__pyx_t_8, 0, ((PyObject *)__pyx_v_I)); __Pyx_GIVEREF(((PyObject *)__pyx_v_I)); __pyx_t_2 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_t_8), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_8)); __pyx_t_8 = 0; __pyx_t_12 = __pyx_PyFloat_AsDouble(__pyx_t_2); if (unlikely((__pyx_t_12 == (double)-1) && PyErr_Occurred())) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 79; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_v_nrm = __pyx_t_12;
80: logL = (1 - 0.5*m)*np.log(2) + (m - 1)*np.log(nrm) - 0.5*nrm**2 - np.log(gamma(0.5*m))
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":80 * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi)) * cdef double nrm = np.linalg.norm(I) * logL = (1 - 0.5*m)*np.log(2) + (m - 1)*np.log(nrm) - 0.5*nrm**2 - np.log(gamma(0.5*m)) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return logL */ __pyx_t_2 = PyFloat_FromDouble((1.0 - (0.5 * __pyx_v_m))); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_8 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_1 = PyObject_GetAttr(__pyx_t_8, __pyx_n_s__log); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyObject_Call(__pyx_t_1, ((PyObject *)__pyx_k_tuple_1), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __pyx_t_1 = PyNumber_Multiply(__pyx_t_2, __pyx_t_8); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyInt_FromLong((__pyx_v_m - 1)); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_2 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_3 = PyObject_GetAttr(__pyx_t_2, __pyx_n_s__log); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_2 = PyFloat_FromDouble(__pyx_v_nrm); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __pyx_t_7 = PyTuple_New(1); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_7)); PyTuple_SET_ITEM(__pyx_t_7, 0, __pyx_t_2); __Pyx_GIVEREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_2 = PyObject_Call(__pyx_t_3, ((PyObject *)__pyx_t_7), NULL); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_3); __pyx_t_3 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_7)); __pyx_t_7 = 0; __pyx_t_7 = PyNumber_Multiply(__pyx_t_8, __pyx_t_2); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __pyx_t_2 = PyNumber_Add(__pyx_t_1, __pyx_t_7); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0; __pyx_t_7 = PyFloat_FromDouble((0.5 * pow(__pyx_v_nrm, 2.0))); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __pyx_t_1 = PyNumber_Subtract(__pyx_t_2, __pyx_t_7); if (unlikely(!__pyx_t_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_1); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0; __pyx_t_7 = __Pyx_GetName(__pyx_m, __pyx_n_s__np); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __pyx_t_2 = PyObject_GetAttr(__pyx_t_7, __pyx_n_s__log); if (unlikely(!__pyx_t_2)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_2); __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0; __pyx_t_7 = __Pyx_GetName(__pyx_m, __pyx_n_s__gamma); if (unlikely(!__pyx_t_7)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_7); __pyx_t_8 = PyFloat_FromDouble((0.5 * __pyx_v_m)); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyObject_Call(__pyx_t_7, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_7); __pyx_t_7 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; __pyx_t_3 = PyTuple_New(1); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_t_3)); PyTuple_SET_ITEM(__pyx_t_3, 0, __pyx_t_8); __Pyx_GIVEREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_t_8 = PyObject_Call(__pyx_t_2, ((PyObject *)__pyx_t_3), NULL); if (unlikely(!__pyx_t_8)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_8); __Pyx_DECREF(__pyx_t_2); __pyx_t_2 = 0; __Pyx_DECREF(((PyObject *)__pyx_t_3)); __pyx_t_3 = 0; __pyx_t_3 = PyNumber_Subtract(__pyx_t_1, __pyx_t_8); if (unlikely(!__pyx_t_3)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(__pyx_t_3); __Pyx_DECREF(__pyx_t_1); __pyx_t_1 = 0; __Pyx_DECREF(__pyx_t_8); __pyx_t_8 = 0; __pyx_v_logL = __pyx_t_3; __pyx_t_3 = 0; /* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":80 * I[i] = np.dot(coeffs, monoms) / np.dot(np.abs(coeffs), np.abs(dxi)) * cdef double nrm = np.linalg.norm(I) * logL = (1 - 0.5*m)*np.log(2) + (m - 1)*np.log(nrm) - 0.5*nrm**2 - np.log(gamma(0.5*m)) #<
<
<
<
<
<
<
<
<
<
<
<
<
<
* return logL */ __pyx_k_tuple_1 = PyTuple_New(1); if (unlikely(!__pyx_k_tuple_1)) {__pyx_filename = __pyx_f[0]; __pyx_lineno = 80; __pyx_clineno = __LINE__; goto __pyx_L1_error;} __Pyx_GOTREF(((PyObject *)__pyx_k_tuple_1)); __Pyx_INCREF(__pyx_int_2); PyTuple_SET_ITEM(__pyx_k_tuple_1, 0, __pyx_int_2); __Pyx_GIVEREF(__pyx_int_2); __Pyx_GIVEREF(((PyObject *)__pyx_k_tuple_1));
81: return logL
/* "_sagenb_servers_sage_notebook_sagenb_sagenb_home_klho_20_code_sage79_spyx_0.pyx":81 * cdef double nrm = np.linalg.norm(I) * logL = (1 - 0.5*m)*np.log(2) + (m - 1)*np.log(nrm) - 0.5*nrm**2 - np.log(gamma(0.5*m)) * return logL #<
<
<
<
<
<
<
<
<
<
<
<
<
<
*/ __Pyx_XDECREF(__pyx_r); __Pyx_INCREF(__pyx_v_logL); __pyx_r = __pyx_v_logL; goto __pyx_L0; __pyx_r = Py_None; __Pyx_INCREF(Py_None); goto __pyx_L0; __pyx_L1_error:; __Pyx_XDECREF(__pyx_t_1); __Pyx_XDECREF(__pyx_t_2); __Pyx_XDECREF(__pyx_t_3); __Pyx_XDECREF(__pyx_t_6); __Pyx_XDECREF(__pyx_t_7); __Pyx_XDECREF(__pyx_t_8); { PyObject *__pyx_type, *__pyx_value, *__pyx_tb; __Pyx_ErrFetch(&__pyx_type, &__pyx_value, &__pyx_tb); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_I); __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_sage79_spyx_0.logL_clust", __pyx_clineno, __pyx_lineno, __pyx_filename); __pyx_r = 0; goto __pyx_L2; __pyx_L0:; __Pyx_SafeReleaseBuffer(&__pyx_bstruct_a); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_I); __Pyx_SafeReleaseBuffer(&__pyx_bstruct_x); __pyx_L2:; __Pyx_XDECREF((PyObject *)__pyx_v_I); __Pyx_XDECREF(__pyx_v_coeffs); __Pyx_XDECREF(__pyx_v_monoms); __Pyx_XDECREF(__pyx_v_dxi); __Pyx_XDECREF(__pyx_v_logL); __Pyx_XGIVEREF(__pyx_r); __Pyx_RefNannyFinishContext(); return __pyx_r; }