Skip to content

Commit

Permalink
Remove commented out code
Browse files Browse the repository at this point in the history
  • Loading branch information
dschwoerer committed Nov 8, 2024
1 parent 925e0d6 commit 23e72bd
Showing 1 changed file with 1 addition and 17 deletions.
18 changes: 1 addition & 17 deletions zoidberg/stencil_dagp_fv.py
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down Expand Up @@ -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
# Λ
Expand Down Expand Up @@ -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
Expand All @@ -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:])

Expand Down Expand Up @@ -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)

0 comments on commit 23e72bd

Please sign in to comment.