wassname/simpeg
1
0
Fork 0
mirror of https://github.com/wassname/simpeg.git synced 2026-06-28 22:21:15 +08:00
Code Issues Packages Projects Releases Wiki Activity

Updates to merging

Browse Source
This commit is contained in:
Rowan Cockett
2015-11-15 17:51:51 -08:00
parent 5e20960335
commit 4ac85e81d6
1 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files
SimPEG/Mesh/TreeMesh.py
+10 -10
View File
@@ -1864,24 +1864,24 @@ class NotBalancedException(Exception):
if __name__ == '__main__':
def function(xc):
r = xc - np.array([0.5*128]*len(xc))
def function(cell):
r = cell.center - np.array([0.5]*len(cell.center))
dist = np.sqrt(r.dot(r))
# if dist < 0.05:
# return 5
if dist < 0.1*128:
if dist < 0.1:
return 4
if dist < 0.3*128:
if dist < 0.3:
return 3
if dist < 1.0*128:
if dist < 1.0:
return 2
else:
return 0
# T = Tree([[(1,128)],[(1,128)],[(1,128)]],levels=7)
# T = Tree([128,128,128],levels=7)
T = Tree([[(1,16)],[(1,16)]],levels=4)
# T = Tree([[(1,128)],[(1,128)]],levels=7)
# T = TreeMesh([[(1,128)],[(1,128)],[(1,128)]],levels=7)
# T = TreeMesh([128,128,128],levels=7)
T = TreeMesh([16,16],levels=4)
# T = TreeMesh([[(1,128)],[(1,128)]],levels=7)
# T.refine(lambda xc:1, balance=False)
# T._index([0,0,0])
# T._pointer(0)
@@ -1894,7 +1894,7 @@ if __name__ == '__main__':
T.plotImage(np.random.rand(T.nC),showIt=True)
print T.getFaceInnerProduct()
# print T.getFaceInnerProduct()
# print T.gridFz