This is C++ code utilzing GSL routines to perform bicubic interpolation on tabular data for the purpose of accurately evaluating micropysical source terms to the equations of gas dynamics. It was specifically developed to utilize the tables here.

The code is self-contained in a .h file so that it can easily be called within a hydro code. For example, to interface with Athena++, create folder named 'user' in /src and place vegas_tables.h there, and then call it from any pgen using

#include "../user/vegas_tables.h"  

The example driver bicubic_speed_test.cpp provides an example of initializing and using a table. It compares the speed of 1000 random evalutions using both bilinear and bicubic interpolation on the sample table Blondin_100by200.dat. Compile and run it using

g++ bicubic_speed_test.cpp -o place_bets -lgsl -lm
./place_bets

For further documentation, refer to the header information in vegas_tables.h

There is a C-wrapper for this code available. Just contact me to get a hold of that.