← Back to branch summary

~ubuntu-branches/ubuntu/saucy/rheolef/saucy

« back to all changes in this revision

Viewing changes to skit/ptst2/vec_dot_tst.cc

  • Committer: Bazaar Package Importer
  • Author(s): Pierre Saramito
  • Date: 2011-02-08 10:00:06 UTC
  • mfrom: (1.1.3 upstream)
  • Revision ID: james.westby@ubuntu.com-20110208100006-jg8rarizvhdm8e82
Tags: 5.92-1
http://bugs.debian.org/607117
http://bugs.debian.org/607184
http://bugs.debian.org/607038
* New upstream release:
  - "rheolef" suffix added to all unix manuals (Closes: #607117)
  - minor source code fixes for g++-4.5 (Closes: #607184)
  - minor reference manual fix (Closes: #607038)
* debian/control: libboost-iostreams-dev dependency added, as required
  by the new upstream version.
 

Show diffs side-by-side

added added

removed removed

Lines of Context:
skit/ptst2/vec_dot_tst.cc
 
1
///
 
2
/// This file is part of Rheolef.
 
3
///
 
4
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
 
5
///
 
6
/// Rheolef is free software; you can redistribute it and/or modify
 
7
/// it under the terms of the GNU General Public License as published by
 
8
/// the Free Software Foundation; either version 2 of the License, or
 
9
/// (at your option) any later version.
 
10
///
 
11
/// Rheolef is distributed in the hope that it will be useful,
 
12
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
14
/// GNU General Public License for more details.
 
15
///
 
16
/// You should have received a copy of the GNU General Public License
 
17
/// along with Rheolef; if not, write to the Free Software
 
18
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 
19
///
 
20
/// =========================================================================
 
21
# include "rheolef/array.h"
 
22
# include "rheolef/parstream.h"
 
23
# include "rheolef/par_inner_product.h"
 
24
using namespace rheolef;
 
25
using namespace std;
 
26
int main(int argc, char**argv)
 
27
{
 
28
    environment parallel(argc, argv);
 
29
    array<Float> x(100, Float(2.0)/Float(3.0));
 
30
    Float norme2_x = par_inner_product (x.begin(),x.end(),x.begin(),x.comm()); 
 
31
    int ndigit;
 
32
    switch(sizeof(Float)) {
 
33
      case 16: ndigit = 20;
 
34
      case 8: 
 
35
      default: ndigit = 10;
 
36
    };
 
37
    cout.precision(ndigit);
 
38
    pcout << "dot(x,x) = " << norme2_x << endl;
 
39
    return 0;
 
40
}