Numerical Recipes 3rd Edition: The Art of Scientific Computing

Category: Software
Author: William H. Press
All Stack Overflow 7
This Month Stack Overflow 3


by anonymous   2019-07-21

Have a look at Boost Graph as it will simplify all your graph work. I am unsure if there is a Markov Chain algiorithm (I am looking for one too) but it should be easy to write once you have the graph -- a concurrent monte carlo approach maybe?

Numerical Recipes has many algorithm and code in both C and C++.

The graphviz tools are all you need to draw graphs.

by anonymous   2019-07-21

As described in Numerical Recipes, the determinant is simply the product of the diagonal elements after an LU decomposition. (This book is available online, and has code for computing the LU decomposition).

by anonymous   2017-08-20

I think it's a very bad idea to calculate two very large numbers and hope that the quotient comes out as something sensible.

Start by taking the natural log:

ln(n!/(q!)^r) = ln(n!) - r*ln(q!)

You can use gammaln() for the two function values, simplify, then take exp() to get the result you want:

value = exp(gammln(n+1) -r*gammln(q+1))

Numerical Recipes has a nice chapter on how to implement functions like gammln().

by anonymous   2017-08-20

On many platforms (such as any modern x86-compatible), many of the maths functions are implemented directly in the floating-point hardware (see for example Not all of them are used, though (as I learnt through comments to other answers here). But, for instance, the sqrt library function is often implemented directly in terms of the SSE hardware instruction.

For some ideas on how the underlying algorithms work, you might try reading Numerical Recipes, which is available as PDFs somewhere.