--Note: carbo star lemons and Xeno suggested the following
--alternate code for create_center.
function p.create_center(diam, bt)
    if diam >= 1 and diam <= 13 then
        local half = (diam - 1) >> 1
        local x1 = 7 - half
        local x2 = 7 + half
        if diam % 2 == 0 then
            x2 = x2 + 1
        end
        create_rect(bt, x1,x1,x1, x2,x2,x2)
    end
end
--This is what create_center does (with a given diameter):
-- 1  -> create_rect(bt, 7,7,7, 7,7,7)
-- 2  -> create_rect(bt, 7,7,7, 8,8,8)
-- 3  -> create_rect(bt, 6,6,6, 8,8,8)
-- 4  -> create_rect(bt, 6,6,6, 9,9,9)
-- 5  -> create_rect(bt, 5,5,5, 9,9,9)
-- 6  -> create_rect(bt, 5,5,5, 10,10,10)
-- 7  -> create_rect(bt, 4,4,4, 10,10,10)
-- 8  -> create_rect(bt, 4,4,4, 11,11,11)
-- 9  -> create_rect(bt, 3,3,3, 11,11,11)
-- 10 -> create_rect(bt, 3,3,3, 12,12,12)
-- 11 -> create_rect(bt, 2,2,2, 12,12,12)
-- 12 -> create_rect(bt, 2,2,2, 13,13,13)
-- 13 -> create_rect(bt, 1,1,1, 13,13,13)

--Creates a tube the the given radius
--along either the x,y or z axis.
--
--To do: Make the code look nice!
--Note a fancy way to do this is to have a virtual
--block coordinate (i,j,k) and turn it into an
--actual block coordinate (x,y,z) after specifying
--some direction.
--However we don't want to impose such a virtual system
--in the base package.
--Credit here goes to carbo star lemons.
function p.create_tube(diam, axis, bt)
    if diam >= 1 and diam <= 6 then
        local half = (diam - 1) >> 1
        local k1 = 7 - half
        local k2 = 7 + half
        
        if diam % 2 == 0 then
            k2 = k2 + 1
        end
        
        if(     axis == "x" ) then create_rect(bt, 0,k1,k1, 15,k2,k2)
        elseif( axis == "y" ) then create_rect(bt, k1,0,k1, k2,15,k2) 
        elseif( axis == "z" ) then create_rect(bt, k1,k1,0, k2,k2,15) end
    end
end

--Creates a "half tube" from one wall
--of the chunk to the center of the chunk.
--See the comment about create_tube on how to rewrite this.
--Credit here goes to carbo star lemons.
function p.create_half_tube(diam, dir, bt)
    if diam >= 1 and diam <= 6 then
        local half = (diam - 1) >> 1
        local i1 = 7 - half
        local i2 = 7 + half
        
        if diam % 2 == 0 then
            i2 = i2 + 1
        end
        
        if(     dir == "x_pos" ) then create_rect(bt, i1,i1,i1, 15,i2,i2) 
        elseif( dir == "x_neg" ) then create_rect(bt, 0,i1,i1, i2,i2,i2) 
        elseif( dir == "y_pos" ) then create_rect(bt, i1,i1,i1, i2,15,i2) 
        elseif( dir == "y_neg" ) then create_rect(bt, i1,0,i1, i2,i2,i2) 
        elseif( dir == "z_pos" ) then create_rect(bt, i1,i1,i1, i2,i2,15) 
        elseif( dir == "z_neg" ) then create_rect(bt, i1,i1,0, i2,i2,i2) end         
    end
end

--Will create the 12 edges of the chunk.
function p.create_edges(bt)
    create_rect(bt, 0,0,0,   15,0,0)
    create_rect(bt, 0,0,15,  15,0,15)
    create_rect(bt, 0,15,0,  15,15,0)
    create_rect(bt, 0,15,15, 15,15,15)
    --
    create_rect(bt, 0,0,0,   0,15,0)
    create_rect(bt, 0,0,15,  0,15,15)
    create_rect(bt, 15,0,0,  15,15,0)
    create_rect(bt, 15,0,15, 15,15,15)
    --
    create_rect(bt, 0,0,0,   0,0,15)
    create_rect(bt, 0,15,0,  0,15,15)
    create_rect(bt, 15,0,0,  15,0,15)
    create_rect(bt, 15,15,0, 15,15,15)
