Hi Ivar, yeah, I know, I thought Andreas was replying to my deleted post. I think Petr already solved the problem with the globals (his gist was apparently not the right context). However, he still reported:
On Thursday, December 11, 2014 6:47:33 PM UTC+1, Petr Krysl wrote: > > One more note: I conjectured that perhaps the compiler was not able to > infer correctly the type of the matrices, so I hardwired (in the actual FE > code) > > Jac = 1.0; gradN = gradNparams[j]/(J); # get rid of Rm for the moment > > About 10% less memory used, runtime about the same. So, no effect really. > Loops are still slower than the vectorized code by a factor of two. > > Petr > > Best, Robert On Friday, December 12, 2014 12:01:44 AM UTC+1, Ivar Nesje wrote: > > https://gist.github.com/anonymous/4ec426096c02faa4354d#comment-1354636