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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkFEMLoadGrav.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 "itkFEMLoadGrav.h"
namespace itk {
namespace fem {
/**
* Read the LoadGravConst object from input stream
*/
void LoadGravConst::Read( std::istream& f, void* info )
{
int n;
/** first call the parent's read function */
LoadGrav::Read(f,info);
/**
* Read and set the force vector
*/
/** first read and set the size of the vector */
this->SkipWhiteSpace(f); f>>n; if(!f) goto out;
Fg_value.set_size(n);
/** then the actual values */
this->SkipWhiteSpace(f); f>>Fg_value; if(!f) goto out;
out:
if( !f )
{
throw FEMExceptionIO(__FILE__,__LINE__,"LoadGravConst::Read()","Error reading FEM load!");
}
}
/**
* Write the LoadGravConst to the output stream
*/
void LoadGravConst::Write( std::ostream& f ) const
{
/** first call the parent's write function */
LoadGrav::Write(f);
/** then write the actual data force vector */
f<<"\t"<<Fg_value.size()<<"\t% Size of the gravity force vector\n";
f<<"\t"<<Fg_value<<"\t% Gravity force vector\n";
/** check for errors */
if (!f)
{
throw FEMExceptionIO(__FILE__,__LINE__,"LoadGravConst::Write()","Error writing FEM load!");
}
}
FEM_CLASS_REGISTER(LoadGravConst)
}} // end namespace itk::fem
|
