See Supplementary material
http://link.springer.com/article/10.1007%2Fs10950-014-9475-7