Description
Currently, we pass the omega and q values as a tuple to all sorts of functions, e.g., the dynamic structure factor. This is very inconvenient, when it comes to broadcast over vectors of q and omega. Therefore, we should replace the tuple input just by two arguments, omega and q.
Validation and testing
ne = 1.0e21u"cm^(-3)"
temp = 10.0u"keV"
om_arr = range(0.1,1.0,100)
q = 0.1
esys = IdealElectronSystem(ne,temp)
old_res = [dynamic_structure_factor(esys,(om,q)) for om in om_arr]
new_res = dynamic_structure_factor.(esys,om_arr,q)
@assert old_res == new_res
Description
Currently, we pass the omega and q values as a tuple to all sorts of functions, e.g., the dynamic structure factor. This is very inconvenient, when it comes to broadcast over vectors of q and omega. Therefore, we should replace the tuple input just by two arguments,
omegaandq.Validation and testing