|
542 | 542 | { |
543 | 543 | "cell_type": "code", |
544 | 544 | "execution_count": 10, |
545 | | - "id": "ce108dbf", |
| 545 | + "id": "39f34f43", |
546 | 546 | "metadata": {}, |
547 | 547 | "outputs": [], |
548 | 548 | "source": [ |
|
555 | 555 | " \n", |
556 | 556 | " U00 = np.sqrt((2*L+1)/(4*np.pi))\n", |
557 | 557 | " T01 = clebsch_gordan(L1, L2, L, 0, 0, 0).evalf()\n", |
558 | | - " P0L1 = pysh.legendre.legendre_lm (L1, 0, cos(theta1), 'unnorm', -1, 1)\n", |
559 | | - " P0L2 = pysh.legendre.legendre_lm (L2, 0, cos(theta2), 'unnorm', -1, 1)\n", |
| 558 | + " P0L1 = pysh.legendre.legendre_lm (L1, 0, cos(theta1), 'unnorm', -1, 0)\n", |
| 559 | + " P0L2 = pysh.legendre.legendre_lm (L2, 0, cos(theta2), 'unnorm', -1, 0)\n", |
560 | 560 | " T1 = T01*P0L1*P0L2\n", |
561 | 561 | " Total += T1\n", |
562 | 562 | " M=1\n", |
563 | 563 | " while (M<=M_max):\n", |
564 | 564 | " T02 = clebsch_gordan(L1, L2, L, M, -M, 0).evalf()\n", |
565 | | - " PmL1 = pysh.legendre.legendre_lm (L1, M, cos(theta1), 'unnorm', -1, 1)\n", |
566 | | - " PmL2 = pysh.legendre.legendre_lm (L2, M, cos(theta2), 'unnorm', -1, 1)\n", |
| 565 | + " PmL1 = pysh.legendre.legendre_lm (L1, M, cos(theta1), 'unnorm', -1, 0)\n", |
| 566 | + " PmL2 = pysh.legendre.legendre_lm (L2, M, cos(theta2), 'unnorm', -1, 0)\n", |
567 | 567 | " T2 = T02*PmL1*PmL2\n", |
568 | 568 | " Total += 2.0*pow((-1),M)*T2*(cos(M*phi))\n", |
569 | 569 | " M+=1\n", |
|
0 commit comments