itkfemelement2dquadratictriangulartest.cxx

DTMK软件开发包,此为开源软件,是一款很好的医学图像开发资源.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkFEMElement2DQuadraticTriangularTest.cxx,v $
  Language:  C++
  Date: $Date: 2003-12-12 14:58:00 $
  Version:   $Revision: 1.4 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/

// disable debug warnings in MS compiler
#ifdef _MSC_VER
#pragma warning(disable: 4786)
#endif

#include "itkFEMElement2DC0QuadraticTriangularStrain.h"
#include "itkFEMElementBase.h"

#include <iostream>


//
int itkFEMElement2DQuadraticTriangularTest(int, char *[])
{
    typedef itk::fem::Node        NodeType;
    typedef itk::fem::Element     ElementType;

    NodeType::Pointer n0,n1,n2;
    ElementType::VectorType pt(2);

    n0=NodeType::New();
    pt[0]=0.;
    pt[1]=0.;
    n0->SetCoordinates(pt);

    n1=NodeType::New();
    pt[0]=1.;
    pt[1]=1.;
    n1->SetCoordinates(pt);

    n2=NodeType::New();
    pt[0]=0.;
    pt[1]=2.;
    n2->SetCoordinates(pt);

    typedef itk::fem::MaterialLinearElasticity   ElasticityType;
    
    ElasticityType::Pointer m = ElasticityType::New();
    
    m->GN=0;
    m->E=300.0;
    m->A=0.02;
    m->I=0.004;

    typedef itk::fem::Element2DC0QuadraticTriangularStrain StrainType;
    StrainType::Pointer e0 = StrainType::New();

    e0->GN=0;
    e0->SetNode(0, &*n0);
    e0->SetNode(1, &*n1);
    e0->SetNode(2, &*n2);
    e0->m_mat=dynamic_cast< ElasticityType * >(&*m);

    pt[0]=0.5;
    pt[1]=0.5;

    std::cout << "# integration points = " << e0->GetNumberOfIntegrationPoints(2) << std::endl;
    std::cout << "shape fxns at " << pt << ":\n" << e0->ShapeFunctions(pt) << std::endl;

    ElementType::MatrixType shapeD;
    e0->ShapeFunctionDerivatives(pt, shapeD);
    std::cout << "shape fxn derivatives:" << std::endl << shapeD << std::endl;

#ifndef FEM_USE_SMART_POINTERS
    delete e0;
    delete m;
    delete n0;
    delete n1;
    delete n2;
#endif

    
    std::cout << "Test PASSED!" << std::endl;
    return EXIT_SUCCESS;
}