SolveEikonal

% Class for solving the Eikonal equation on a 2-D triangulation
classdef SolveEikonal
    properties %(SetAccess='private',GetAccess='private')
        TM        % triangle mesh object
        Bdy       % data struct containing boundary vertices and dirichlet solution data
        Param     % optional parameters
    end
    methods
        function obj = SolveEikonal(varargin)
            
            if nargin ~= 3
                disp('Requires 3 arguments!');
                disp('First is a struct containing a triangulation.');
                disp('Second is a struct containing the boundary data.');
                disp('Third contains some parameters for controlling the algorithm.');
                error('Check the input data!');
            end
            
            LINE_str = '-------------------------------';
            OUT_str  = '| Init SolveEikonal Object... |';
            disp(LINE_str);
            disp(OUT_str);
            disp(LINE_str);
            
            obj.TM        = varargin{1};
            obj.Bdy       = varargin{2};
            obj.Param     = varargin{3};
            
            % check input
            fields_TF = isfield(obj.TM,{'Vtx','DoFmap','NegMask'});
            if (min(fields_TF)~=1)
                error('Invalid mesh struct...');
            end
            if ~isempty(obj.TM.NegMask)
                if (length(obj.TM.NegMask)~=size(obj.TM.Vtx,1))
                    error('Length of negative mask must be the same as the number of mesh vertices.');
                end
            end
            
            %%%%%
            
            fields_TF = isfield(obj.Bdy,{'Nodes','Data'});
            if (min(fields_TF)~=1)
                error('Invalid boundary data struct...');
            end
            if (length(obj.Bdy.Nodes)~=length(obj.Bdy.Data))
                error('Given Dirichlet data is not consistent: lengths are different.');
            end
            if (max(obj.Bdy.Nodes)>size(obj.TM.Vtx,1))
                error('Given Dirichlet data is not consistent: more Dirichlet nodes than mesh vertices.');
            end
            
            %%%%%
            
            fields_TF = isfield(obj.Param,{'NumGaussSeidel','INF_VAL','TOL'});
            if (min(fields_TF)~=1)
                error('Invalid parameter list...');
            end
            
            %%%%%

            % init neighboring vertex/triangle data
            obj.TM.TriStarData = Get_Triangle_Neighbors(size(obj.TM.Vtx,1),obj.TM.DoFmap);
            
        end
    end
end

% END %