scipy - Take Values inside of odeint in python -


my question if there's way take values in function not integrated in odeint.

exemple: if have derivative dy(x)/dt = a*x+ln(x) , before equation computed throught of intermediate equation a = b*d . take a's value during process.

more detailed (only exemple):

def func(y,t)     k = y[0]     b = 3     = cos(t**2) + b     dy/dt = a*t+ln(t)     return [dy/dt]

can take a's values of function?

the answer josh karpel

the code that:

def reaction(state,t):    # integrate results    p = state[0]    t = state[1]    # function determine enthalpy of system    f1(t,p) = enthalpy    # function determine specific volume of system    f2(t,p) = specific volume    # function determine heat release reactions    f3(t,p,t) = heat release reactions

# derivatives

 dp/dt = f(t,p,enthalpy,specific volume,heat release reactions)   dt/dt = f(t,p,enthalpy,specific volume,heat release reactions)

the real code bigger that. but, know if there way store values of f1 (enthalpy), f2 (specific volume), f3 (heat release) vector or tuple during process of solution of odeint same size of p , t.

it's not entirely clear want, sounds need pass value function you're integrating over. there 2 options can think of:

  1. scipy.integrate.odeint takes args argument contains arguments passed integrand function, have signature y(t, a).
  2. you use functools.partial construct new function has argument a integrand function y(t, a) set.

Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj...
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer ...
Read more

c# - Asp.net web api : redirect unauthorized requst to forbidden page -

im trying redirect unauthorized request forbidden page instead i'm getting forbidden page in response body , how can fix ? here's startup class : app.createperowincontext(storecontext.create); app.createperowincontext<applicationusermanager>(applicationusermanager.create); app.createperowincontext<applicationsigninmanager>(applicationsigninmanager.create); app.usecookieauthentication(new cookieauthenticationoptions { authenticationtype = defaultauthenticationtypes.applicationcookie, expiretimespan = timespan.fromdays(30), }); app.useexternalsignincookie(defaultauthenticationtypes.externalcookie); app.usetwofactorsignincookie(defaultauthenticationtypes.twofactorcookie, timespan.fromminutes(5)); app.usetwofactorrememberbrowsercookie(defaultauthenticationtypes.twofactorrememberbrowsercookie); the method im trying reach : [httpget] [authorize(roles = "admin")] public string getcurrentusername() { return userman...
Read more