main_runUnitTests.m

%% Licensing
%
% License:         BSD License
%                  cane Multiphysics default license: cane/license.txt
%
% Main authors:    Andreas Apostolatos
%
%% Script documentation
% 
% Task : Run unitTest cases for the following type of analysis,
%        - Numerical quadrature
%        - Isogeometric beam analysis
%        - Isogeometric membrane analysis
%        - Isogeometric Kichhoff-Love shell analysis
%        - Finite element formulation of thermal conduction analysis
%        - Finite element formulation of plate in membrane action analysis
%        - Finite element formulation of the frictionless Signorini problem
%        - Stabilized isogeometric analysis for incompressible flows
%        - Stabilized finite element analysis for the Navier-Stokes problem
%        - Finite element methods for Fluid-Structure Interaction analysis
%        - Shape optimization analysis
%
% Date : 16.05.2020
%
%% Clear memnory and command window
clc;
clear;

%% Includes

% Add transient analysis functions
addpath('../../transientAnalysis/');

% Add general math functions
addpath('../../generalMath/');

% Add general auxiliary functions
addpath('../../auxiliary/');

% Add all fast computation functions
addpath('../../efficientComputation/');

% Add the equation system solvers
addpath('../../equationSystemSolvers/');

% Add the classical finite element basis functions
addpath('../../basisFunctions/');

% Add all functions related to the Computer-Aided Geometric Design (GACD) 
% kernel
addpath('../../CAGDKernel/CAGDKernel_basisFunctions/',...
        '../../CAGDKernel/CAGDKernel_geometryResolutionRefinement/',...
        '../../CAGDKernel/CAGDKernel_baseVectors/',...
        '../../CAGDKernel/CAGDKernel_graphics/',...
        '../../CAGDKernel/CAGDKernel_BSplineCurve/',...
        '../../CAGDKernel/CAGDKernel_BSplineSurface/');

% Add all functions related to the isogeometric beam analysis
addpath('../../isogeometricBeamAnalysis/stiffnessMatrices/',...
        '../../isogeometricBeamAnalysis/graphics/',...
        '../../isogeometricBeamAnalysis/loads/',...
        '../../isogeometricBeamAnalysis/postprocessing/',...
        '../../isogeometricBeamAnalysis/solvers/',...
        '../../isogeometricBeamAnalysis/math/',...
        '../../isogeometricBeamAnalysis/auxiliary/',...
        '../../isogeometricBeamAnalysis/errorComputation/');
    
% Add all functions related to the isogeometric Kirchhoff-Love shell 
% formulation
addpath('../../isogeometricThinStructureAnalysis/graphicsSinglePatch/',...
        '../../isogeometricThinStructureAnalysis/graphicsMultipatches/',...
        '../../isogeometricThinStructureAnalysis/loads/',...
        '../../isogeometricThinStructureAnalysis/solutionMatricesAndVectors/',...
        '../../isogeometricThinStructureAnalysis/solvers/',...
        '../../isogeometricThinStructureAnalysis/metrics/',...
        '../../isogeometricThinStructureAnalysis/auxiliary/',...
        '../../isogeometricThinStructureAnalysis/postprocessing/',...
        '../../isogeometricThinStructureAnalysis/BOperatorMatrices/',...
        '../../isogeometricThinStructureAnalysis/penaltyDecompositionKLShell/',...
        '../../isogeometricThinStructureAnalysis/penaltyDecompositionMembrane/',...
        '../../isogeometricThinStructureAnalysis/lagrangeMultipliersDecompositionKLShell/',...
        '../../isogeometricThinStructureAnalysis/nitscheDecompositionMembrane/',...
        '../../isogeometricThinStructureAnalysis/errorComputation/',...
        '../../isogeometricThinStructureAnalysis/transientAnalysis/',...
        '../../isogeometricThinStructureAnalysis/initialConditions/',...
        '../../isogeometricThinStructureAnalysis/weakDBCMembrane/',...
        '../../isogeometricThinStructureAnalysis/formFindingAnalysis/');

