In [1]:

    

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


    

    
Out[1]:





HPFEM

In [2]:

    

b = H.Basis1d(5,7);


    

In [3]:

    

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

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


    

In [4]:

    

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


    

In [5]:

    

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


    

In [6]:

    

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


    

    




======================================
PROCESSING ELEMENT 
-0.013888888888889006
======================================
PROCESSING ELEMENT 
-0.013888888888889006
======================================
PROCESSING ELEMENT 
-0.01388888888888895

In [7]:

    

solver.Abb


    

    
Out[7]:





HPFEM.BBSymTri{Float64}(4,4,[0.06944444444444431,0.13888888888888878,0.1388888888888889,0.06944444444444448],[-0.013888888888889006,-0.013888888888889006,-0.01388888888888895])

In [ ]: