import numpy as np

# Contains All Two-Body Matrix Routines

#----------------------------------------------------------------------------
def propagate_two_body(wf_up,wf_down,size,n_up,n_down,U=0.0, delta_tau=0.1):
   # Propagates By Two Body Terms As Separated by HS Transform

   # Form One Body Constant In Front of HS Transformation
  
   # Form Constant That Goes Into Two-Body Exponential

   # Run Through Sites Sampling Fields and Forming One Body HS Terms*/
   #for sites in range(size):

    # Generate Gaussian-Distributed Random Field   

    # Store Exponential Two Body Terms

    # Multiply Constants Into Wave Functions

   pass
#----------------------------------------------------------------------------