% Add all functions related to heat transfer analysis
addpath('../../FEMThermalConductionAnalysis/solvers/',...
        '../../FEMThermalConductionAnalysis/solutionMatricesAndVectors/',...
        '../../FEMThermalConductionAnalysis/loads/',...
        '../../FEMThermalConductionAnalysis/graphics/',...
        '../../FEMThermalConductionAnalysis/output/',...
        '../../FEMThermalConductionAnalysis/initialConditions/',...
        '../../FEMThermalConductionAnalysis/postprocessing/', ...
        '../../FEMThermalConductionAnalysis/transientAnalysis/');
    
% Add all functions related to plate in membrane action analysis
addpath('../../FEMPlateInMembraneActionAnalysis/solvers/',...
        '../../FEMPlateInMembraneActionAnalysis/solutionMatricesAndVectors/',...
        '../../FEMPlateInMembraneActionAnalysis/loads/',...
        '../../FEMPlateInMembraneActionAnalysis/graphics/',...
        '../../FEMPlateInMembraneActionAnalysis/output/',...
        '../../FEMPlateInMembraneActionAnalysis/postprocessing/',...
        '../../FEMPlateInMembraneActionAnalysis/initialConditions/',...
        '../../FEMPlateInMembraneActionAnalysis/errorComputation/');

% Add all functions related to Signorini frictionless contact problem
addpath('../../FEMContactMechanicsAnalysis/graphics/',...
        '../../FEMContactMechanicsAnalysis/solvers',...
        '../../FEMContactMechanicsAnalysis/auxiliary/',...
        '../../FEMContactMechanicsAnalysis/contactSegments/',...
        '../../FEMContactMechanicsAnalysis/postprocessing/');
    
% Add all functions related to the isogeometric Computational Fluid
% Dynamics problems
addpath('../../isogeometricComputationalFluidDynamicsAnalysis/solutionMatricesAndVectors/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/solvers/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/neumannBoundaryConditions/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/graphics/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/postProcessing/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/inhomogeneousDirichletBoundaryConditions/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/initialConditions/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/transientAnalysis/',...
        '../../isogeometricComputationalFluidDynamicsAnalysis/errorComputation/');
    
% Add all functions related to the Finite Element Methods for Computational
% Fluid Dynamics problems
addpath('../../FEMComputationalFluidDynamicsAnalysis/solutionMatricesAndVectors/',...
        '../../FEMComputationalFluidDynamicsAnalysis/boundaryConditions/',...
        '../../FEMComputationalFluidDynamicsAnalysis/initialConditions',...
        '../../FEMComputationalFluidDynamicsAnalysis/solvers/',...
        '../../FEMComputationalFluidDynamicsAnalysis/loads/',...
        '../../FEMComputationalFluidDynamicsAnalysis/ALEMotion/',...
        '../../FEMComputationalFluidDynamicsAnalysis/transientAnalysis/',...
        '../../FEMComputationalFluidDynamicsAnalysis/postProcessing/');
    
% Add all functions related to the Finite Element Methods for
% Fluid-Structure Interaction
addpath('../../FEMComputationalFluidStructureInteractionAnalysis/ALEMotion/', ...
        '../../FEMComputationalFluidStructureInteractionAnalysis/solvers/');
    
% Add all functions related to parsing
addpath('../../parsers/');

% Add all unit test functions and classes
addpath('../../unitTest/');
    
%% Global variables for the unit tests
isLight = true;

%% Import the matlab test suite
import matlab.unittest.TestSuite;
import matlab.unittest.selectors.HasName;
import matlab.unittest.constraints.EndsWithSubstring;

%% Run the unit test cases for the quadratures
suiteQuadrature = TestSuite.fromClass(?testQuadrature);
resultQuadrature = run(suiteQuadrature);

%% Run the unit test cases for the utilities functions
suiteClassUtilityFunctions = TestSuite.fromClass(?testUtilityFunctions);
resultIGAUtilityFunctions = run(suiteClassUtilityFunctions);

