nalc/vector_extras
1
0
Fork 0
forked from mtcontrib/vector_extras
Code Issues Pull Requests Releases Wiki Activity

edit code

Browse Source
This commit is contained in:
HybridDog 2015-04-25 08:35:23 +02:00
parent c75042050a
commit 02291fcb65
1 changed files with 12 additions and 55 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

67
init.lua
View File

@ -282,57 +282,14 @@ function vector.area(ps)
local vn = vector.normalize(vector.crossproduct(B, C)) local vn = vector.normalize(vector.crossproduct(B, C))
--[[bB.x*C.x+bB.y*C.y+C.z = 0 local nAB = vector.normalize(B)
bB.x*vn.x+bB.y*vn.y+vn.z = 0 local nAC = vector.normalize(C)
local angle_BAC = math.acos(vector.scalarproduct(nAB, nAC))
bB.x = -(bB.y*C.y+C.z)/*C.x local nBA = vector.multiply(nAB, -1)
bB.x = -(bB.y*vn.y+vn.z)/vn.x local nBC = vector.normalize(vector.subtract(C, B))
local angle_ABC = math.acos(vector.scalarproduct(nBA, nBC))
(bB.y*C.y+C.z)/*C.x = (bB.y*vn.y+vn.z)/vn.x
(bB.y*C.y+C.z)*vn.x = C.x*(bB.y*vn.y+vn.z)
bB.y*C.y*vn.x+C.z*vn.x = C.x*bB.y*vn.y+C.x*vn.z
bB.y*C.y*vn.x-C.x*bB.y*vn.y = C.x*vn.z-C.z*vn.x
bB.y*(C.y*vn.x-C.x*vn.y) = C.x*vn.z-C.z*vn.x
bB.y = (C.x*vn.z-C.z*vn.x)/(C.y*vn.x-C.x*vn.y)
bB.y = -(bB.x*C.x+C.z)/C.y
bB.y = -(bB.x*vn.x+vn.z)/vn.y]]
--[[
local v3 = vector.subtract(C, B)
local bB = {
x = (C.y*vn.z-C.z*vn.y)/(C.x*vn.y-C.y*vn.x),
y = (C.x*vn.z-C.z*vn.x)/(C.y*vn.x-C.x*vn.y),
z = 1,
}
local cC = {
x = (B.y*vn.z-B.z*vn.y)/(B.x*vn.y-B.y*vn.x),
y = (B.x*vn.z-B.z*vn.x)/(B.y*vn.x-B.x*vn.y),
z = 1,
}
local aA = {
x = (v3.y*vn.z-v3.z*vn.y)/(v3.x*vn.y-v3.y*vn.x),
y = (v3.x*vn.z-v3.z*vn.x)/(v3.y*vn.x-v3.x*vn.y),
z = 1,
}]]
local B2 = vector.mirror(B, vn, C)
local C2 = vector.mirror(C, vn, B)
local BC = vector.subtract(C, B)
local BA = vector.multiply(B, -1)
local A2 = vector.add(vector.mirror(BA, vn, BC), B)
local nB = vector.normalize(B)
local nC = vector.normalize(C)
local angle_BC = math.acos(vector.scalarproduct(nB, nC))
local BA = vector.normalize(vector.multiply(B, -1))
local BC = vector.normalize(vector.subtract(C, B))
local angle_AC = math.acos(vector.scalarproduct(BA, BC))
local bB = {z=1}
local area = {} local area = {}
for z = cube_p1.z, cube_p2.z do for z = cube_p1.z, cube_p2.z do
for y = cube_p1.y, cube_p2.y do for y = cube_p1.y, cube_p2.y do
@ -364,15 +321,15 @@ function vector.area(ps)
end end
end--]] end--]]
local nep = vector.normalize(ep) local nep = vector.normalize(ep)
local angle_Bep = math.acos(vector.scalarproduct(nB, nep)) local angle_BAep = math.acos(vector.scalarproduct(nAB, nep))
local angle_Cep = math.acos(vector.scalarproduct(nC, nep)) local angle_CAep = math.acos(vector.scalarproduct(nAC, nep))
local angldif = angle_BC - (angle_Bep+angle_Cep) local angldif = angle_BAC - (angle_BAep+angle_CAep)
if math.abs(angldif) < 0.01 then if math.abs(angldif) < 0.01 then
ep = vector.subtract(ep, B) ep = vector.subtract(ep, B)
nep = vector.normalize(ep) nep = vector.normalize(ep)
local angle_Aep = math.acos(vector.scalarproduct(BA, nep)) local angle_ABep = math.acos(vector.scalarproduct(nBA, nep))
local angle_Cep = math.acos(vector.scalarproduct(BC, nep)) local angle_CBep = math.acos(vector.scalarproduct(nBC, nep))
local angldif = angle_AC - (angle_Aep+angle_Cep) local angldif = angle_ABC - (angle_ABep+angle_CBep)
if math.abs(angldif) < 0.01 then if math.abs(angldif) < 0.01 then
table.insert(area, vector.add(pos, p)) table.insert(area, vector.add(pos, p))
end end