Implemented analytic integration.
Including Aaron's spin-boson code for comparison.
Set up Aaron's FSSH code and mine for comparison.
Aaron's code seems to start on the excited state?   I think I might be missing something here.
I get idential bath dynamics on the excited state
About to try changing the integrator.
About to modify bath so that v*d = v*d12 always
Analytic wavefunction integration now agrees.
Added hopping prob. according to Aaron's scheme.
Where do small differences in p dot d come from?
Did I change something I shouldn't have?
About to remove a bunch of print statements
Bath energies and A11s agree *exactly*!!
A22 also agrees exactly
V dot d is identical as well
Off-diagonal elements of the density disagree!
Found a compiler specific bug when using real() -- specifically for arguments of type double complex.  It only happens when using default gfortran.
Analzying Aaron's code introduced a bug into it.  But now our hopping probabilities agree!  Basically, I commented out something that I didn't think should have make a difference, but it really did.  Now that's fixed.
Updated bath: it now computes initial state.
Added ability to get diabatic wavefun. coefficient