---Removed all loops over theta and picked a specific theta value to
calculate at. This gives significant speedup and reduces memory
footprint. Time correlation still needs work. Look at page 385 and eqn
11.101 of B&P.