Thanks all. These solutions look great. Lots of things to learn from. I consider myself a master Jedi in the “art” of MATLAB (I am the developer of the GIBBON project), but I am slowly realizing that converting super optimized MATLAB code nearly directly into Julia code quite often does not equal super optimized Julia code… at all. Then by optimizing the Julia code, and embracing the “Julian approach” it ends up faster than I thought was possible. Sadly for me though none of the structures/approaches/syntax to speed things up come naturally to me yet, leaving me feeling demoted to a Padawan learner. In MATLAB one avoids for-loops like the plague, with Julia one seems to have to embrace them like the cure. However the promise of unrivaled speed, and a fully open source software project, keep me motivated 
.
2 Likes