blob: 3a0de102edd2ff10da6696934a750e1bc4a81927 (
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
|
// Licensed to the .NET Foundation under one or more agreements.
// The .NET Foundation licenses this file to you under the MIT license.
// See the LICENSE file in the project root for more information.
//
// Integration by corrected trapezoid rule adapted from Conte and de Boor
using Microsoft.Xunit.Performance;
using System;
using System.Runtime.CompilerServices;
using Xunit;
[assembly: OptimizeForBenchmarks]
[assembly: MeasureInstructionsRetired]
namespace Benchstone.BenchF
{
public static class Trap
{
#if DEBUG
public const int Iterations = 1;
#else
public const int Iterations = 240000;
#endif
[MethodImpl(MethodImplOptions.NoInlining)]
private static bool Bench()
{
int nm1, idbg;
double t2, cortrp, trap, a, b, h;
trap = 0.0;
cortrp = 0.0;
idbg = 0;
for (int j = 1; j <= Iterations; j++)
{
a = 0;
b = 1;
if (idbg != 0)
{
System.Console.WriteLine("trapazoid sum corr.trap sum \n");
}
for (int n = 10; n <= 15; n++)
{
h = (b - a) / n;
nm1 = n - 1;
trap = (F(a) + F(b)) / 2;
for (int i = 1; i <= nm1; i++)
{
t2 = a + i * h;
trap = trap + F(t2);
}
trap = trap * h;
cortrp = trap + h * h * (FPrime(a) - FPrime(b)) / 12;
if (idbg != 0)
{
System.Console.WriteLine("{0}, {1}, {2}\n", n, trap, cortrp);
}
}
}
return true;
}
private static double F(double x)
{
return (System.Math.Exp(-(x) * (x)));
}
private static double FPrime(double x)
{
return ((-2) * (x) * (F(x)));
}
[Benchmark]
public static void Test()
{
foreach (var iteration in Benchmark.Iterations)
{
using (iteration.StartMeasurement())
{
Bench();
}
}
}
private static bool TestBase()
{
bool result = Bench();
return result;
}
public static int Main()
{
bool result = TestBase();
return (result ? 100 : -1);
}
}
}
|