Need to debug taylor_recur start order

Author: hirish99

test/biharmonic.ipynb

@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 1,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -23,7 +23,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 10,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -32,13 +32,14 @@
     "abs_dist = sp.sqrt((var[0]-var_t[0])**2 + (var[1]-var_t[1])**2)\n",
     "w = make_identity_diff_op(2)\n",
     "\n",
-    "partial_1x = DerivativeIdentifier((4,0), 0)\n",
-    "partial_1y = DerivativeIdentifier((0,4), 0)\n",
-    "biharmonic_op = {partial_1x: 1, partial_1y: 1}\n",
+    "partial_4x = DerivativeIdentifier((4,0), 0)\n",
+    "partial_4y = DerivativeIdentifier((0,4), 0)\n",
+    "partial_2x2y = DerivativeIdentifier((2,2), 0)\n",
+    "biharmonic_op = {partial_4x: 1, partial_4y: 1, partial_2x2y:2}\n",
     "list_pde = immutabledict(biharmonic_op)\n",
     "\n",
     "biharmonic_pde = LinearPDESystemOperator(2, (list_pde,))\n",
-    "g_x_y = abs_dist**2 * (sp.log(abs_dist)-1)\n",
+    "g_x_y = abs_dist**2 * (sp.log(abs_dist))\n",
     "\n",
     "n_init, _, r = get_reindexed_and_center_origin_on_axis_recurrence(biharmonic_pde)\n",
     "\n",
@@ -55,7 +56,7 @@
     "s = sp.Function(\"s\")\n",
     "\n",
     "# pylint: disable-next=not-callable\n",
-    "subs_dict = {s(0): derivs[0], s(1): derivs[1], s(2): derivs[1], s(3): derivs[1]}\n",
+    "subs_dict = {s(0): derivs[0], s(1): derivs[1], s(2): derivs[2], s(3): derivs[3]}\n",
     "check = []\n",
     "\n",
     "assert n_init == 4\n",
@@ -68,7 +69,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": 11,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -81,52 +82,48 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 14,
+   "execution_count": 18,
    "metadata": {},
    "outputs": [
     {
      "data": {
       "text/latex": [
-       "$\\displaystyle 16.9964566798618$"
+       "$\\displaystyle -0.00388859878925869$"
       ],
       "text/plain": [
-       "16.9964566798618"
+       "-0.00388859878925869"
       ]
      },
-     "execution_count": 14,
+     "execution_count": 18,
      "metadata": {},
      "output_type": "execute_result"
     }
    ],
    "source": [
-    "r.subs(n, 4).subs(subs_dict).subs({var[0]: 1.2, var[1]: 2.3})"
+    "r.subs(n, 4).subs(subs_dict).subs({var[0]: 23, var[1]: 2.3})"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": 19,
    "metadata": {},
    "outputs": [
     {
      "data": {
+      "text/latex": [
+       "$\\displaystyle - \\frac{\\left(- 2 x_{0}^{2} + 6 x_{1}^{2}\\right) s{\\left(3 \\right)}}{x_{0}^{3} + x_{0} x_{1}^{2}} - \\frac{\\left(- x_{0}^{4} + 6 x_{0}^{2} x_{1}^{2} + 15 x_{1}^{4}\\right) s{\\left(1 \\right)}}{x_{0}^{7} + 2 x_{0}^{5} x_{1}^{2} + x_{0}^{3} x_{1}^{4}} + \\frac{\\left(x_{0}^{4} - 6 x_{0}^{2} x_{1}^{2} - 15 x_{1}^{4}\\right) s{\\left(2 \\right)}}{x_{0}^{6} + 2 x_{0}^{4} x_{1}^{2} + x_{0}^{2} x_{1}^{4}}$"
+      ],
       "text/plain": [
-       "{s(0): (x0**2 + x1**2)*(log(sqrt(x0**2 + x1**2)) - 1),\n",
-       " s(1): -2*x0*(log(sqrt(x0**2 + x1**2)) - 1) - x0,\n",
-       " s(2): -2*x0*(log(sqrt(x0**2 + x1**2)) - 1) - x0,\n",
-       " s(3): -2*x0*(log(sqrt(x0**2 + x1**2)) - 1) - x0,\n",
-       " s(4): 2*(-24*x0**4/(x0**2 + x1**2)**2 + 8*x0**2*(4*x0**2/(x0**2 + x1**2) - 3)/(x0**2 + x1**2) + 12*x0**2/(x0**2 + x1**2) + 3)/(x0**2 + x1**2),\n",
-       " s(5): -4*x0*(-24*x0**4/(x0**2 + x1**2)**2 + 40*x0**2/(x0**2 + x1**2) - 15)/(x0**2 + x1**2)**2,\n",
-       " s(6): 12*(-320*x0**6/(x0**2 + x1**2)**3 + 360*x0**4/(x0**2 + x1**2)**2 + 24*x0**2*(16*x0**4/(x0**2 + x1**2)**2 - 20*x0**2/(x0**2 + x1**2) + 5)/(x0**2 + x1**2) - 60*x0**2/(x0**2 + x1**2) - 5)/(x0**2 + x1**2)**2,\n",
-       " s(7): -48*x0*(-160*x0**6/(x0**2 + x1**2)**3 + 336*x0**4/(x0**2 + x1**2)**2 - 210*x0**2/(x0**2 + x1**2) + 35)/(x0**2 + x1**2)**3}"
+       "-(-2*x0**2 + 6*x1**2)*s(3)/(x0**3 + x0*x1**2) - (-x0**4 + 6*x0**2*x1**2 + 15*x1**4)*s(1)/(x0**7 + 2*x0**5*x1**2 + x0**3*x1**4) + (x0**4 - 6*x0**2*x1**2 - 15*x1**4)*s(2)/(x0**6 + 2*x0**4*x1**2 + x0**2*x1**4)"
       ]
      },
-     "execution_count": 6,
+     "execution_count": 19,
      "metadata": {},
      "output_type": "execute_result"
     }
    ],
    "source": [
-    "subs_dict"
+    "r.subs(n, 4)"
    ]
   },
   {