13 February 2009

While autotooling a project that uses GSL, I wanted GSL to first attempt to detect and use a system BLAS. If none can be found, then I wanted GSL to fall back to its own gslcblas implementation. The rational being that a savvy user could eke out some performance, but that for everyone else the code just builds and runs without being too fussy.

which a quick Google tells me is because -lgfortran isn't present in the linking line. I tried to obtain the required gfortran libraries in ${FLIBS} by using AC_FC_LIBRARY_LDFLAGS and then use them via

export GSL_CBLAS_LIB="${BLAS_LIBS} ${FLIBS}"

but later in the build I started hitting up against duplicate main issues because -lgfortranbegin winds up in ${FLIBS}.

More random Googling found the answer buried in lapack++'s configure script. The total solution looks like:

09 February 2009

The Austin Chapter of SIAM is proud to host the 3rd annual Texas Applied
Mathematics Meeting for Students (TAMMS) on March 27–28th. Attendees
will have an opportunity both to present their own research and to meet
fellow students from other Texas institutions. Though targeted at SIAM
members, this meeting is open to anyone with a mathematical leaning.
Please see http://www.ices.utexas.edu/siam
for more details.

Compared to the last time I wrote about the site, I've gotten the back button and search engine problems fixed. Or, I've got the back button working fine on Firefox. I'm unsure if dojo.back can be made to work on IE. Though not perfect, the site is usable sans JavaScript. I spent some time looking at a reCaptcha/PHP-based contact form but opted to forgo the contact form in the end.

It took me longer than I care to admit to track down and work through all the details/assumptions underneath the derivation of the conservative, compressible Navier-Stokes model for a perfect gas. This is on page two of every thesis I read, but is apparently such common knowledge that no one bothers to stick it in an appendix or give a good reference.

04 February 2009

For a turbulence warm up assignment, we needed to compute some correlation coefficients for trajectory behavior on the Lorenz attractor.
Specifically, compute the mean, variance, and cross-correlations of the three solution components. I'd skimmed through the Boost Accumulator Framework a couple of months back and thought it looked interesting:

Compared with the alternative loops/bookkeeping/etc., this is a pretty slick alternative. The only downside was that the accumulator declarations took me awhile to decipher from the documentation.

As you probably noticed, I am reading/writing plain text data files and computing the statistics at every time step. This is overkill and slow. I'd probably have fixed that by now if it weren't for the fact that I ran across Pipe Viewer in the last two days. Quite a neat shell tool that.