d3q Code

Examples updated

Example updated

Removed twfcollect and lkpoint_dir

New checkpoint because I rolled back to the wrong version

Making a check point before rolling back to previous version, because there is something not working properly

Also note that everything can be done 2x faster, but using much more RAM if we use the fact that A_in is symmetric (implemented here, but not in parallel)

Looks like my brilliant idea was not that brillians, fixing previous commit, now two matrices are used again (but not for SMA) also looks stable but I have to test

New form for variational tk (as in PRB paper) seems to work much better than the previous one (that was from the thermalk code), it is also twice as fast as the same matrix can be reused, but some more testing is needed. In particular, I have some doubt about the index to use in freqm1

recompute-sma modified to use only files from tk.x with store_lw=.true.: using hand-made lineiwidths is too complicated and you'ld have to do it by hand anyway
also, recompute-sma-mfp shows how to recompute sma with a cut-off on the mean free path

I did not put the 0.5 factor after all..

Implemented a smart mechanism to stop the code immediately if the step end is after the wall limit from input, to save CPU time