-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmain.py
66 lines (49 loc) · 1.44 KB
/
main.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
import numpy as np
from meshHandle.multiscaleMesh import FineScaleMeshMS as msh
import time
import pdb
import geoUtil.geoTools as gtool
from math import pi, sqrt
from pymoab import core, types, rng, topo_util
# docker run -t -it -v /home/arturcastiel/projetos:/pytest desenvolvimento:latest bash -c "cd /pytest; bash"
# docker run -it -v /home/arturcastiel/projetos:/pytest desenvolvimento:latest bash -c "cd /pytest; bash"
# %load_ext autoreload
# %autoreload 2
start = time. time()
M = msh("malha03.msh", dim = 2)
vec = np.arange(len(M.alma)).astype(int)
# M.core.print()
end = time. time()
print("Execution time for a {2}d mesh with {0} elements: {1} seconds".format(len(M), end-start, M.dim ))
start = time. time()
#M.core.print()
end = time. time()
print("Printing Routine Time: {}".format(end-start))
start = time. time()
M.faces.center[:]
end = time. time()
print("Printing Routine Time: {}".format(end-start))
# start = time. time()
# M.faces.normal[:]
# end = time. time()
# print("Printing Routine Time: {}".format(end-start))
# vec2 = np.array([vec,vec,vec]).T
# #M.alma.set_data(vec2)
# M.alma[:] = vec2
# vec3 = np.array([10,12,13,14,30,35])
#
# vec1 = np.arange(8)
# p = vec2[2:10,:]
#
#
#
# M.core.print()
# cumaru = M.core.mb.tag_get_handle("GEOM_DIMENSION")
# M.core.handleDic["GEOM_DIMENSION"] = cumaru
#
#
#
#
# loook = M.core.mb.get_entities_by_type_and_tag(
# 0, types.MBENTITYSET, np.array(
# (cumaru,)), np.array((None,)))