DeletingP1D.m

classdef DeletingP1D < handle
    
    properties (Access = public)
        
    end
    
    properties (Access = private)
        
    end
    
    properties (Access = private)
       mesh
       fun
       funP1
    end
    
    methods (Access = public)
        
        function obj = DeletingP1D()
            obj.init()
            obj.createMesh();
            obj.createAnalyticalFunciton();
            obj.createP1ContinousFun();
            obj.createP1DiscontinousFun();
        end
        
    end
    
    methods (Access = private)
        
        function init(obj)
            
        end
        
        function createMesh(obj)
            m = TriangleMesh(1,1,20,20);
            obj.mesh = m;
        end
        
        function createAnalyticalFunciton(obj)
            s.fHandle = @(x) [cos(10*x(1,:,:)+x(2,:,:));sin(10*x(1,:,:)+x(2,:,:))];
            s.ndimf  = 2;
            s.mesh   = obj.mesh;
            obj.fun = AnalyticalFunction(s);
        end

        function createP1ContinousFun(obj)
            obj.funP1 = obj.fun.project('P1');
        end

        function createP1DiscontinousFun(obj)
            funP1DC = obj.funP1.project('P1D');

            obj.funP1.computeL2norm()
            funP1DC.computeL2norm()

            obj.funP1.plot()
            funP1DC.plot()

            gradP1DC = SymGrad(funP1DC).project('P1');
            gradP1DC.plot()

            gradP1 = SymGrad(obj.funP1).project('P1');
            gradP1.plot()

            dif = gradP1DC - gradP1;
            dif.computeL2norm()/gradP1.computeL2norm()
        end        
        
        
        
        
        
    end
    
end