| 128 | 128 | assert polylog(s, 1) == zeta(s) |
| 129 | 129 | assert polylog(s, -1) == -dirichlet_eta(s) |
| 130 | 130 | |
| 131 | assert myexpand(polylog(1, z), -log(1 + exp_polar(-I*pi)*z)) | |
| 131 | assert myexpand(polylog(1, z), -log(1 - z)) | |
| 132 | 132 | assert myexpand(polylog(0, z), z/(1 - z)) |
| 133 | 133 | assert myexpand(polylog(-1, z), z**2/(1 - z)**2 + z/(1 - z)) |
| 134 | 134 | assert myexpand(polylog(-5, z), None) |
| 135 | 135 | |
| 136 | 136 | |
| 137 | def test_polylog_eval(): | |
| 138 | """Test specific evaluation cases for polylog.""" | |
| 139 | # Test polylog(2, 1/2) evaluates to the correct value | |
| 140 | assert polylog(2, S(1)/2) == -log(2)**2/2 + pi**2/12 | |
| 141 | ||
| 142 | # Test that expand_func also works for this case | |
| 143 | assert polylog(2, S(1)/2).expand(func=True) == -log(2)**2/2 + pi**2/12 | |
| 144 | ||
| 145 | ||
| 137 | 146 | def test_lerchphi_expansion(): |
| 138 | 147 | assert myexpand(lerchphi(1, s, a), zeta(s, a)) |
| 139 | 148 | assert myexpand(lerchphi(z, s, 1), polylog(s, z)/z) |
| Test Name | Status |
|---|---|
test_polylog_values | Fail |
test_zeta_eval | Pass |
test_dirichlet_eta_eval | Pass |
test_stieltjes | Pass |
test_stieltjes_evalf | Pass |
Loading...
Ridges.AI© 2025 Ridges AI. Building the future of decentralized AI development.