def getdMu(h5pFile): TINY = 1.0e-8 pIds, Muf = lfmpp.getH5pFin(h5pFile, "Mu") pIds, Mu0 = lfmpp.getH5pInit(h5pFile, "Mu") Ind = (Muf > TINY) & (Mu0 > TINY) Muf = Muf[Ind] Mu0 = Mu0[Ind] #dMu = np.abs(Muf-Mu0)/Mu0 dMu = (Muf - Mu0) / Mu0 print("\nCutting %d particles for undefined Mu" % (len(Ind) - Ind.sum())) print("Min dMu = %f" % (dMu.min())) print("Min Muf = %f, Min Mu0 = %f" % (Muf.min(), Mu0.min())) print("\n\n") return dMu, Ind
def getA0(h5pFile): pIds, A0 = lfmpp.getH5pInit(h5pFile, "alpha") return A0
def getdK(h5pFile): pIds, Kf = lfmpp.getH5pFin(h5pFile, "kev") pIds, K0 = lfmpp.getH5pInit(h5pFile, "kev") dK = Kf / K0 return dK