Developer Documentation
OVMPropertyVisualizerVectorFieldDifferenceT_impl.hh
1 /*===========================================================================*\
2 * *
3 * OpenFlipper *
4  * Copyright (c) 2001-2015, RWTH-Aachen University *
5  * Department of Computer Graphics and Multimedia *
6  * All rights reserved. *
7  * www.openflipper.org *
8  * *
9  *---------------------------------------------------------------------------*
10  * This file is part of OpenFlipper. *
11  *---------------------------------------------------------------------------*
12  * *
13  * Redistribution and use in source and binary forms, with or without *
14  * modification, are permitted provided that the following conditions *
15  * are met: *
16  * *
17  * 1. Redistributions of source code must retain the above copyright notice, *
18  * this list of conditions and the following disclaimer. *
19  * *
20  * 2. Redistributions in binary form must reproduce the above copyright *
21  * notice, this list of conditions and the following disclaimer in the *
22  * documentation and/or other materials provided with the distribution. *
23  * *
24  * 3. Neither the name of the copyright holder nor the names of its *
25  * contributors may be used to endorse or promote products derived from *
26  * this software without specific prior written permission. *
27  * *
28  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS *
29  * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED *
30  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A *
31  * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER *
32  * OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, *
33  * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, *
34  * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR *
35  * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF *
36  * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING *
37  * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS *
38  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. *
39 * *
40 \*===========================================================================*/
41 
42 
43 
44 
45 #define OVM_PROPERTY_VISZUALIZTER_VECTOR_FIELD_DIFFERENCE_CC
46 
47 #include "OVMPropertyVisualizerVectorFieldDifference.hh"
48 
49 template <typename MeshT>
50 OVMPropertyVisualizerVectorFieldDifference<MeshT>::OVMPropertyVisualizerVectorFieldDifference(MeshT* _mesh, int objectID, PropertyInfo _propertyInfo1, PropertyInfo _propertyInfo2)
51  : OVMPropertyVisualizer<MeshT>(_mesh, objectID, _propertyInfo1),
52  propertyInfo2(_propertyInfo2)
53 {
54  if (PropertyVisualizer::widget) delete PropertyVisualizer::widget;
55  VectorFieldDifferenceWidget* w = new VectorFieldDifferenceWidget();
56  w->paramVectorFieldDifference->setTitle(QString("3D Vector Field Difference Parameters of ").append(PropertyVisualizer::propertyInfo.propName().c_str()));
57  PropertyVisualizer::widget = w;
58 }
59 
60 namespace OVMPVVFD{ //OVMPVVFD = OpenMeshProeprtyVisualizerVectorFieldDifference
61 
62 template <typename MeshT>
63 float scalarFn_norm_of_diff(const ACG::Vec3d &a, const ACG::Vec3d &b) {
64  return (a - b).norm();
65 }
66 
67 template <typename MeshT>
68 float scalarFn_diff_of_norms(const ACG::Vec3d &a, const ACG::Vec3d &b) {
69  return std::fabs(a.norm() - b.norm());
70 }
71 
72 template <typename MeshT>
73 float scalarFn_4_symm_diff(const ACG::Vec3d &a, const ACG::Vec3d &b) {
74  double alpha = std::acos((a|b) / a.norm() / b.norm());
75  alpha -= std::floor((alpha + M_PI_4) / M_PI_2) * M_PI_2;
76  return std::fabs(alpha);
77 }
78 
79 
80 template<typename Prop1, typename Prop2, float (*ScalarFn)(const typename Prop1::value_type &, const typename Prop2::value_type &)>
81 class ScalarAssigner {
82  public:
83  ScalarAssigner(const Prop1 &prop1, const Prop2 &prop2) :
84  prop1(prop1), prop2(prop2) { }
85 
86  template<typename Handle>
87  float operator() (const Handle &handle) const {
88  return ScalarFn(prop1[handle], prop2[handle]);
89  }
90 
91  protected:
92  const Prop1 &prop1;
93  const Prop2 &prop2;
94 };
95 
96 template<typename MeshT, typename IteratorT, typename PropHandle, float (*ScalarFn)(const ACG::Vec3d &, const ACG::Vec3d &)>
97 void colorElements(int objectID,
98  PropHandle prop1,
99  PropHandle prop2,
100  IteratorT primitivesBegin,
101  IteratorT primitivesEnd)
102 {
103 
104  std::vector<float> scalars;
105  std::transform(primitivesBegin, primitivesEnd, std::back_inserter(scalars),
106  ScalarAssigner<PropHandle, PropHandle, ScalarFn >(prop1, prop2));
107 
108  const float min = *std::min_element(scalars.begin(), scalars.end());
109  const float max = *std::max_element(scalars.begin(), scalars.end());
110 
111  VolumeMeshObject<MeshT>* object;
112  PluginFunctions::getObject(objectID, object);
113 
114  ACG::ColorCoder colCod(min,max,false);
115  for (std::vector<float>::iterator i = scalars.begin(); i != scalars.end(); ++i)
116  object->colors()[*(primitivesBegin++)] = colCod(*i);
117 }
118 
119 }
120 
121 
122 template <typename MeshT>
123 void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeCellProp(bool _setDrawMode)
124 {
125  MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
126 
127  OpenVolumeMesh::CellPropertyT<ACG::Vec3d> prop1 = mesh->template request_cell_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
128  OpenVolumeMesh::CellPropertyT<ACG::Vec3d> prop2 = mesh->template request_cell_property<ACG::Vec3d>(propertyInfo2.propName());
129 
130  VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
131 
132  if (w->vecFieldDiff_4symm_rb->isChecked())
133  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::CellIter, OpenVolumeMesh::CellPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->cells_begin(), mesh->cells_end());
134  if (w->vecFieldDiff_diff_norm_rb->isChecked())
135  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::CellIter, OpenVolumeMesh::CellPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->cells_begin(), mesh->cells_end());
136  if (w->vecFieldDiff_norm_diff_rb->isChecked())
137  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::CellIter, OpenVolumeMesh::CellPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->cells_begin(), mesh->cells_end());
138  if (_setDrawMode)
139  {
140  VolumeMeshObject<MeshT>* object;
141  PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
142  object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.cellsColoredPerCell);
143  }
144 }
145 
146 template <typename MeshT>
147 void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeFaceProp(bool _setDrawMode)
148 {
149  MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
150 
151  OpenVolumeMesh::FacePropertyT<ACG::Vec3d> prop1 = mesh->template request_face_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
152  OpenVolumeMesh::FacePropertyT<ACG::Vec3d> prop2 = mesh->template request_face_property<ACG::Vec3d>(propertyInfo2.propName());
153 
154  VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
155 
156  if (w->vecFieldDiff_4symm_rb->isChecked())
157  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::FaceIter, OpenVolumeMesh::FacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->faces_begin(), mesh->faces_end());
158  if (w->vecFieldDiff_diff_norm_rb->isChecked())
159  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::FaceIter, OpenVolumeMesh::FacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->faces_begin(), mesh->faces_end());
160  if (w->vecFieldDiff_norm_diff_rb->isChecked())
161  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::FaceIter, OpenVolumeMesh::FacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->faces_begin(), mesh->faces_end());
162  if (_setDrawMode)
163  {
164  VolumeMeshObject<MeshT>* object;
165  PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
166  object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.facesColoredPerFace);
167  }
168 }
169 
170 template <typename MeshT>
171 void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeHalffaceProp(bool _setDrawMode)
172 {
173  MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
174 
175  OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d> prop1 = mesh->template request_halfface_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
176  OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d> prop2 = mesh->template request_halfface_property<ACG::Vec3d>(propertyInfo2.propName());
177 
178  VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
179 
180  if (w->vecFieldDiff_4symm_rb->isChecked())
181  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfFaceIter, OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halffaces_begin(), mesh->halffaces_end());
182  if (w->vecFieldDiff_diff_norm_rb->isChecked())
183  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfFaceIter, OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halffaces_begin(), mesh->halffaces_end());
184  if (w->vecFieldDiff_norm_diff_rb->isChecked())
185  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfFaceIter, OpenVolumeMesh::HalfFacePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halffaces_begin(), mesh->halffaces_end());
186  if (_setDrawMode)
187  {
188  VolumeMeshObject<MeshT>* object;
189  PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
190  object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.halffacesColoredPerHalfface);
191  }
192 }
193 
194 template <typename MeshT>
195 void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeEdgeProp(bool _setDrawMode)
196 {
197  MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
198 
199  OpenVolumeMesh::EdgePropertyT<ACG::Vec3d> prop1 = mesh->template request_edge_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
200  OpenVolumeMesh::EdgePropertyT<ACG::Vec3d> prop2 = mesh->template request_edge_property<ACG::Vec3d>(propertyInfo2.propName());
201 
202  VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
203 
204  if (w->vecFieldDiff_4symm_rb->isChecked())
205  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::EdgeIter, OpenVolumeMesh::EdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->edges_begin(), mesh->edges_end());
206  if (w->vecFieldDiff_diff_norm_rb->isChecked())
207  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::EdgeIter, OpenVolumeMesh::EdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->edges_begin(), mesh->edges_end());
208  if (w->vecFieldDiff_norm_diff_rb->isChecked())
209  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::EdgeIter, OpenVolumeMesh::EdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->edges_begin(), mesh->edges_end());
210  if (_setDrawMode)
211  {
212  VolumeMeshObject<MeshT>* object;
213  PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
214  object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.edgesColoredPerEdge);
215  }
216 }
217 
218 template <typename MeshT>
219 void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeHalfedgeProp(bool _setDrawMode)
220 {
221  MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
222 
223  OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d> prop1 = mesh->template request_halfedge_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
224  OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d> prop2 = mesh->template request_halfedge_property<ACG::Vec3d>(propertyInfo2.propName());
225 
226  VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
227 
228  if (w->vecFieldDiff_4symm_rb->isChecked())
229  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfEdgeIter, OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halfedges_begin(), mesh->halfedges_end());
230  if (w->vecFieldDiff_diff_norm_rb->isChecked())
231  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfEdgeIter, OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halfedges_begin(), mesh->halfedges_end());
232  if (w->vecFieldDiff_norm_diff_rb->isChecked())
233  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::HalfEdgeIter, OpenVolumeMesh::HalfEdgePropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->halfedges_begin(), mesh->halfedges_end());
234  if (_setDrawMode)
235  {
236  VolumeMeshObject<MeshT>* object;
237  PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
238  object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.halfedgesColoredPerHalfedge);
239  }
240 }
241 
242 template <typename MeshT>
243 void OVMPropertyVisualizerVectorFieldDifference<MeshT>::visualizeVertexProp(bool _setDrawMode)
244 {
245  MeshT* mesh = OVMPropertyVisualizer<MeshT>::mesh;
246 
247  OpenVolumeMesh::VertexPropertyT<ACG::Vec3d> prop1 = mesh->template request_vertex_property<ACG::Vec3d>(PropertyVisualizer::propertyInfo.propName());
248  OpenVolumeMesh::VertexPropertyT<ACG::Vec3d> prop2 = mesh->template request_vertex_property<ACG::Vec3d>(propertyInfo2.propName());
249 
250  VectorFieldDifferenceWidget* w = static_cast<VectorFieldDifferenceWidget*>(PropertyVisualizer::widget);
251 
252  if (w->vecFieldDiff_4symm_rb->isChecked())
253  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::VertexIter, OpenVolumeMesh::VertexPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_4_symm_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->vertices_begin(), mesh->vertices_end());
254  if (w->vecFieldDiff_diff_norm_rb->isChecked())
255  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::VertexIter, OpenVolumeMesh::VertexPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_diff_of_norms<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->vertices_begin(), mesh->vertices_end());
256  if (w->vecFieldDiff_norm_diff_rb->isChecked())
257  OVMPVVFD::colorElements< MeshT, OpenVolumeMesh::VertexIter, OpenVolumeMesh::VertexPropertyT<ACG::Vec3d>, OVMPVVFD::scalarFn_norm_of_diff<MeshT> >(OVMPropertyVisualizer<MeshT>::mObjectID, prop1, prop2, mesh->vertices_begin(), mesh->vertices_end());
258  if (_setDrawMode)
259  {
260  VolumeMeshObject<MeshT>* object;
261  PluginFunctions::getObject(OVMPropertyVisualizer<MeshT>::mObjectID, object);
262  object->setObjectDrawMode(OVMPropertyVisualizer<MeshT>::drawModes.verticesColored);
263  }
264 }
265 
PluginFunctions::getObject
bool getObject(const int _identifier, BaseObject *&_object)
Get the object which has the given identifier.
Definition: PluginFunctions.cc:298
PropertyInfo
Cellection of information about a property.
Definition: Utils.hh:109
OpenMesh::VectorT::norm
auto norm() const -> decltype(std::sqrt(std::declval< VectorT< S, DIM >>().sqrnorm()))
compute euclidean norm
Definition: Vector11T.hh:433
ACG::ColorCoder
Class for generating nice colors for doubles.
Definition: ColorCoder.hh:68