python - Constrained maximization of function and not-successful attempts in Python3 -


firstly, present problem. latex description of problem

my attempt

i tried solve problem scipy.optimize package. had pass initial vector (initial guess). , value, given algorithm dependent on initial guess. tried take random initial guesses , tried find maximal values among more , more initial guesses

import numpy np scipy.optimize import minimize  def j(data):     l = len(data)     if (l%2==0 | l<2):         raise exception('wrong length or arguments!')      phi = data[0:l//2]     p = data[l//2:l-1]     theta = data[l-1]      phi = np.array(phi)         p = np.array(p)      w = (1 + np.cos(phi)*np.cos(theta))*p     lam = (1+np.cos(phi-theta))/(2*(1+np.cos(phi)*np.cos(theta)))      return -(1 - sum(w * ((-lam*np.log2(lam)) - (1-lam)*np.log2(1-lam))))  def probability_constraint(x):     return sum(x[len(x)//2:len(x)-1])-1  def lambda_w_constraint(data):     l = len(data)     phi = data[0:l//2]     p = data[l//2:-1]     theta = data[l-1]      phi = np.array(phi)         p = np.array(p)      w = (1 + np.cos(phi)*np.cos(theta))*p     lam = (1+np.cos(phi-theta))/(2*(1+np.cos(phi)*np.cos(theta)))         return sum(w*lam)-1/2  def lambda_greater_than_zero(data):     l = len(data)     phi = data[0:l//2]     p = data[l//2:-1]     theta = data[l-1]      lam = (1+np.cos(phi-theta))/(2*(1+np.cos(phi)*np.cos(theta)))         return lam  def probability_positive(x):     return x  cons = ({'type': 'ineq', 'fun' : lambda x: lambda_greater_than_zero(x)},         {'type': 'eq', 'fun' : lambda x: probability_constraint(x)},         {'type': 'eq', 'fun' : lambda x: lambda_w_constraint(x)},         {'type':'ineq',  'fun' : lambda x: probability_positive(x)}  )  def maximize(numiter):     maxj = -1000     maxx = []      in range(0, numiter):         x0 = np.random.rand(5)*10         res = minimize(j, x0, constraints=cons)         r = j(res.x)         if (r>maxj):             maxj = r             maxx = res.x      return maxj, maxx  maximize(10)  maximize(100)

i know, method naive , incorrect.

so question is: know how solve problem? can in python3, other languages (matlab?) me.


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

jquery - Responsive Navbar with Sub Navbar -

good day, forgive me nav vocabulary not best. if question confusing apologize. i'm working on navigation bar needs centered on page @ 768px , above. items in nav open (or perhaps toggle) subnav lives directly below. these nav should appear centered on page. i have been able work through responsive portion of nav bar, having main items align vertically , sub nav menu of each item roll out below, , gave them indent separation. working well. however, issue not responsive side, behavior @ width greater 768px. happens toggling menu's sub navigation pushes other main items down page below sub nav. here's code: $(document).ready(function () { //stick in fixed 100% height behind navbar don't wrap $('#slide-nav.navbar-inverse').after($('<div class="inverse" id="navbar-height-col"></div>')); $('#slide-nav.navbar-default').after($('<div id="navbar-height-col"></d
Read more