blob: b90127b87434e0fe620655db5ad84c1093827805 (
plain)
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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
|
// Copyright (c) Microsoft. All rights reserved.
// Licensed under the MIT license. See LICENSE file in the project root for full license information.
//
using Microsoft.Xunit.Performance;
using System;
using System.Runtime.CompilerServices;
using Xunit;
[assembly: OptimizeForBenchmarks]
[assembly: MeasureInstructionsRetired]
public static class DMath
{
#if DEBUG
public const int Iterations = 1;
#else
public const int Iterations = 100000;
#endif
const double Deg2Rad = 57.29577951;
static volatile object VolatileObject;
static void Escape(object obj) {
VolatileObject = obj;
}
static double Fact(double n)
{
double res;
res = 1.0;
while (n > 0.0) {
res *= n;
n -= 1.0;
}
return res;
}
static double Power(double n, double p)
{
double res;
res = 1.0;
while (p > 0.0) {
res *= n;
p -= 1.0;
}
return res;
}
[MethodImpl(MethodImplOptions.NoInlining)]
static bool Bench(int loop)
{
double[] sines = new double[91];
double angle, radians, sine, worksine, temp, k;
double diff;
for (int iter = 1; iter <= loop; iter++) {
for (angle = 0.0; angle <= 90.0; angle += 1.0) {
radians = angle / Deg2Rad;
k = 0.0;
worksine = 0.0;
do {
sine = worksine;
temp = (2.0 * k) + 1.0;
worksine += (Power(-1.0, k) / Fact(temp)) * Power(radians, temp);
k += 1.0;
diff = Math.Abs(sine - worksine);
} while (diff > 1E-8);
sines[(int)angle] = worksine;
}
}
// Escape sines array so that its elements appear live-out
Escape(sines);
return true;
}
[Benchmark]
public static void Test() {
foreach (var iteration in Benchmark.Iterations) {
using (iteration.StartMeasurement()) {
Bench(Iterations);
}
}
}
static bool TestBase() {
bool result = Bench(Iterations);
return result;
}
public static int Main() {
bool result = TestBase();
return (result ? 100 : -1);
}
}
|