From 23e72bd2b2c563c2c0590fa9f5a6be5aab8065d2 Mon Sep 17 00:00:00 2001 From: David Bold Date: Fri, 8 Nov 2024 10:16:20 +0100 Subject: [PATCH] Remove commented out code --- zoidberg/stencil_dagp_fv.py | 18 +----------------- 1 file changed, 1 insertion(+), 17 deletions(-) diff --git a/zoidberg/stencil_dagp_fv.py b/zoidberg/stencil_dagp_fv.py index f2cb219..cbdb183 100755 --- a/zoidberg/stencil_dagp_fv.py +++ b/zoidberg/stencil_dagp_fv.py @@ -61,12 +61,6 @@ def load(fn): def doit(pols, plot=False): - # with DF(fn) as f: - # RZ = [f[k] for k in "RZ"] - # log("opened") - # RZ = np.array(RZ) - # log("read") - RZs = [] ### Calculate Volume of the cell @@ -188,8 +182,6 @@ def doit(pols, plot=False): coefsX = np.linalg.solve(dxzR, dRr.transpose(1, 2, 3, 0)) log("done") - # In[154]: - # AreaZplus # Z # Λ @@ -238,12 +230,11 @@ def doit(pols, plot=False): test(RZ, volume, coefsX, coefsZ, plot=plot) return write(RZ, volume, coefsX, coefsZ) - # In[155]: def test(RZ, volume, coefsX, coefsZ, plot=False): inp = np.sin(RZ[1]) - ana = -np.sin(RZ[1]) # + np.cos(RZ[0])/RZ[0] + ana = -np.sin(RZ[1]) if 0: inp = RZ[1] ** 2 ana = RZ[1] ** 0 @@ -254,12 +245,6 @@ def test(RZ, volume, coefsX, coefsZ, plot=False): inp = RZ[0] ** 3 ana = 6 * np.sin(RZ[0]) - # print("Fuck it up!!!") - # coefsX[..., 1] *= -1 - # coefsZ[..., 0] *= -1 - - # In[156]: - def xp(ijk): return (ijk[0] + 1, *ijk[1:]) @@ -373,5 +358,4 @@ def write(RZ, volume, coefsX, coefsZ): for fn in sys.argv[1:]: print(fn) - # doit(fn) test(*load(fn), plot=plot)