In [201]:

    

using Jacobi
using PyPlot
include("../src/HPFEM.jl");


    

    




WARNING: replacing module HPFEM

In [202]:

    

M=50
Q=M
lagr = HPFEM.Lagrange1d(M);
quad = HPFEM.QuadType(Q);
base = HPFEM.Basis1d(lagr, quad);


    

In [203]:

    

elemento = HPFEM.Element1d(1, -1.0, 1.0, base);
Mass = HPFEM.mass_matrix(base, elemento);


    

In [204]:

    

fun(x) = sin(2*pi*x)
x = linspace(-1,1,101)
#plot(x,fun(x));


    

In [205]:

    

Fe = zeros(Q)
fe = fun(elemento.ξ)
HPFEM.add_rhs!(base, elemento, fe, Fe)
uh = Mass\Fe
#plot(uh)
uh - fe


    

    
Out[205]:





50-element Array{Float64,1}:
 -4.93038e-32
 -3.46945e-18
  1.38778e-17
  0.0        
  0.0        
  0.0        
  0.0        
  0.0        
  0.0        
  1.11022e-16
  0.0        
  0.0        
  0.0        
  ⋮          
  0.0        
  0.0        
 -1.11022e-16
  0.0        
  0.0        
  0.0        
  0.0        
  0.0        
  0.0        
 -1.38778e-17
  3.46945e-18
  4.93038e-32

In [206]:

    

Q2 = 101
x =linspace(-1,1,Q2)


    

    
Out[206]:





linspace(-1.0,1.0,101)

In [207]:

    

#u = zeros(Q)
test = zeros(101,M)
for i in 1 :Q
    test[:,i] = base(x,i)
end
plot(test*uh,"o")
plot(fun(x))


    

    













    
Out[207]:





1-element Array{Any,1}:
 PyObject <matplotlib.lines.Line2D object at 0x7f29d8be1f10>