L-Nafaryus/anisotropy
L-Nafaryus/anisotropy
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fix: normal vector

Browse Source
This commit is contained in:
L-Nafaryus 2021-04-02 12:22:55 +05:00
parent cc3e03ab27
commit 1fe34c9c48
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
samples/__init__.py
View File

@ -33,7 +33,7 @@ def genMesh(stype, theta, flowdirection, saveto):
if flowdirection == [1, 1, 1]: if flowdirection == [1, 1, 1]:
geometry = geometry2 geometry = geometry2
norm = [-1, -1, 1] norm = [-1, 1, 0]
bcount = 6 bcount = 6
# initial angle # initial angle
@ -46,7 +46,8 @@ def genMesh(stype, theta, flowdirection, saveto):
# math.sin(angle * 0.5) * flowdirection[2]) # math.sin(angle * 0.5) * flowdirection[2])
#qvec = (ax * vec * ax.inverse).vector #qvec = (ax * vec * ax.inverse).vector
#normvec = [qvec[0], qvec[1], qvec[2]] #normvec = [qvec[0], qvec[1], qvec[2]]
normvec = Quaternion(axis = flowdirection, angle = angle).rotate(norm) v1 = Quaternion(axis = norm, angle = math.pi / 2).rotate(flowdirection)
normvec = Quaternion(axis = flowdirection, angle = angle).rotate(v1)
direction = [1, 1, 1] direction = [1, 1, 1]
#direction = fd([1, 1, 1], [1, -1, 1], [1, -1, -1]) #direction = fd([1, 1, 1], [1, -1, 1], [1, -1, -1])