File: timeLinearize.h

package info (click to toggle)
gtsam 4.2.0%2Bdfsg-3
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 46,108 kB
  • sloc: cpp: 127,191; python: 14,312; xml: 8,442; makefile: 252; sh: 119; ansic: 101
file content (56 lines) | stat: -rw-r--r-- 1,740 bytes parent folder | download | duplicates (2)
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
/* ----------------------------------------------------------------------------

 * GTSAM Copyright 2010, Georgia Tech Research Corporation,
 * Atlanta, Georgia 30332-0415
 * All Rights Reserved
 * Authors: Frank Dellaert, et al. (see THANKS for the full author list)

 * See LICENSE for the license information

 * -------------------------------------------------------------------------- */

/**
 * @file    timeLinearize.h
 * @brief   time linearize
 * @author  Frank Dellaert
 * @date    October 10, 2014
 */

#pragma once

#include <gtsam/nonlinear/NonlinearFactorGraph.h>
#include <gtsam/linear/GaussianFactorGraph.h>

#include <time.h>
#include <iostream>
#include <iomanip>

static const int n = 1000000;

void timeSingleThreaded(const std::string& str,
    const gtsam::NonlinearFactor::shared_ptr& f, const gtsam::Values& values) {
  long timeLog = clock();
  gtsam::GaussianFactor::shared_ptr gf;
  for (int i = 0; i < n; i++)
    gf = f->linearize(values);
  long timeLog2 = clock();
  double seconds = (double) (timeLog2 - timeLog) / CLOCKS_PER_SEC;
  std::cout << std::setprecision(3);
  std::cout << str << ((double) seconds * 1000000 / n) << " musecs/call"
      << std::endl;
}

void timeMultiThreaded(const std::string& str,
    const gtsam::NonlinearFactor::shared_ptr& f, const gtsam::Values& values) {
  gtsam::NonlinearFactorGraph graph;
  for (int i = 0; i < n; i++)
    graph.push_back(f);
  long timeLog = clock();
  gtsam::GaussianFactorGraph::shared_ptr gfg = graph.linearize(values);
  long timeLog2 = clock();
  double seconds = (double) (timeLog2 - timeLog) / CLOCKS_PER_SEC;
  std::cout << std::setprecision(3);
  std::cout << str << ((double) seconds * 1000000 / n) << " musecs/call"
      << std::endl;
}