1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkFEMLoadElementBase.cxx
Language: C++
Date: $Date$
Version: $Revision$
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 "itkFEMLoadElementBase.h"
namespace itk {
namespace fem {
/**
* Read the LoadElement object from input stream
*/
void LoadElement::Read( std::istream& f, void* info )
{
int n;
/**
* Convert the info pointer to a usable objects
*/
ReadInfoType::ElementArrayPointer elements=static_cast<ReadInfoType*>(info)->m_el;
/** first call the parent's read function */
Superclass::Read(f,info);
/**
* read and set pointers to element that we're applying the load to
*/
/** first we read number of pointers in a list */
this->SkipWhiteSpace(f); f>>n; if(!f) goto out;
if (n<=0)
{
/** if this is <= 0, the load applies on all elements in a system */
el.clear();
}
else
{
/**
* otherwise we read all the element numbers.
* there should be n of them
*/
for(int i=0;i<n;i++)
{
int m;
this->SkipWhiteSpace(f); f>>m; if(!f) goto out;
Element::ConstPointer e;
try
{
e=elements->Find(m);
}
catch ( FEMExceptionObjectNotFound exc )
{
throw FEMExceptionObjectNotFound(__FILE__,__LINE__,"LoadElementBase::Read()",exc.m_baseClassName,exc.m_GN);
}
el.push_back(e);
}
}
out:
if( !f )
{
throw FEMExceptionIO(__FILE__,__LINE__,"LoadElementBase::Read()","Error reading FEM load!");
}
}
/**
* Write the LoadElement to the output stream
*/
void LoadElement::Write( std::ostream& f ) const
{
/**
* first call the parent's write function
*/
Superclass::Write(f);
/** Write the list of element global numbers */
if (!el.empty())
{
f << "\t" <<static_cast<int>((el.size()));
f << "\t% # of elements on which the load acts" << std::endl;
f << "\t";
for(ElementPointersVectorType::const_iterator i=el.begin(); i != el.end(); i++)
{
f<<((*i)->GN)<<" ";
}
f << "\t% GNs of elements" << std::endl;
}
else
{
f << "\t-1\t% Load acts on all elements" << std::endl;
}
/** check for errors */
if (!f)
{
throw FEMExceptionIO(__FILE__,__LINE__,"LoadElement::Write()","Error writing FEM load!");
}
}
FEM_CLASS_REGISTER(LoadElement)
}} // end namespace itk::fem
|
