Re: How to add an if statement in an objective function in scipy.optimize
Travis Vaught <[email protected]> Tue, 22 Nov 2022 10:31:44 -0600
| Newsgroups | gmane.comp.python.scientific.user |
|---|---|
| Message-ID | <[email protected]> |
DQoNCj4gT24gTm92IDIyLCAyMDIyLCBhdCA5OjQwIEFNLCBaZWluYSBBYnUtQWlzaGVoIDx6ZWlu YS5hYnVhaXNoZWhAZ21haWwuY29tPiB3cm90ZToNCj4gDQo+IEkgaGF2ZSBkYXRhICh4KSB0aGF0 IHJlcHJlc2VudHMgc29tZSBhbmdsZXMgdGhhdCBsb29rcyBsaWtlOiBbMjAsIDMwLCA0OCwgNDYs IDUwLCA3MCwgMTAwXSwgYW5kIGEgY2F0ZWdvcmljYWwgZ3JvdW5kIHRydXRoIHdoaWNoIGhhcyB0 d28gY2xhc3NlcyAoMCBhbmQgMSkgZS5nLiBbMCwgMCwgMCwgMCwgMSwgMSwgMV0uIEkgd2FudCB0 byBmaW5kIGEgdGhyZXNob2xkIHRoYXQgZGVjaWRlcyB3aGVuIGFuIGFuZ2xlIGJlbG9uZ3MgdG8g dGhlIGNsYXNzIDAgb3IgMS4NCj4gDQo+IEkgd2FzIGRvaW5nIHRoZSBmb2xsb3dpbmc6DQo+IA0K PiBmcm9tIHNjaXB5Lm9wdGltaXplIGltcG9ydCBjdXJ2ZV9maXQNCj4gDQo+IGRlZiBvYmplY3Rp dmUoeCwgdGhyZXNob2xkKToNCj4gIGlmIHggPiB0aHJlc2hvbGQ6DQo+ICAgIHJldHVybiAxDQo+ ICBlbHNlOg0KPiAgICByZXR1cm4gMA0KPiANCj4gIyBjdXJ2ZSBmaXQNCj4gdGhyZXNob2xkLCBf ID0gY3VydmVfZml0KG9iamVjdGl2ZSwgeCwgeSkNCj4gDQo+IGhvd2V2ZXIsIHRoaXMgY29kZSBk aWRuJ3Qgd29yayBhbmQgSSBnb3QgdGhpcyBlcnJvcjoNCj4gDQo+IFZhbHVlRXJyb3I6IFRoZSB0 cnV0aCB2YWx1ZSBvZiBhbiBhcnJheSB3aXRoIG1vcmUgdGhhbiBvbmUgZWxlbWVudCBpcyBhbWJp Z3VvdXMuIFVzZSBhLmFueSgpIG9yIGEuYWxsKCkNCj4gDQo+IGl0IHNlZW1zIHRoYXQgdGhlIHBy b2JsZW0gd2FzIGluIGhhdmluZyB0aGUgaWYgc3RhdGVtZW50IGluIHRoZSBvYmplY3RpdmUgZnVu Y3Rpb24sIGhvdyBjb3VsZCBJIHVzZSBvcHRpbWlzZSB0byBnZXQgdGhlIGJlc3QgdGhyZXNob2xk IGluIHN1Y2ggYSBjYXNlPw0KPiANCj4gVGhhbmtzDQo+IF9fX19fX19fX19fX19fX19fX19fX19f X19fX19fX19fX19fX19fX19fX19fX19fDQo+IFNjaVB5LVVzZXIgbWFpbGluZyBsaXN0IC0tIHNj aXB5LXVzZXJAcHl0aG9uLm9yZw0KPiBUbyB1bnN1YnNjcmliZSBzZW5kIGFuIGVtYWlsIHRvIHNj aXB5LXVzZXItbGVhdmVAcHl0aG9uLm9yZw0KPiBodHRwczovL21haWwucHl0aG9uLm9yZy9tYWls bWFuMy9saXN0cy9zY2lweS11c2VyLnB5dGhvbi5vcmcvDQo+IE1lbWJlciBhZGRyZXNzOiB0cmF2 aXNAdmF1Z2h0Lm5ldA0KDQoNClVubGVzcyB0aGlzIGlzIGp1c3QgYSBob21ld29yayBwcm9ibGVt LCBJ4oCZbSBub3Qgc3VyZSB3aHkgeW91IHdvdWxkIG5lZWQgY3VydmVfZml0IGZvciBhIGNhdGVn b3JpY2FsIChjbGFzc2lmaWNhdGlvbikgcHJvYmxlbeKAlGVzcGVjaWFsbHkgYSBiaW5hcnkgb25l LiAgV291bGQgaXQgYmUgcG9zc2libGUgdG8ganVzdCBkbyBzb21ldGhpbmcgbGlrZSB0aGlzPw0K DQoNCmltcG9ydCBudW1weSBhcyBucA0KDQp4ID0gbnAuYXJyYXkoWzIwLCAzMCwgNDgsIDQ2LCA1 MCwgNzAsIDEwMF0pDQp5ID0gbnAuYXJyYXkoWzAsIDAsIDEsIDAsIDEsIDEsIDFdKSAjIGV4YW1w bGUgY2hhbmdlZCB0byB0ZXN0IHVuc29ydGVkIGRhdGENCg0KeHkgPSBucC5hcnJheShbeCx5XSkN Cnh5LnNvcnQoKQ0KeDEgPSB4eVswXQ0KeTEgPSBucC5pbnZlcnQoeHlbMV0uYXN0eXBlKGJvb2wp KQ0KDQp0ID0geDFbeTFdWy0xXSAjIGZpbHRlciB0aGUgc29ydGVkIGFuZ2xlcyB3aXRoIHRoZSBz b3J0ZWQgY2F0ZWdvcmllcyBhbmQgZ2V0IHRoZSBtYXggKGp1c3QgcmV0dXJuIGluZGV4IFstMV0g c2luY2UgaXTigJlzIHNvcnRlZCkNCnByaW50KHQpDQoNCjQ2DQoNCg0KX19fX19fX19fX19fX19f X19fX19fX19fX19fX19fX19fX19fX19fX19fX19fX18KU2NpUHktVXNlciBtYWlsaW5nIGxpc3Qg LS0gc2NpcHktdXNlckBweXRob24ub3JnClRvIHVuc3Vic2NyaWJlIHNlbmQgYW4gZW1haWwgdG8g c2NpcHktdXNlci1sZWF2ZUBweXRob24ub3JnCmh0dHBzOi8vbWFpbC5weXRob24ub3JnL21haWxt YW4zL2xpc3RzL3NjaXB5LXVzZXIucHl0aG9uLm9yZy8KTWVtYmVyIGFkZHJlc3M6IGdjcHN1LXNj aXB5LXVzZXJAbS5nbWFuZS1teC5vcmcK