Notes made for separating 'Co' and 'lamc' from eigenvalues.  Expect three changes in 'KLresearch' when creating Eigenvalue values, 2 in 'KLmeanadjust' when using them (one evaluating at a point, one evaluating an integral), and 1 in 'radtransMC' when evaluating at a point, duplicate function of one in 'KLmeanadjust', someday should make only one.  Also fixed a spelling error in a comment and added to 'codingnotes'.
Made actual 6 changes (plus 'use'ing variables) to pull 'CoExp' and 'lamc' out of Eigenvalues.  Previous commit flags where changes are to be made, they are actually made here.  Transport values seem as expected.  In reactor like case here, LP and atomic mix actually do very well, although you can see convergence upon correct values going from 1 to 2 to 3 to 4 eigenvalues using this method.
Added comments for upcoming addition of the new material centric xi distributions, which will allow noncongruent input scattering ratios for the two materials in a binary Markov case.  I updated the comments for KL_collect, and commented where changes will be made.  Specifically with that: 'KLresearch', 'KL_collect' update this sub, 'radtransMC' 'MCtransport' use new function for evaluation of scattering ratio and add new function 'KLWood_actscatrat' to module, 'variables' 'KLvars' add 'Coscat' and 'Coabs', 'Loadcase' 'testinputstoc' allow dissimilar scattering ratios.
Added flag 'KLxigentype' with values of either 'totxs' or 'material' to differentiat old implementation and new implementation here.  Also modified and verified test in 'testinputstoc' to allow dissimilar cross sections when in 'material' mode, but not 'totxs' mode.
Added new variables in 'variables'.  Used flag to implement different generations of xi from common coding depending on input.  Old method still produces same results, new method runs, and should get correct results after implementing rest of method.
'Save commit', not verified at all.  I have added 'KLWood_actscatrat' function, expanded 'KLrxi_point2' function which the previous calls, added precomputing of 'Coscat' and 'Coabs'.  I need to run verification tests on all of this work, and consider precomputing the average scattering xs and average absorption xs now used in 'KLrxi_point2'.
I think new xi gen scheme implemented correctly.  But I have found that I did not remove Co and lamc from their places correctly.  My work with that mainained the mean, but reduced the effective variance so much that the problem is nearly atomic mix.  Must fix that, then finish verifying this.
Sidebar: moved 'Coscat' and 'Coabs' to be 'genRealzvars' like they should be.  Created 'sigscatave' and 'sigabsave' as 'genRealzvars', and precomputed them.  Results are as before this commit, with an error in the way I pulled out 'Co' and 'lamc' remaining.
Sidebar: moved 'sigave' from 'KLvars' to 'genRealzvars'.  Everything else is as before.
Tried to make Co pulling out work.  Not there yet.  Variance of reconstructed cross sections always greatly reduces, which should not happen.  If and when I figure it out, it would be nice to get rid of this commit.  I really should learn to stash.
Sidebar: Added 'compare.sh'.  It is a bash script which will run the currenty copy of the code side-by-side with the same code which is located at '../stochastic2/stochastic' relative to this copy.  For outputs to compare, include the string 'flag' in print statements.  This is to be used for side-by-side debugging.
Progress: Finished pulling the Co term out right.  Actually, since there is 'sqrt(Eig)' in the denominator of the xi values, you don't need one in front of the summation.  You can pull anything out you want without affecting the reconstructions, since the eigenvalue cancels itself out anyway.
Added option in 'KLrxi_point2' to specify number of eigenvalues to use.
Sidebar: Made input to 'KLrxi_point2' optional which specifies whether to create total, scattering, or absorption cross section.
Sidebar: Created 'Eigenvalue', a function to solve for values of eigenvalues.  Currently there are three times I solve this, common coding is good.  Although now that I have taken variance out of the eigenvalue solve I can probably reduce this to one instance.
Sidebar: Replaced all instances of 'KLrxi_point' with 'KLrxi_point2' and got rid of 'KLrxi_point'.
Sidebar: Changed all instances of 'KLrxi_point2' to 'KLrxi_point'.
Sidebar: Moved 'KLrxi_point' from 'radtransMC' to 'KLreconstruct' where it belongs.
Sidebar: Moved 'Eigfunc' from 'KLmeanadjust' to 'KLreconstruct' where it belongs.
Sidebar: Moved 'CoExp' from 'KLvars' to 'genRealzvars' where it belongs.