end

--Will create the outer shell of the chunk.
function p.create_shell(bt)
    create_rect(bt, 0,0,0,  0,15,15)
    create_rect(bt, 15,0,0, 15,15,15)
    create_rect(bt, 0,0,0,  15,0,15)
    create_rect(bt, 0,15,0, 15,15,15)
    create_rect(bt, 0,0,0,  15,15,0)
    create_rect(bt, 0,0,15, 15,15,15)
end

--Can be use in conjunction with create_shell
--to make a room with doors in it.
--Thanks to Xeno who found a problem with this function.
function p.create_2x2_door(dir, rim_type, hole_type)
    if(dir == "x_pos") then
        create_rect(rim_type, 15,6,6, 15,9,9)
        create_rect(hole_type, 15,7,7, 15,8,8)
    elseif(dir == "x_neg") then
        create_rect(rim_type, 0,6,6, 0,9,9)
        create_rect(hole_type, 0,7,7, 0,8,8)
    elseif(dir == "y_pos") then
        create_rect(rim_type, 6,15,6, 9,15,9)
        create_rect(hole_type, 7,15,7, 8,15,8)
    elseif(dir == "y_neg") then
        create_rect(rim_type, 6,0,6, 9,0,9)
        create_rect(hole_type, 7,0,7, 8,0,8)
    elseif(dir == "z_pos") then
        create_rect(rim_type, 6,6,15, 9,9,15)
        create_rect(hole_type, 7,7,15, 8,8,15)
    elseif(dir == "z_neg" ) then
        create_rect(rim_type, 6,6,0, 9,9,0)
        create_rect(hole_type, 7,7,0, 8,8,0)
    end
end

--VECTORS:

function p.vec(x_value, y_value, z_value)
    return {x=x_value, y=y_value, z=z_value}
end

function p.vec_zero()
    return p.vec(0.0, 0.0, 0.0)
end

function p.dot(a,b)
    return a.x*b.x + a.y*b.y + a.z*b.z
end

function p.length(v)
    return math.sqrt(v.x*v.x + v.y*v.y + v.z*v.z)
end

function p.normalize(v)
    local len = p.length(v)
    if( len < 0.00001 ) then return p.vec_zero() end
    local inv_len = 1.0 / len
    p.vec_scale_eq(v,inv_len)
end

function p.vec_add(v1, v2)
    return p.vec(v1.x+v2.x, v1.y+v2.y, v1.z+v2.z)
end

function p.vec_add_eq(v1, v2)
    v1.x = v1.x + v2.x
    v1.y = v1.y + v2.y
    v1.z = v1.z + v2.z
    return v1
end

function p.vec_sub(v1, v2)
    return p.vec(v1.x-v2.x, v1.y-v2.y, v1.z-v2.z)
end

function p.dist(v1, v2)
    local diff = p.vec_sub(v1, v2)
    return p.length(diff)
end

function p.vec_sub_eq(v1, v2)
    v1.x = v1.x - v2.x
    v1.y = v1.y - v2.y
    v1.z = v1.z - v2.z
    return v1
end

function p.vec_scale(v, c)
    return p.vec(v.x*c, v.y*c, v.z*c)
end

function p.vec_scale_eq(v, c)
    v.x = v.x * c
    v.y = v.y * c
    v.z = v.z * c
    return v
end

function p.angle_between_vec(a, b)
    local arg =
        p.dot(a,b) / ( p.length(a) * p.length(b) )
    local screw_acos = 0.9999
    if( arg > screw_acos ) then arg = screw_acos end
    return math.acos(arg)
    --Note: acos(1) is essentially undefined.
end