Large difference in output of optimize_scalar and brent
"ashwin .D" <[email protected]> Tue, 24 May 2022 19:39:31 +0530
| Newsgroups | gmane.comp.python.scientific.user |
|---|---|
| Message-ID | <CAH0LXy75rFnYrXSXUDPSkp8vbRrgwFX2SG43bvPLB9Yp=sgYAA@mail.gmail.com> |
--00000000000011671f05dfc2829c
Content-Type: multipart/alternative; boundary="00000000000011671c05dfc28299"
--00000000000011671c05dfc28299
Content-Type: text/plain; charset="UTF-8"
Hello,
I have enclosed a minimum reproducible example with data (as SO
demands) and I am trying to use two methods in the scipy optimize package -
https://docs.scipy.org/doc/scipy/reference/optimize.html
One is minimize_scalar and the other is Brent's method. But both the
outputs have large differences. Can anybody explain why ?
Also I get a runtime error while running Brent-
test.py:29: RuntimeWarning: invalid value encountered in power
ls = np.sum((x1[1:] - x1[:-1]*np.sign(a)*a**dt)**2)
Additionally the Brent one has another one if I specify a triple interval
ala - dum = brent(func,brack=(-1,0,1),tol=3.0e-8)
ValueError : Not a bracketing interval.
https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.brent.html
But docs say that I can specify a triplet in the bracketing interval.
I could have asked this on SO (and I still can) but this is more of a
discussion oriented question so hence posted here. If a bug please advise
appropriately as I can open on github.
Regards,
Ashwin.
--00000000000011671c05dfc28299
Content-Type: text/html; charset="UTF-8"
Content-Transfer-Encoding: quoted-printable
<div dir=3D"ltr">Hello,<div>=C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0 =C2=A0I have=
enclosed a minimum reproducible example with data (as SO demands) and I am=
trying to use two methods in the scipy optimize package -=C2=A0<a href=3D"=
https://docs.scipy.org/doc/scipy/reference/optimize.html">https://docs.scip=
y.org/doc/scipy/reference/optimize.html</a></div><div><br></div><div>One is=
minimize_scalar and the other is Brent's method. But both the outputs=
=C2=A0have large differences. Can anybody explain why ?</div><div><br></div=
><div>Also I get a runtime error while running Brent-=C2=A0</div><div><br><=
/div>test.py:29: RuntimeWarning: invalid value encountered in power<br><div=
>=C2=A0 ls =3D np.sum((x1[1:] - x1[:-1]*np.sign(a)*a**dt)**2)=C2=A0</div><d=
iv><br></div><div>Additionally the=C2=A0Brent one has another one if I spec=
ify a triple interval</div><div>ala -=C2=A0 dum =3D brent(func,brack=3D(-1,=
0,1),tol=3D3.0e-8)</div><div>ValueError=C2=A0 : Not a bracketing interval.=
=C2=A0</div><div><br></div><div><a href=3D"https://docs.scipy.org/doc/scipy=
/reference/generated/scipy.optimize.brent.html">https://docs.scipy.org/doc/=
scipy/reference/generated/scipy.optimize.brent.html</a> But docs say that I=
can specify a triplet in the bracketing interval.=C2=A0<br></div><div><br>=
</div><div>I could have asked this on SO (and I still can) but this is more=
of a discussion oriented question so hence posted here. If a bug please ad=
vise appropriately as I can open on github.=C2=A0</div><div><br></div><div>=
Regards,</div><div>Ashwin.=C2=A0</div></div>
--00000000000011671c05dfc28299--
--00000000000011671f05dfc2829c
Content-Type: application/octet-stream; name="x.dat"
Content-Disposition: attachment; filename="x.dat"
Content-Transfer-Encoding: base64
Content-ID: <f_l3k8a08l1>
X-Attachment-Id: f_l3k8a08l1
MC4wMDAwMDBFKzAwICAwLjIxMTU0NkUrMDENCjAuMTQ4NTY1RSswMSAgMC44NjE1NzBFKzAwDQow
LjI2MDQ3MkUrMDEgIC0uMjU2NjIzRSswMA0KMC41NjY0MDNFKzAxICAwLjQzNjg2OEUrMDENCjAu
NjIwODQ3RSswMSAgMC4yMzQ1MTJFKzAxDQowLjk1NTgwNkUrMDEgIDAuMTkxNjA1RSswMQ0KMC4x
MDQ3ODVFKzAyICAwLjM0Mzc1NUUrMDENCjAuMTA3NjgwRSswMiAgMC4zNDY0OTVFKzAxDQowLjEw
OTg3MEUrMDIgIDAuMTg5NjY3RSswMQ0KMC4xMTI0NTBFKzAyICAwLjE1MjYwMUUrMDENCjAuMTQy
MzA0RSswMiAgLS42NzY1ODBFKzAwDQowLjE2Mjg3NUUrMDIgIDAuNzQ4NTkzRSswMA0KMC4xNzAx
NjFFKzAyICAtLjEwMDY2NUUrMDENCjAuMTc5NTc1RSswMiAgLS4xNzM1MDdFKzAxDQowLjIwMTg3
NUUrMDIgIDAuMTgzMTA5RSswMQ0KMC4yMTM1NzhFKzAyICAwLjIzMDk1NUUrMDENCjAuMjE1MTc0
RSswMiAgMC44MzM4MzBFKzAwDQowLjIxNzY4NUUrMDIgIC0uMjgwODAyRS0wMQ0KMC4yMzY1MzlF
KzAyICAtLjE4MTAyNUUrMDENCjAuMjQ1MjgzRSswMiAgMC4xNzYwODdFKzAxDQowLjI0Nzc1MkUr
MDIgIDAuMTUzMTQzRSswMQ0KMC4yNTAzNjNFKzAyICAwLjEwNzAyNUUrMDENCjAuMjcyNDY2RSsw
MiAgLS43OTE2MTFFKzAwDQowLjI3NTE4NkUrMDIgIC0uMTE2ODMzRSswMQ0KMC4yOTExMzlFKzAy
ICAtLjc2NzA4MkUrMDANCjAuMjkzMjI2RSswMiAgMC41NzE2ODRFKzAwDQowLjMyMDU0N0UrMDIg
IDAuNTA5MDY3RSswMA0KMC4zNDI2OTNFKzAyICAwLjM2Mzg1OUUrMDANCjAuMzUzMzk2RSswMiAg
MC4xMjIzNTVFKzAwDQowLjM3MjkwN0UrMDIgIC0uMzM2MjcyRSswMA0KMC40MDA1MjBFKzAyICAw
LjE0ODgwOUUrMDENCjAuNDE0MjAxRSswMiAgMC45OTA5ODZFKzAwDQowLjQyNTgzMUUrMDIgIDAu
MjIzODU5RSswMA0KMC40MjgyNjNFKzAyICAwLjE2NTQ5MEUrMDANCjAuNDU0NDQyRSswMiAgMC4y
OTU4ODdFKzAxDQowLjQ3MDU0OUUrMDIgIDAuOTE1NDM5RSswMA0KMC40NzE2ODNFKzAyICAwLjI4
MjMyMEUrMDANCjAuNDc1NDgyRSswMiAgLS4zMTA3MzFFKzAwDQowLjQ4MjkxMEUrMDIgIC0uMTE5
MjY0RSswMQ0KMC40ODQzMzNFKzAyICAtLjQ5OTg1M0UtMDENCjAuNDg4MTExRSswMiAgMC4yMTE1
ODhFKzAwDQowLjUwNDc3MEUrMDIgIDAuMjIyMzAzRSswMQ0KMC41MDkwNTNFKzAyICAwLjE2Nzk1
MEUrMDENCjAuNTExNTUxRSswMiAgMC4yNjAyNzZFKzAxDQowLjUxODY3MEUrMDIgIDAuNzM5NDM3
RSswMA0KMC41MjU0NzJFKzAyICAtLjExMDA2NEUrMDENCjAuNTQwNjQ2RSswMiAgLS43NTgwMDdF
KzAwDQowLjU1Mzc5N0UrMDIgIC0uODk0NzQwRSswMA0KMC41NjE3NDdFKzAyICAtLjE0OTQyNEUr
MDENCjAuNTY2MzI0RSswMiAgLS45NzQyMzlFKzAwDQowLjU3NTc0MUUrMDIgIC0uMjQ4MDkzRSsw
MQ0KMC41OTQ5NTNFKzAyICAwLjE5MDE3OEUrMDENCjAuNjEwMTU0RSswMiAgLS42MzQ4MDlFKzAw
DQowLjYyNTE3MEUrMDIgIC0uMTUzMjE1RSswMQ0KMC42NzU1NzRFKzAyICAtLjE4MjE0MUUrMDEN
CjAuNjkzNDc2RSswMiAgMC40Nzk5NjlFKzAwDQowLjcwOTg0M0UrMDIgIDAuMTkzNTQ4RSswMQ0K
MC43MTI5MTFFKzAyICAwLjEzMjI5N0UrMDANCjAuNzE0ODcwRSswMiAgLS41NDY2MDdFKzAwDQow
LjcxNTk0NEUrMDIgIC0uNzkyMjk4RSswMA0KMC43MzQ2NDFFKzAyICAwLjIwOTIyMUUrMDANCjAu
NzM2NTc4RSswMiAgLS4zMDQwNzVFLTAxDQowLjc0MjMxN0UrMDIgIDAuNzE5OTcwRSswMA0KMC43
NDg2MTJFKzAyICAwLjE2NTUwMEUrMDENCjAuNzg5NjQxRSswMiAgLS4zMjA2OTZFKzAwDQowLjgx
MjE5MEUrMDIgIDAuMjUxNzI1RSswMQ0KMC44MTU0NDhFKzAyICAwLjI3OTMxOUUrMDENCjAuODI4
MjkyRSswMiAgLS4xMzgzODVFKzAxDQowLjg1NTkwN0UrMDIgIDAuMjUzNjM2RSswMQ0KMC44NjI3
NTNFKzAyICAwLjExMzQxMEUrMDANCjAuODY0Mjc3RSswMiAgMC4xNzA1OTNFKzAwDQowLjg2Njk4
NkUrMDIgIC0uMTI0NTMyRSswMQ0KMC44ODQxNzhFKzAyICAtLjY5MjA1OEUrMDANCjAuODg2MDgz
RSswMiAgLS45NjkxNzlFKzAwDQowLjkwMjY1MkUrMDIgIDAuMTE0NzA5RSswMQ0KMC45MDQwMzlF
KzAyICAwLjU3MzkyM0UrMDANCjAuOTIyNjI3RSswMiAgLS4xNjEyNTVFKzAxDQowLjkyNTgyNUUr
MDIgIC0uMTkyMjI5RSswMQ0KMC45NDc4NzNFKzAyICAtLjE0ODY3MEUrMDENCjAuOTUzMDc3RSsw
MiAgMC4xNTUyODdFKzAxDQowLjk2MTYwNEUrMDIgIDAuMTgxNTk4RSswMQ0KMC45NzQ3NThFKzAy
ICAtLjI3OTY4N0UrMDANCjAuOTc5MzY5RSswMiAgLS4xOTg3MzJFKzAxDQowLjEwMDk1OEUrMDMg
IDAuMTk5NjUwRSswMQ0KMC4xMDE5NTBFKzAzICAwLjczNzM4NkUrMDANCjAuMTAyMzE5RSswMyAg
MC41NDk0MDhFKzAwDQowLjEwMzE5MEUrMDMgIC0uNTgyNjY5RSswMA0KMC4xMDMzNjNFKzAzICAt
LjUzNDI1OUUrMDANCjAuMTAzNTMwRSswMyAgMC4zNjAzNzlFKzAwDQowLjEwNDY5NEUrMDMgIDAu
MTgyMzg0RSswMQ0KMC4xMDYzMjNFKzAzICAwLjIyMzkzOEUrMDENCjAuMTA2OTM4RSswMyAgLS41
MDkzMzZFKzAwDQowLjEwODQ2M0UrMDMgIC0uMTIxMjA0RSswMQ0KMC4xMDk1MTZFKzAzICAwLjE1
MDI2N0UrMDENCjAuMTEyMDE3RSswMyAgLS40NTcxNzFFKzAwDQowLjExMzA4OEUrMDMgIC0uMjA0
NTQ1RSswMQ0KMC4xMTM2ODdFKzAzICAtLjE4NDE3OEUrMDENCjAuMTE0MTc0RSswMyAgLS4zMzAw
MzhFKzAwDQowLjExNTUwOUUrMDMgIDAuMjA4MzU0RSswMQ0KMC4xMTYxOThFKzAzICAtLjE4NDI5
NkUrMDANCjAuMTE4MTE1RSswMyAgLS4yODUyMDNFKzAxDQowLjExODQzOEUrMDMgIC0uMzE4MjYz
RSswMQ0KMC4xMTkxNzVFKzAzICAtLjM5ODYzNUUrMDANCjAuMTE5NDg5RSswMyAgMC41MTA4NzZF
KzAwDQowLjEyMDc0NkUrMDMgIDAuMTAzNjM1RSswMQ0KMC4xMjEwOTRFKzAzICAwLjU5OTYyN0Ur
MDANCjAuMTIzNTc1RSswMyAgLS4yMTEwOTZFKzAxDQowLjEyMzg2NEUrMDMgIC0uMTQxMzEwRSsw
MQ0KMC4xMjUwNTNFKzAzICAtLjY2MjY1N0UrMDANCjAuMTI2MDA4RSswMyAgMC43NTM4MjhFKzAw
DQowLjEyNjk1MUUrMDMgIC0uMTYyOTAwRSswMQ0KMC4xMjg1OTBFKzAzICAtLjQ4MTgxOEUrMDAN
CjAuMTI4ODAzRSswMyAgLS41NzY1NjNFKzAwDQowLjEyOTc1MUUrMDMgIDAuMTc2ODAxRSswMQ0K
MC4xMzAwNjBFKzAzICAwLjIyMjg0NkUrMDENCjAuMTMxMTgzRSswMyAgLS4zMjc3OTZFKzAwDQow
LjEzMTM3M0UrMDMgIDAuMzY1NzY4RSswMA0KMC4xMzI3MTRFKzAzICAwLjMwMDYzOUUrMDANCjAu
MTMzNDQxRSswMyAgMC42MjQ4MzRFKzAwDQowLjEzNjMxM0UrMDMgIDAuMTU5MTMzRSswMQ0KMC4x
MzgxNDJFKzAzICAtLjY1MzAzNUUrMDANCjAuMTM5MjI2RSswMyAgMC4yNjU2MjlFKzAwDQowLjE0
MTA3MkUrMDMgIDAuMTY2MjA5RSswMQ0KMC4xNDE0NTlFKzAzICAwLjEzNDM2MEUrMDENCjAuMTQz
NDU4RSswMyAgLS43MTM0MzZFKzAwDQowLjE0NDE0N0UrMDMgIDAuMTU1MTk2RSswMQ0KMC4xNDUy
MDBFKzAzICAwLjI1NzU5NEUrMDENCjAuMTQ1ODA5RSswMyAgMC4xNjk0MjNFKzAxDQowLjE0NjY5
MUUrMDMgIDAuNzMyOTEyRSswMA0KMC4xNDY4NDhFKzAzICAwLjcyNDc0OUUtMDINCjAuMTQ4Njgw
RSswMyAgLS4xMjI3MTVFKzAxDQowLjE1MDgwNEUrMDMgIDAuMTE3MjAyRSswMQ0KMC4xNTEwMzRF
KzAzICAwLjIwMDE4OUUrMDENCjAuMTUxMzY3RSswMyAgMC4xOTU5NDdFKzAxDQowLjE1MjcxOUUr
MDMgIC0uODIxMzA0RSswMA0KMC4xNTM2NjdFKzAzICAtLjIxNzY3NkUrMDENCjAuMTU1MDEyRSsw
MyAgLS44MjQzOTlFKzAwDQowLjE1NTMyNUUrMDMgIC0uMzUwOTEzRSswMA0KMC4xNTY2MDlFKzAz
ICAtLjMyMTQ5NEUrMDENCjAuMTU5MDg5RSswMyAgLS4yMjAzMjNFKzAxDQowLjE1OTI4NkUrMDMg
IC0uMTQwMjMyRSswMQ0KMC4xNTk2MDFFKzAzICAtLjcyNDQyNUUrMDANCjAuMTYwNDU3RSswMyAg
MC45NjQ3NDFFKzAwDQowLjE2MTM3NEUrMDMgIDAuMzEwMTk0RSswMA0KMC4xNjE0ODZFKzAzICAw
Ljc1NzQ5MEUrMDANCjAuMTYyOTY1RSswMyAgLS4zNjA4NDJFKzAwDQowLjE2NTM4MEUrMDMgIDAu
MTMyNjkyRSswMQ0KMC4xNjU1NTFFKzAzICAwLjE1MTkyNUUtMDENCjAuMTY3MTk4RSswMyAgLS4x
ODU2MzlFKzAxDQowLjE2ODMxMUUrMDMgIC0uODcyMTM5RSswMA0KMC4xNjg2MTBFKzAzICAtLjUy
ODE1OEUrMDANCjAuMTY5MDYwRSswMyAgMC4xMTM1MzJFKzAxDQowLjE3MTIzMUUrMDMgIDAuNTA1
Mjc0RSswMA0KMC4xNzE3NDRFKzAzICAwLjExNzQ5NUUrMDENCjAuMTcyMzc0RSswMyAgMC4xNTYy
NDhFKzAxDQowLjE3MjUyN0UrMDMgIDAuMTU0Nzg5RSswMQ0KMC4xNzI2MzRFKzAzICAwLjEwMjkx
NkUrMDENCjAuMTczNjE3RSswMyAgMC4yOTMxMTFFKzAwDQowLjE3NDExM0UrMDMgIDAuMzIyNTk3
RSswMA0KMC4xNzQ5NDJFKzAzICAwLjUwMTg2MEUrMDANCjAuMTc1ODA1RSswMyAgMC4zMjc1ODJF
KzAwDQowLjE4MDcyM0UrMDMgIDAuMjQ1Nzg0RSswMQ0KMC4xODE1NDJFKzAzICAwLjExMDQ2NkUr
MDENCjAuMTgxNjk5RSswMyAgMC41OTY4MzZFKzAwDQowLjE4MjM4N0UrMDMgIC0uMTc4MDUyRSsw
MQ0KMC4xODMxOTJFKzAzICAtLjU0NTg2OEUrMDANCjAuMTgzNDE1RSswMyAgLS42MDY3MTBFKzAw
DQowLjE4NDc3NkUrMDMgIDAuOTA3MjQ3RSswMA0KMC4xODU5ODFFKzAzICAwLjUwMDM4OUUtMDEN
CjAuMTg3ODUyRSswMyAgLS4yNDMzNTFFKzAxDQowLjE5MDYwN0UrMDMgIDAuMjI5ODE2RSswMQ0K
MC4xOTI4MTZFKzAzICAtLjM2ODk0NkUrMDENCjAuMTkzMDY0RSswMyAgLS4zNDY2MDNFKzAxDQow
LjE5NDg5NUUrMDMgIC0uOTM3NTk5RSswMA0KMC4xOTYyNzZFKzAzICAtLjY0MTcyMUUrMDANCjAu
MTk3MTQ1RSswMyAgLS41OTQzODVFKzAwDQowLjE5Nzc5NUUrMDMgIC0uMjQ3MzI4RSswMQ0KMC4x
OTgwNzFFKzAzICAtLjI1MTQ3NEUrMDENCjAuMTk5MTE2RSswMyAgMC45NzQwNzJFKzAwDQowLjE5
OTQ4MEUrMDMgIDAuODQwNjM5RSswMA0KMC4yMDAzOTZFKzAzICAwLjI0MDgxNUUrMDENCjAuMjAw
NTI4RSswMyAgMC4yODc4MTdFKzAxDQowLjIwMDY3OEUrMDMgIDAuMjg0MTE3RSswMQ0KMC4yMDEy
NTNFKzAzICAwLjE3NDI5MEUrMDENCjAuMjAxNzI3RSswMyAgMC4zMzM3ODBFLTAxDQowLjIwMjI0
NEUrMDMgIC0uMjQ3MTI5RSswMA0KMC4yMDIzOThFKzAzICAtLjc2NzM4NUUrMDANCjAuMjAyNzAz
RSswMyAgLS43MDcyMzRFKzAwDQowLjIwMjg0NUUrMDMgIC0uOTkwODcyRS0wMQ0KMC4yMDM1MzZF
KzAzICAwLjE1MTYxM0UrMDENCjAuMjA0MzQ2RSswMyAgMC4yMzczNTRFKzAxDQowLjIwNDkzNEUr
MDMgIDAuMjAzNTQ2RSswMQ0KMC4yMDY1MDVFKzAzICAwLjE4NTczMkUrMDENCjAuMjA2NzYyRSsw
MyAgMC4xMzgxODhFKzAxDQowLjIwODI3NkUrMDMgIDAuMzg1MjEzRSswMA0KMC4yMDk2MDRFKzAz
ICAwLjE1NjA5MEUrMDENCjAuMjExOTExRSswMyAgLS44MDg5NTZFKzAwDQowLjIxNDY0MUUrMDMg
IC0uMTEzNjIwRSswMQ0KMC4yMTUxNDlFKzAzICAwLjI3ODgyOEUrMDANCjAuMjE2MjQ4RSswMyAg
MC4xMzgzNTdFKzAxDQowLjIxNjc1NUUrMDMgIDAuMjIxNTU5RSswMQ0KMC4yMTc5MzdFKzAzICAw
LjM1Njk1M0UrMDANCjAuMjE5NTU3RSswMyAgLS4zMzMwMjJFLTAyDQowLjIyMDA5MUUrMDMgIDAu
Mzg5OTU1RSswMA0KMC4yMjA2MTZFKzAzICAwLjE3NjYyNUUrMDENCjAuMjIxMjM0RSswMyAgMC4x
NTYxNDdFKzAxDQowLjIyMjY3NkUrMDMgIC0uMjM5NTk1RSswMQ0KMC4yMjM1NDZFKzAzICAtLjE5
MTI2N0UrMDENCjAuMjI0MzE1RSswMyAgLS4xNjA4MTNFKzAwDQowLjIyNzMxMUUrMDMgIC0uMjQ0
ODgzRSswMQ0KMC4yMjc0MzdFKzAzICAtLjI0ODc2MUUrMDENCjAuMjI4MzM5RSswMyAgLS45NTgy
MjBFKzAwDQowLjIyOTc3NkUrMDMgIDAuMzI3NjYwRSswMA0KMC4yMzAzODRFKzAzICAwLjEyODUz
NkUrMDANCjAuMjMwOTI4RSswMyAgLS44MTI3NThFLTAxDQowLjIzMTIwMEUrMDMgIDAuNjMxOTg5
RSswMA0KMC4yMzE2MzZFKzAzICAwLjY5ODM4MUUrMDANCjAuMjMxNzU4RSswMyAgMC4zNTg0MTVF
KzAwDQowLjIzMTg3NUUrMDMgIDAuMjcwOTc2RSswMA0KMC4yMzQ5NzJFKzAzICAwLjk5MDE1NkUr
MDANCjAuMjM2NTg2RSswMyAgMC4yOTQ5NDlFKzAxDQowLjIzNzMyOEUrMDMgIDAuODU3NDQyRSsw
MA0KMC4yMzc2NTBFKzAzICAtLjIyMzAwM0UrMDANCjAuMjM4MDQwRSswMyAgLS41MjAzODBFKzAw
DQowLjIzOTE4NkUrMDMgIDAuMTQ3MjgxRSswMQ0KMC4yNDIxMTZFKzAzICAtLjk5ODY5OEUtMDEN
CjAuMjQyNTA1RSswMyAgLS4xNjYzNjFFKzAwDQowLjI0MzI2OEUrMDMgIC0uMTQ5MjQ0RSswMA0K
MC4yNDQ0MzNFKzAzICAwLjEzNTkwMkUrMDENCjAuMjQ0NjM1RSswMyAgMC4xMjQxMjdFKzAxDQow
LjI0NDg2MkUrMDMgIDAuMjIyMTAzRSswMQ0KMC4yNDY3NzFFKzAzICAtLjUzMjU1OUUrMDANCjAu
MjQ3MDk5RSswMyAgLS4xOTQ1MTFFKzAxDQowLjI0NzI5M0UrMDMgIC0uMTg0NzU4RSswMQ0KMC4y
NDg0NDdFKzAzICAwLjUyOTM2MEUrMDANCjAuMjQ5MzIyRSswMyAgLS4xMTUzMTVFKzAwDQowLjI0
OTk5M0UrMDMgIDAuMTg3MjcyRSswMQ0KMC4yNTA2MDdFKzAzICAwLjI5Nzg4OEUrMDENCjAuMjUx
MTM1RSswMyAgMC44MTgyODdFKzAwDQowLjI1MTg5NkUrMDMgIC0uNzQzMjc1RSswMA0KMC4yNTIx
OTFFKzAzICAtLjk4MjQwNkUrMDANCjAuMjUzNzY0RSswMyAgLS4xMTMwMzZFKzAxDQowLjI1NTYx
OEUrMDMgIDAuMTYzMTc0RSswMQ0KMC4yNTYwNDlFKzAzICAwLjczMzAwOEUrMDANCjAuMjU2MzA1
RSswMyAgMC4yMzQ4OTdFKzAwDQowLjI1Nzk2NUUrMDMgIC0uMjM1MTU4RSswMQ0KMC4yNTkyNTVF
KzAzICAtLjY1Nzg0MUUrMDANCjAuMjYwMTI5RSswMyAgLS40OTMyMzFFKzAwDQowLjI2MDk4NkUr
MDMgIDAuMjEyNzU0RSswMA0KMC4yNjMzNTJFKzAzICAtLjE4NzMwNUUrMDENCjAuMjYzOTIxRSsw
MyAgLS43NDUwMTBFKzAwDQowLjI2NDQ5OUUrMDMgIC0uNTI5MjgzRS0wMQ0KMC4yNjQ3NzJFKzAz
ICAtLjI1MzA2MkUrMDANCjAuMjY4MzAxRSswMyAgLS4xNDAwMTNFKzAxDQowLjI2ODc1M0UrMDMg
IC0uNzUxMjUyRSswMA0KMC4yNzAzNDFFKzAzICAwLjI0ODQwMkUrMDENCjAuMjcwNTc2RSswMyAg
MC4yOTgzNzFFKzAxDQowLjI3MDc2NkUrMDMgIDAuMjY4Njk5RSswMQ0KMC4yNzE5MThFKzAzICAw
LjE4MDU1NEUrMDANCjAuMjczOTQ5RSswMyAgMC41NjA5NjhFKzAwDQowLjI3NDQ2MEUrMDMgIDAu
MTI3MzE4RSswMQ0KMC4yNzY3NjdFKzAzICAwLjIyMDI0MUUrMDENCjAuMjc3OTU0RSswMyAgMC45
MTY4NTRFKzAwDQowLjI3ODM3NkUrMDMgIC0uNzIzNTk3RSswMA0KMC4yNzk5NzhFKzAzICAwLjIy
MjU2NkUrMDENCjAuMjgxOTc4RSswMyAgMC45MjI1NDZFKzAwDQowLjI4Mzk0OUUrMDMgIDAuOTQ2
MTc4RSswMA0KMC4yODUzMDJFKzAzICAwLjIzODExMEUrMDENCjAuMjg1NzYzRSswMyAgMC4yMjc5
MzdFKzAxDQowLjI4NjYyNEUrMDMgIC0uMTY4MTUwRSswMA0KMC4yODg4OTZFKzAzICAtLjE4Nzg0
OUUrMDENCjAuMjg5MTU4RSswMyAgLS43MjAzMzlFKzAwDQowLjI5MTI5M0UrMDMgIDAuNzczNjg5
RSswMA0KMC4yOTE2OTJFKzAzICAtLjEwNDkyNkUrMDENCjAuMjkzMDA0RSswMyAgLS4xNjQ1NjFF
KzAxDQowLjI5MzM4MkUrMDMgIC0uMjU5NzEwRSswMQ0KMC4yOTM3OTNFKzAzICAtLjE3MTI5N0Ur
MDENCjAuMjk3NDk3RSswMyAgLS4yNTkzNDVFKzAxDQowLjI5ODg0NUUrMDMgIC0uMTM2NTQ4RSsw
MQ0KMC4zMDAxMjlFKzAzICAwLjE0NzEyMEUrMDENCjAuMzAwMjczRSswMyAgMC4xOTMwMjdFKzAx
DQowLjMwMjI5M0UrMDMgIC0uNTIxNDA0RSswMA0KMC4zMDI1NTNFKzAzICAtLjk2MDI5MkUrMDAN
CjAuMzAzMjMyRSswMyAgLS4yODYzNjJFKzAxDQowLjMwNDUxMkUrMDMgIDAuMTMyODM1RSswMQ0K
MC4zMDY3ODdFKzAzICAwLjcxOTA1MEUrMDANCjAuMzA3Njg2RSswMyAgLS44MjY0MjVFKzAwDQow
LjMwNzg4OEUrMDMgIC0uMTMxOTkzRSswMA0KMC4zMDk2NzBFKzAzICAwLjE5NDgyMEUrMDENCjAu
MzEwODEyRSswMyAgMC4yOTMzMjRFKzAxDQowLjMxMTM1NkUrMDMgIDAuMjAwNzA4RSswMQ0KMC4z
MTE2NThFKzAzICAwLjIyNDI5MEUrMDENCjAuMzEyODY4RSswMyAgMC43MjQwNzJFKzAwDQowLjMx
MzU5N0UrMDMgIDAuMzc4MzY3RSswMA0KMC4zMTM5ODlFKzAzICAtLjQzOTU3N0UrMDANCjAuMzE0
NjcwRSswMyAgMC4xMTY3NzFFKzAxDQowLjMxNTY0OUUrMDMgIDAuMTMwNDQ4RSswMQ0KMC4zMTg2
MTZFKzAzICAtLjE4MTMyMEUrMDENCjAuMzE4NzE5RSswMyAgLS4xNzEwNjRFKzAxDQowLjMxOTMw
MUUrMDMgIC0uMTM2NzE1RSswMQ0K
--00000000000011671f05dfc2829c
Content-Type: text/x-python; charset="US-ASCII"; name="test.py"
Content-Disposition: attachment; filename="test.py"
Content-Transfer-Encoding: base64
Content-ID: <f_l3k8a07p0>
X-Attachment-Id: f_l3k8a07p0
aW1wb3J0IG51bXB5IGFzIG5wCmZyb20gc2NpcHkub3B0aW1pemUgaW1wb3J0IGJyZW50LG1pbmlt
aXplX3NjYWxhcgppbXBvcnQgbWF0aAoKZGVmIG1haW4oKToKCiAgICBkYXRhID0gbnAubG9hZHR4
dCgieC5kYXQiLHNraXByb3dzPTApCiAgICB0ID0gZGF0YVs6LDBdIAogICAgeCA9IGRhdGFbOiwx
XQogICAgZHVtID0gbWlubHModCx4KQogICAgCmRlZiBtaW5scyh0LHgpOgogICAgYSA9IDEvbWF0
aC5lCiAgICBkdCA9IG5wLmRpZmYodCkKICAgIGRlZiBmdW5jKGEpOgogICAgICAgIHRyeToKICAg
ICAgICAgICAgbHMgPSBucC5zdW0oKHhbMTpdIC0geFs6LTFdKmEqKmR0KSoqMikKICAgICAgICBl
eGNlcHQgV2FybmluZzoKICAgICAgICAgICAgcHJpbnQodHJhY2ViYWNrLmZvcm1hdF9leGMoKSkK
ICAgICAgICByZXR1cm4gbHMKCiAgICBkdW0gPSBicmVudChmdW5jLGJyYWNrPSgwLDEpLHRvbD0z
LjBlLTgpCiAgICBkdW0xID0gbWluaW1pemVfc2NhbGFyKGZ1bmMsIGJvdW5kcz1bMCwgMV0sIG1l
dGhvZD0nYm91bmRlZCcpLngKICAgIHByaW50KGR1bSxkdW0xKQogICAgcmV0dXJuIGR1bQoKbWFp
bigpCg==
--00000000000011671f05dfc2829c
Content-Type: text/plain; charset="us-ascii"
MIME-Version: 1.0
Content-Transfer-Encoding: 7bit
Content-Disposition: inline
_______________________________________________
SciPy-User mailing list -- [email protected]
To unsubscribe send an email to [email protected]
https://mail.python.org/mailman3/lists/scipy-user.python.org/
Member address: [email protected]
--00000000000011671f05dfc2829c--