runValids.m

Go to the documentation of this file.

function runValids(varargin)
% function runValids(varargin)
%   Run validation tests for Mass (#validMass3DP1), Stiff (#validStiffP1),
%   StiffElas (#validStiffElas3DP1) and MassVF (#validMassVF3DP1) matrices.
%   For each assembly matrix, we make 3 tests 
%     - Test 1 : computation of the assembly Matrix using all the versions giving 
%       errors and computation times.
%     - 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 figures for 'Test 3' obtained on our reference machine. 
%  - validMass3DP1() figure
%  \image html images/validMass3DP1.png "figure : Mass Matrix order validation" 
%  - validStiff3DP1() figure
%  \image html images/validStiff3DP1.png "figure : Stiff Matrix order validation"
%  - validStiffElas3DP1(0) figure
%  \image html images/validStiffElas3DP1.png "figure : StiffElas Matrix order validation"
%  - validMassVF3DP1(0) figure
%  \image html images/validMassVF3DP1.png "figure : MassVF Matrix order validation"
%
% See also:
%   #validMass3DP1, #validStiff3DP1, #validStiffElas3DP1, #validMassVF3DP1, #SaveFigure
% Copyright:
%   See \ref license
InitOptFEM3D();

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)
validMass3DP1();
SaveFigure(p.Results.save,'validMass3DP1',p.Results.percent)
figure(2)
validMassVF3DP1(0);
SaveFigure(p.Results.save,'validMassVF3DP1',p.Results.percent)
figure(3)
validStiff3DP1();
SaveFigure(p.Results.save,'validStiff3DP1',p.Results.percent)
figure(4)
validStiffElas3DP1(0);
SaveFigure(p.Results.save,'validStiffElas3DP1',p.Results.percent)