Create photon_scat_gauss_2.c

Exceedingly fast, gcc version of theCashwell-Carter / Prahl combined with Kirk ATTEN5 program in c
	 @ Copyright (c), 2014 
         Donald L Walters Sept 17, 2014, W9DKI
         Hilbert Space Technical Arts, LLC
         Monterey, California, 93940

	 Has been used to follow 10^12 individual photons in the turbid ocean/atomosphere 
	 and stores their trajectory lengths and exit angles efficiently.   
	 Based on the Cashwell-Carter LLNL neutron scattering code highly modified and adapted
	 
	 ***
         CAUTION
	 This is intended to be used with a water cooled, i7 to avoid CPU over damage. 
	 Use at your own risk otherwise! You may very well fry you CPU since it
         loads 1-2 threads at 100% for up to 30 minutes depending on the # of photons 
          photons = 100000000;  (default)  // one billion photons tracked typically runs in ~1 minute
          on a 2600K i7 ,turbo extreme mode  in a ASUS P8P67 with intel BXTS13X water cooler
          with a Titan K-20 GPU
	 ****
  Simulates a Gaussian laser beam initial size wo and includes diffraction in water or air
	 SSE2 Meresenne twister random number generator
	 and the very fast, -pi to pi sine and cosine polonomial algorithm
	 by Charles K Garrett, Feb 2012

   standard arrays with option base 0
   calloc is needed for arrays > 100000
   Petzold data must not be CSV

    compile with:
    gcc -Ofast -msse2 -frename-registers -fno-strict-aliasing -DHAVE_SSE2=1 -DDSFMT_MEXP=19937 -o photon_scat_gauss dSFMT.c photon_scat_gauss.c -lrt -lm -lgsl -lgslcblas
    
 Uses: dsfmt.c and dsfmt.h frp, Saito and Matsumoto  
 contained in dsfmt-ser-2.2.3
    Copyright (C) 2007, 2008 Mutsuo Saito, Makoto Matsumoto and
    Hiroshima University. All rights reserved.
    Copyright (C) 2012 Mutsuo Saito, Makoto Matsumoto,
    Hiroshima University and The University of Tokyo.
    All rights reserved. see  BSD License in dsfmt.h