%% Run the unit test cases for the isogeometric beam analysis
suiteClassIGABeam = TestSuite.fromClass(?testIGABeamAnalysis);
resultIGABeam = run(suiteClassIGABeam);

%% Run the unit test cases for the isogeometric membrane analysis
suiteClassIGAMembrane = TestSuite.fromClass(?testIGAMembraneAnalysis);
if isLight
    suiteClassIGAMembrane = suiteClassIGAMembrane.selectIf...
        (HasName('testIGAMembraneAnalysis/testPagewiseComputationTangStiffMtxIGAMembrane') | ...
        HasName('testIGAMembraneAnalysis/testNitscheWeakDBCPlateMembraneAnalysis') | ...
        HasName('testIGAMembraneAnalysis/testTransientNonlinearMembraneAnalysis') | ...
        HasName('testIGAMembraneAnalysis/testTransientNitscheWeakDBCPlateMembraneAnalysis') | ...
        HasName('testIGAMembraneAnalysis/testDDMNitschePlateMembraneAnalysis') | ...
        HasName('testIGAMembraneAnalysis/testFoFiMembraneAnalysis') | ...
        HasName('testIGAMembraneAnalysis/testDDMFourPointSailAnalysis'));
end
resultIGAMembrane = run(suiteClassIGAMembrane);

%% Run the unit test cases for the isogeometric Kichhoff-Love shell analysis
warning('Test testLinearKirchoffLoveShellMultipatchAnalysis needs to be fixed');
suiteClassIGAKLShell = TestSuite.fromClass(?testIGAKirchhoffLoveShellAnalysis);
if isLight
    suiteClassIGAKLShell = suiteClassIGAKLShell.selectIf...
        (HasName('testIGAKirchhoffLoveShellAnalysis/testLinearKirchoffLoveShellAnalysis') | ...
        HasName('testIGAKirchhoffLoveShellAnalysis/testNonlinearKirchoffLoveShellAnalysis') | ...
        HasName('testIGAKirchhoffLoveShellAnalysis/testLinearKirchoffLoveShellMultipatchAnalysis'));
end
resultIGAKLShell = run(suiteClassIGAKLShell);

%% Run the unit test cases for the finite element formulation of the thermal conduction analysis
suiteClassFEMThermalConduction = TestSuite.fromClass(?testFEMThermalConductionAnalysis);
resultFEMThermalConduction = run(suiteClassFEMThermalConduction);

%% Run the unit test cases for the finite element formulation of the plate in membrane action analysis
suiteClassFEMPlateInMembraneAction = TestSuite.fromClass(?testFEMPlateInMembraneActionAnalysis);
resultFEMPlateInMembraneAction = run(suiteClassFEMPlateInMembraneAction);

%% Run the unit test cases for the finite element formulation of the contact mechanics analysis
suiteClassFEMContactMechanics = TestSuite.fromClass(?testFEMContactMechanicsAnalysis);
resultContactMechanicsAnalysis = run(suiteClassFEMContactMechanics);

%% Run the unit test cases for the stabilized isogeometric incompressible flow equations
suiteClassIGA4CFD = TestSuite.fromClass(?testIGAComputationalFluidDynamicsAnalysis);
resultIGA4CFD = run(suiteClassIGA4CFD);

%% Run the unit test cases for stabilized finite element formulation for the Navier-Stokes problem
suiteClassFEM4CFD = TestSuite.fromClass(?testFEMComputationalFluidDynamicsAnalysis);
resultFEM4CFD = run(suiteClassFEM4CFD);

%% Run the unit test cases for finite element methods for Fluid-Structure Interaction analysis
suiteClassFEM4FSI = TestSuite.fromClass(?testFEMComputationalFluidStructureInteractionAnalysis);
resultFEM4FSI = run(suiteClassFEM4FSI);

%% Run the unit test cases for shape optimization analysis
suiteClassShapeOptimization = TestSuite.fromClass(?testShapeOptimizationAnalysis);
resultShapeOptimization = run(suiteClassShapeOptimization);

%% END OF SCRIPT