Archive for March, 2015

Tuesday, March 31st, 2015

A use for common blocks

I have long implored others to abandon the use of common blocks for a myriad of reasons. They create a web of interconnections between the routines that use the common block (and which routines, exactly, are those? You can find out by grepping, but it is not immediately obvious). It is not clear at all […]

No Comments » - Posted in IDL by

Tuesday, March 3rd, 2015

Batched matrix-vector multiplication

After spending awhile last Friday trying to vectorize a loop of a small matrix-vector multiplication for every pixel of an image, I gave up and decided to just write it as a DLM. For my image sizes of 1024 by 1024 pixels (actually two images of that size), the run time went from 3.15 seconds […]

No Comments » - Posted in HPC,IDL,mglib by