runValids.m

Go to the documentation of this file.

function runValids(varargin)
% function runValids(varargin)
%   Run validation tests for Mass (#validMassP1), MassW (#validMassWP1),
%   Stiff (#validStiffP1) and StiffElas (#validStiffElasP1) matrices.
%   For each assembly matrices, we make 3 tests 
%     - Test 1 : computation of the assembly Matrix using all the versions giving 
%       errors and cputimes.
%     - Test 2 : comparison of the associated integral and its `P_1`-Lagrange approximation.
%     - Test 3 : retrieving of numerical order for `P_1`-Lagrange approximation.
%
% Parameters:
%   save    :  set to true to save figures in 'image' directory.
%   percent :  value for resizing the figure (only png format). See #SaveFigure
%
% Example:
%  @verbatim
%     runValids('save',true,'percent',50)@endverbatim
% 
% Results under Matlab R2012b:
%   Here are the 'Test 3' figures obtained on our reference machine . 
%  - validMassP1() figure
%  \image html images/validMassP1.png "figure : Mass Matrix order validation" 
%  - validMassWP1() figure
%  \image html images/validMassWP1.png "figure : MassW Matrix order validation"
%  - validStiffP1() figure
%  \image html images/validStiffP1.png "figure : Stiff Matrix order validation"
%  - validStiffElasP1() figure
%  \image html images/validStiffElasP1.png "figure : StiffElas Matrix order validation"
%
% See also:
%   #validMassP1, #validMassWP1, #validStiffP1, #validStiffElasP1, #SaveFigure
%
% OptFEM2DP1 [V1.2b1] -    Copyright (C) 2013  CJS (LAGA)
%
%   This file is part of OptFEM2DP1.
%   OptFEM2DP1 is free software: you can redistribute it and/or modify
%   it under the terms of the GNU General Public License as published by
%   the Free Software Foundation, either version 3 of the License, or
%   (at your option) any later version.
%
%   OptFEM2DP1 is distributed in the hope that it will be useful,
%   but WITHOUT ANY WARRANTY; without even the implied warranty of
%   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
%   GNU General Public License for more details.
%
%   You should have received a copy of the GNU General Public License
%   along with this program.  If not, see <http://www.gnu.org/licenses/>.
InitOptFEM2D();

p = inputParser; 
  
if isOctave()
  p=p.addParamValue('save', false, @islogical );
  p=p.addParamValue('percent', 50 , @(t) ((t>0)&&(t<=100)) );
  p=p.parse(varargin{:});
else % Matlab
  p.addParamValue('save', false, @islogical );
  p.addParamValue('percent', 50 , @(t) ((t>0)&&(t<=100)) );
  p.parse(varargin{:});
end
close all

figure(1)
validMassP1();
SaveFigure(p.Results.save,'validMassP1',p.Results.percent)
figure(2)
validMassWP1();
SaveFigure(p.Results.save,'validMassWP1',p.Results.percent)
figure(3)
validStiffP1();
SaveFigure(p.Results.save,'validStiffP1',p.Results.percent)
figure(4)
validStiffElasP1(0);
SaveFigure(p.Results.save,'validStiffElasP1',p.Results.percent)