python - Symbolically computing the input of an interpolation function? -


i have rather complicated function h(x), , i'm trying solve value of x such h(x) = constant. interpolation object generated discrete interval , corresponding output of h(interval), other inputs held constant. denote interpolation object f.

my problem call function of interpolation object accepts array_like, passing symbol f(x) use sage's solver method out of question. ideas of how around this?

i have interpolation function f. solve equation f(x) == sageconstant forx.

    scipy.interpolate import interpolatedunivariatespline iuspline     import numpy np      #generating interpolation object     xint = srange(30,200,step=.1)     val = [h(i,1,.1,0,.2,.005,40) in srange(30,299,step=.1)]     f = iuspline(xint,val,k=4)      #this yield sage constant     eq_g(x) = freeb - x       #relation solve     eq_m(x) = eq_g(39.9) == f(x)     m = solve(eq_m(x),x)

the above code (f(x) more specific) generates

"typeerror: cannot cast array data dtype('0') dtype('float64') according rule 'safe'.

edit: function h(x) result in same error, hence doesn't matter h(x) is. simplicity (i wasn't kidding when said h complicated), try h(x) = x. block read:

    scipy.interpolate import interpolatedunivariatespline iuspline     import numpy np      #generating interpolation object     xint = srange(30,200,step=.1)     h(x) = x     val = [h(i) in srange(30,299,step=.1)]     f = iuspline(xint,val,k=4)      #this yield sage constant     eq_g(x) = freeb - x       #relation solve     eq_m(x) = eq_g(39.9) == f(x)     m = solve(eq_m(x),x)

when working numpy , scipy, prefer python types sage types.

instead of sage integers , reals, use python ints , floats.

maybe can fix code this.

from scipy.interpolate import interpolatedunivariatespline iuspline import numpy np  # generate interpolation object xint = srange(30,200,step=.1) xint = [float(x) x in xint] val = [float(h(i,1,.1,0,.2,.005,40)) in srange(30,299,step=.1)] f = iuspline(xint,val,k=4)  # yield sage constant eq_g(x) = freeb - x   # relation solve eq_m(x) = eq_g(39.9) == f(x) m = solve(eq_m(x),x)

Comments

Post a Comment

Popular posts from this blog

java - Spring Data JPA: Why findOne(id) executing delete query internally? -

i have 1 application using @manytomany annotation link 2 different entities 1 parent entity. not sure if have implemented properly. when execute .findone() on parent entity's repository find it, @ first fetch data along executing delete queries internally delete record mapping table (a table jpa creates after adding @manytomany annotation.). following generated sql queries , code snippet: java class @entity @table(name = "inventory") public class inventory { @id long id; // many other properties @manytomany(fetch=fetchtype.lazy) private list<stream> streams; @manytomany(fetch=fetchtype.lazy) private list<subject> subjects; } there no other annotation i've put in stream , subject entity class link inventory class. when run app, create following tables: inventory stream subject inventory_stream inventory_subject (inventory_stream , inventory_subject, 2 tables created automatically, don't k...
Read more

python - Mongodb How to add addtional information when aggregating? -

i beginner , wrote line pipeline works, wanna add other information output, screen name , or number of tweets.i tried add under $group gave me syntax error everytime here pipeline: def make_pipeline(): # complete aggregation pipeline pipeline = [ { '$match': { "user.statuses_count": {"$gt":99 }, "user.time_zone": "brasilia" } }, { "$group": { "_id": "$user.id", "followers": { "$max": "$user.followers_count" } } }, { "$sort": { "followers": -1 } }, { "$limit" : 1 } ]; i using on example : { "_id" : objectid("5304e2e3cc9e684aa98bef97"), "text" : "first week of school on :p", "in_reply_to_status_id" : null, "retweet_count" : null, ...
Read more

java - Incorrect order of records in M-M relationship in hibernate -

Image
i have 4 records in bridge table in mysql database. using @manytomany in hibernate, accurate order of employees incorrect incorrect of projects. following required material figure out doing mistake. this source. @entity @table(name="project") public class projectbean { @id @generatedvalue @column(name="project_id") private int projectid; @column(name="title") private string projecttitle; @manytomany(mappedby="projects") private collection<employeebean> employees; public projectbean() { employees = new arraylist<employeebean>(); } //getters & setters } @entity @table(name="employee") public class employeebean { @id @generatedvalue @column(name="employee_id") private int employee_id; @column(name="first_name") private string ...
Read more