Skip to content

Instantly share code, notes, and snippets.

/heat.c

Created January 22, 2017 18:06
Show Gist options
  • Save anonymous/e54964121b212e5f783fb10c696ed9e2 to your computer and use it in GitHub Desktop.
Save anonymous/e54964121b212e5f783fb10c696ed9e2 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
heat.c
__pyx_t_9 = __pyx_v_iter;
for (__pyx_t_10 = 0; __pyx_t_10 < __pyx_t_9; __pyx_t_10+=1) {
__pyx_v_count = __pyx_t_10;
__pyx_t_11 = (__pyx_v_M - 1);
for (__pyx_t_12 = 1; __pyx_t_12 < __pyx_t_11; __pyx_t_12+=1) {
__pyx_v_i = __pyx_t_12;
__pyx_t_13 = (__pyx_v_N - 1);
for (__pyx_t_14 = 1; __pyx_t_14 < __pyx_t_13; __pyx_t_14+=1) {
__pyx_v_j = __pyx_t_14;
__pyx_t_15 = (__pyx_v_i - 1);
__pyx_t_16 = __pyx_v_j;
if (__pyx_t_15 < 0) __pyx_t_15 += __pyx_pybuffernd_cur.diminfo[0].shape;
if (__pyx_t_16 < 0) __pyx_t_16 += __pyx_pybuffernd_cur.diminfo[1].shape;
__pyx_t_17 = (__pyx_v_i + 1);
__pyx_t_18 = __pyx_v_j;
if (__pyx_t_17 < 0) __pyx_t_17 += __pyx_pybuffernd_cur.diminfo[0].shape;
if (__pyx_t_18 < 0) __pyx_t_18 += __pyx_pybuffernd_cur.diminfo[1].shape;
__pyx_t_19 = __pyx_v_i;
__pyx_t_20 = (__pyx_v_j - 1);
if (__pyx_t_19 < 0) __pyx_t_19 += __pyx_pybuffernd_cur.diminfo[0].shape;
if (__pyx_t_20 < 0) __pyx_t_20 += __pyx_pybuffernd_cur.diminfo[1].shape;
__pyx_t_21 = __pyx_v_i;
__pyx_t_22 = (__pyx_v_j + 1);
if (__pyx_t_21 < 0) __pyx_t_21 += __pyx_pybuffernd_cur.diminfo[0].shape;
if (__pyx_t_22 < 0) __pyx_t_22 += __pyx_pybuffernd_cur.diminfo[1].shape;
__pyx_t_23 = __pyx_v_i;
__pyx_t_24 = __pyx_v_j;
if (__pyx_t_23 < 0) __pyx_t_23 += __pyx_pybuffernd_cur.diminfo[0].shape;
if (__pyx_t_24 < 0) __pyx_t_24 += __pyx_pybuffernd_cur.diminfo[1].shape;
__pyx_v_step = (((((*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_cur.rcbuffer->pybuffer.buf, __pyx_t_15, __pyx_pybuffernd_cur.diminfo[0].strides, __pyx_t_16, __pyx_pybuffernd_cur.diminfo[1].strides)) + (*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_cur.rcbuffer->pybuffer.buf, __pyx_t_17, __pyx_pybuffernd_cur.diminfo[0].strides, __pyx_t_18, __pyx_pybuffernd_cur.diminfo[1].strides))) + (*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_cur.rcbuffer->pybuffer.buf, __pyx_t_19, __pyx_pybuffernd_cur.diminfo[0].strides, __pyx_t_20, __pyx_pybuffernd_cur.diminfo[1].strides))) + (*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_cur.rcbuffer->pybuffer.buf, __pyx_t_21, __pyx_pybuffernd_cur.diminfo[0].strides, __pyx_t_22, __pyx_pybuffernd_cur.diminfo[1].strides))) - (4.0 * (*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_cur.rcbuffer->pybuffer.buf, __pyx_t_23, __pyx_pybuffernd_cur.diminfo[0].strides, __pyx_t_24, __pyx_pybuffernd_cur.diminfo[1].strides))));
__pyx_t_25 = __pyx_v_i;
__pyx_t_26 = __pyx_v_j;
if (__pyx_t_25 < 0) __pyx_t_25 += __pyx_pybuffernd_cur.diminfo[0].shape;
if (__pyx_t_26 < 0) __pyx_t_26 += __pyx_pybuffernd_cur.diminfo[1].shape;
__pyx_t_27 = (__pyx_v_dt * __pyx_v_step);
__pyx_t_28 = pow(__pyx_v_dx, 2.0);
if (unlikely(__pyx_t_28 == 0)) {
PyErr_SetString(PyExc_ZeroDivisionError, "float division");
__PYX_ERR(0, 17, __pyx_L1_error)
}
__pyx_t_29 = __pyx_v_i;
__pyx_t_30 = __pyx_v_j;
if (__pyx_t_29 < 0) __pyx_t_29 += __pyx_pybuffernd_next.diminfo[0].shape;
if (__pyx_t_30 < 0) __pyx_t_30 += __pyx_pybuffernd_next.diminfo[1].shape;
*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_next.rcbuffer->pybuffer.buf, __pyx_t_29, __pyx_pybuffernd_next.diminfo[0].strides, __pyx_t_30, __pyx_pybuffernd_next.diminfo[1].strides) = ((*__Pyx_BufPtrStrided2d(double *, __pyx_pybuffernd_cur.rcbuffer->pybuffer.buf, __pyx_t_25, __pyx_pybuffernd_cur.diminfo[0].strides, __pyx_t_26, __pyx_pybuffernd_cur.diminfo[1].strides)) + (__pyx_t_27 / __pyx_t_28));
}
}
__pyx_t_31 = ((PyObject *)__pyx_v_next);
__pyx_t_32 = ((PyObject *)__pyx_v_cur);
{
__Pyx_BufFmt_StackElem __pyx_stack[1];
__Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_cur.rcbuffer->pybuffer);
__pyx_t_12 = __Pyx_GetBufferAndValidate(&__pyx_pybuffernd_cur.rcbuffer->pybuffer, (PyObject*)((PyArrayObject *)__pyx_t_31), &__Pyx_TypeInfo_double, PyBUF_FORMAT| PyBUF_STRIDES, 2, 0, __pyx_stack);
if (unlikely(__pyx_t_12 < 0)) {
PyErr_Fetch(&__pyx_t_33, &__pyx_t_34, &__pyx_t_35);
if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_cur.rcbuffer->pybuffer, (PyObject*)__pyx_v_cur, &__Pyx_TypeInfo_double, PyBUF_FORMAT| PyBUF_STRIDES, 2, 0, __pyx_stack) == -1)) {
Py_XDECREF(__pyx_t_33); Py_XDECREF(__pyx_t_34); Py_XDECREF(__pyx_t_35);
__Pyx_RaiseBufferFallbackError();
} else {
PyErr_Restore(__pyx_t_33, __pyx_t_34, __pyx_t_35);
}
}
__pyx_pybuffernd_cur.diminfo[0].strides = __pyx_pybuffernd_cur.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_cur.diminfo[0].shape = __pyx_pybuffernd_cur.rcbuffer->pybuffer.shape[0]; __pyx_pybuffernd_cur.diminfo[1].strides = __pyx_pybuffernd_cur.rcbuffer->pybuffer.strides[1]; __pyx_pybuffernd_cur.diminfo[1].shape = __pyx_pybuffernd_cur.rcbuffer->pybuffer.shape[1];
if (unlikely(__pyx_t_12 < 0)) __PYX_ERR(0, 18, __pyx_L1_error)
}
__pyx_v_cur = ((PyArrayObject *)__pyx_t_31);
__pyx_t_31 = 0;
{
__Pyx_BufFmt_StackElem __pyx_stack[1];
__Pyx_SafeReleaseBuffer(&__pyx_pybuffernd_next.rcbuffer->pybuffer);
__pyx_t_12 = __Pyx_GetBufferAndValidate(&__pyx_pybuffernd_next.rcbuffer->pybuffer, (PyObject*)((PyArrayObject *)__pyx_t_32), &__Pyx_TypeInfo_double, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack);
if (unlikely(__pyx_t_12 < 0)) {
PyErr_Fetch(&__pyx_t_31, &__pyx_t_35, &__pyx_t_34);
if (unlikely(__Pyx_GetBufferAndValidate(&__pyx_pybuffernd_next.rcbuffer->pybuffer, (PyObject*)__pyx_v_next, &__Pyx_TypeInfo_double, PyBUF_FORMAT| PyBUF_STRIDES| PyBUF_WRITABLE, 2, 0, __pyx_stack) == -1)) {
Py_XDECREF(__pyx_t_31); Py_XDECREF(__pyx_t_35); Py_XDECREF(__pyx_t_34);
__Pyx_RaiseBufferFallbackError();
} else {
PyErr_Restore(__pyx_t_31, __pyx_t_35, __pyx_t_34);
}
}
__pyx_pybuffernd_next.diminfo[0].strides = __pyx_pybuffernd_next.rcbuffer->pybuffer.strides[0]; __pyx_pybuffernd_next.diminfo[0].shape = __pyx_pybuffernd_next.rcbuffer->pybuffer.shape[0]; __pyx_pybuffernd_next.diminfo[1].strides = __pyx_pybuffernd_next.rcbuffer->pybuffer.strides[1]; __pyx_pybuffernd_next.diminfo[1].shape = __pyx_pybuffernd_next.rcbuffer->pybuffer.shape[1];
if (unlikely(__pyx_t_12 < 0)) __PYX_ERR(0, 18, __pyx_L1_error)
}
__pyx_v_next = ((PyArrayObject *)__pyx_t_32);
__pyx_t_32 = 0;
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment