ITK  4.4.0
Insight Segmentation and Registration Toolkit
itkFEMLoadGrav.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright Insight Software Consortium
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * http://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 
19 #ifndef __itkFEMLoadGrav_h
20 #define __itkFEMLoadGrav_h
21 
22 #include "itkFEMLoadElementBase.h"
23 #include "vnl/vnl_vector.h"
24 
25 namespace itk
26 {
27 namespace fem
28 {
39 class LoadGrav : public LoadElement
40 {
41 public:
42 
44  typedef LoadGrav Self;
48 
50  itkTypeMacro(LoadGrav, LoadElement);
51 
52  virtual vnl_vector<Float> GetGravitationalForceAtPoint(vnl_vector<Float> ) = 0;
53 
54 protected:
55  virtual void PrintSelf(std::ostream& os, Indent indent) const
56  {
57  Superclass::PrintSelf(os, indent);
58  }
59 
60 };
61 
70 class LoadGravConst : public LoadGrav
71 {
72 public:
73 
79 
81  itkSimpleNewMacro(Self);
82 
84  itkTypeMacro(LoadGravConst, LoadGrav);
85 
88  virtual::itk::LightObject::Pointer CreateAnother(void) const;
89 
90  virtual vnl_vector<Float> GetGravitationalForceAtPoint(vnl_vector<Float> )
91  {
92  return m_GravityForce;
93  }
94 
98  void SetForce(const vnl_vector<itk::fem::Element::Float> force);
99 
103  vnl_vector<itk::fem::Element::Float> & GetForce();
104  const vnl_vector<itk::fem::Element::Float> & GetForce() const;
106 
108  virtual void ApplyLoad(Element::ConstPointer element, Element::VectorType & Fe);
109 
110 protected:
111  virtual void PrintSelf(std::ostream& os, Indent indent) const;
112 
113  vnl_vector<Float> m_GravityForce;
114 };
115 
116 }
117 } // end namespace itk::fem
118 
119 #endif // #ifndef __itkFEMLoadGrav_h
120