recursion - how to make double for loops faster in R -


i trying below calculation using r. function recursive , uses double loop calculate values of "result" matrix. there method replace loops or achieve if condition faster? 

x<-rnorm(2400,0, 3) y<-rnorm(400,0,3) no_row<-length(x) no_col<-length(y) input<-matrix(data=1,nrow = no_row, ncol = no_col) result<-matrix(nrow = no_row, ncol = no_col) calculation<-function(x,y) {   for(i in 1:no_row)  {   for(j in 1:no_col)   {     z<-exp(x[i]-y[j])     result[i,j]<-(z/1+z)    }  }   new_x<-x-1  new_y<-y-1  residual<-input-result      sq_sum_residulas<-sum((rowsums(residual, na.rm = t))^2)  if(sq_sum_residulas>=1){calculation(new_x,new_y)}   else(return(residual)) } output<-calculation(x,y)

the outer function tool looking for.

compare these 2 functions generate result matrix

x<-rnorm(100,0, 3) y<-rnorm(100,0,3)  calculation<-function(x,y) {     result <- matrix(nrow = length(x), ncol = length(y))   for(i in seq_along(x))   {     for(j in seq_along(y))     {       z<-exp(x[i]-y[j])       result[i,j]<-(z/1+z)     }   }    result }   calculation2 <- function(x, y){   result <- outer(x, y, function(x, y) { z <- exp(x - y); z / 1 + z})   result }  library(microbenchmark) microbenchmark(   calculation(x, y),   calculation2(x, y) )  unit: microseconds                expr     min       lq      mean   median        uq      max neval   calculation(x, y) 1862.40 1868.119 1941.5523 1871.490 1876.1825 8375.666   100  calculation2(x, y)  466.26  469.192  515.3696  471.392  480.9225 4481.371   100

that discrepancy in time seems grow length of vectors increases.

note, solve speed double for loop, there seem other issues in function. isn't clear me trying do, or why calling calculation within itself. have written, there no changes x , y before gets calling again, stuck in loop forever, if worked @ (it doesn't on machine)


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