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
|
// Copyright 2017 Global Phasing Ltd.
// Microbenchmark of gemmi::Residue::has_standard_pdb_name().
// Requires the google/benchmark library. It can be built manually:
// c++ -Wall -O2 -I../include -I$GB/include resinfo.cpp $GB/src/libbenchmark.a -pthread
#include "gemmi/resinfo.hpp"
#include <benchmark/benchmark.h>
#include <stdlib.h> // for rand
static const std::string residue_names[10] =
{ "GLY", "ASN", "ASN", "HOH", "LEU", "DC", "G", "GLU", "ALA", "SER" };
static void find_tabulated_residue_x10(benchmark::State& state) {
std::string names[10];
for (int i = 0; i != 10; ++i)
names[i] = rand() % 1000 == 0 ? "GLN" : residue_names[i];
while (state.KeepRunning()) {
int nh = 0;
for (int i = 0; i != 10; ++i)
nh += gemmi::find_tabulated_residue(names[i]).hydrogen_count;
benchmark::DoNotOptimize(nh);
}
}
BENCHMARK(find_tabulated_residue_x10);
BENCHMARK_MAIN();
/* Output from my laptop:
------------------------------------------------------------------
Benchmark Time CPU Iterations
------------------------------------------------------------------
find_tabulated_residue_x10 52 ns 52 ns 13458745
*/
|
