~ubuntu-branches/ubuntu/trusty/rheolef/trusty-proposed

« back to all changes in this revision

Viewing changes to nfem/lib/compute-bbox.h

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme
  • Date: 2010-06-12 09:08:59 UTC
  • Revision ID: james.westby@ubuntu.com-20100612090859-8gpm2gc7j3ab43et
Tags: upstream-5.89
ImportĀ upstreamĀ versionĀ 5.89

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
#ifndef _RHEO_COMPUTE_BBOX_H
 
2
#define _RHEO_COMPUTE_BBOX_H
 
3
///
 
4
/// This file is part of Rheolef.
 
5
///
 
6
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
 
7
///
 
8
/// Rheolef is free software; you can redistribute it and/or modify
 
9
/// it under the terms of the GNU General Public License as published by
 
10
/// the Free Software Foundation; either version 2 of the License, or
 
11
/// (at your option) any later version.
 
12
///
 
13
/// Rheolef is distributed in the hope that it will be useful,
 
14
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
 
15
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
16
/// GNU General Public License for more details.
 
17
///
 
18
/// You should have received a copy of the GNU General Public License
 
19
/// along with Rheolef; if not, write to the Free Software
 
20
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
21
/// 
 
22
/// =========================================================================
 
23
//
 
24
// used by cad and geo classes
 
25
//
 
26
template <
 
27
        class PointIterator,
 
28
        class Point,
 
29
        class Size>
 
30
void
 
31
compute_bbox (
 
32
        PointIterator iter_p,
 
33
        PointIterator last_p,
 
34
        Size          dim,
 
35
        Point&        xmin,
 
36
        Point&        xmax)
 
37
{
 
38
  typedef typename Point::float_type                          float_type;
 
39
  typedef typename iterator_traits<PointIterator>::value_type Point1;
 
40
 
 
41
  for (Size j = 0; j < dim; j++) {
 
42
    xmin[j] =  std::numeric_limits<float_type>::max();
 
43
    xmax[j] = -std::numeric_limits<float_type>::max();
 
44
  }
 
45
  for (Size j = dim+1; j < 3; j++) {
 
46
    xmin [j] = xmax [j] = 0;
 
47
  }
 
48
  while (iter_p != last_p) {
 
49
    const Point1& p = *iter_p++;
 
50
    for (Size j = 0 ; j < dim; j++) {
 
51
      xmin[j] = min(p[j], xmin[j]);
 
52
      xmax[j] = max(p[j], xmax[j]);
 
53
    }
 
54
  }
 
55
}
 
56
#endif // _RHEO_BBOX_H