In [ ]:

    

using PyPlot
include("../src/HPFEM.jl")
H = HPFEM


    

In [ ]:

    

b = H.Basis1d(5,7);


    

In [ ]:

    

Nel = 3
pts = [linspace(0, 5, Nel+1);]

elems = [H.Element1d(1, pts[i], pts[i+1], b) for i in 1:Nel];


    

In [ ]:

    

dof = H.DofMap1d(H.nmodes(b), Nel+1, []);
solver = H.CholeskySC(dof, H.BBSymTri);


    

In [ ]:

    

ii = H.interior_idx(b)
ib = H.bndry_idx(b);


    

In [ ]:

    

for e = 1:Nel
    println("======================================")
    println("PROCESSING ELEMENT ", )
    M = H.mass_matrix(b, elems[e])
    Abb = M[ib,ib]
    Aii = M[ii,ii]
    Abi = M[ib,ii]
    println(H.add_local_matrix(solver, e, Abb, Abi, Aii))
end


    

In [ ]:

    

solver.Abb


    

In [ ]: