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
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
|
// 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.
// Adapted from mandelbrot C# .NET Core #4 program
// http://benchmarksgame.alioth.debian.org/u64q/program.php?test=mandelbrot&lang=csharpcore&id=4
// Best-scoring C# .NET Core version as of 2017-09-01
/* The Computer Language Benchmarks Game
http://benchmarksgame.alioth.debian.org/
started with Java #2 program (Krause/Whipkey/Bennet/AhnTran/Enotus/Stalcup)
adapted for C# by Jan de Vaan
simplified and optimised to use TPL by Anthony Lloyd
*/
using System;
using System.Threading.Tasks;
using System.IO;
using System.Runtime.CompilerServices;
using System.Security.Cryptography;
using Microsoft.Xunit.Performance;
using Xunit;
[assembly: OptimizeForBenchmarks]
namespace BenchmarksGame
{
public class MandelBrot_4
{
[MethodImpl(MethodImplOptions.AggressiveInlining)]
static byte getByte(double[] Crb, double Ciby, int x, int y)
{
int res = 0;
for (int i = 0; i < 8; i += 2)
{
double Crbx = Crb[x + i], Crbx1 = Crb[x + i + 1];
double Zr1 = Crbx, Zr2 = Crbx1;
double Zi1 = Ciby, Zi2 = Ciby;
int b = 0;
int j = 49;
do
{
double nZr1 = Zr1 * Zr1 - Zi1 * Zi1 + Crbx;
Zi1 = Zr1 * Zi1 + Zr1 * Zi1 + Ciby;
Zr1 = nZr1;
double nZr2 = Zr2 * Zr2 - Zi2 * Zi2 + Crbx1;
Zi2 = Zr2 * Zi2 + Zr2 * Zi2 + Ciby;
Zr2 = nZr2;
if (Zr1 * Zr1 + Zi1 * Zi1 > 4) { b |= 2; if (b == 3) break; }
if (Zr2 * Zr2 + Zi2 * Zi2 > 4) { b |= 1; if (b == 3) break; }
} while (--j > 0);
res = (res << 2) + b;
}
return (byte)(res ^ -1);
}
public static int Main(String[] args)
{
var n = args.Length > 0 ? Int32.Parse(args[0]) : 80;
var data = DoBench(n);
Console.Out.WriteLine("P4\n{0} {0}", n);
Console.OpenStandardOutput().Write(data, 0, data.Length);
if (!MatchesChecksum(data, "3B-EF-65-05-1D-39-7F-9B-96-8D-EF-98-BF-06-CE-74"))
{
return -1;
}
return 100;
}
// Commented out data left in source to provide checksums for each case
[Benchmark(InnerIterationCount = 7)]
//[InlineData(1000, "B2-13-51-CE-B0-29-2C-4E-75-5E-91-19-18-E4-0C-D9")]
//[InlineData(2000, "5A-21-55-9B-7B-18-2F-34-9B-33-C5-F9-B5-2C-40-56")]
//[InlineData(3000, "E5-82-85-0A-3C-89-69-B1-A8-21-63-52-75-B3-C8-33")]
[InlineData(4000, "C7-E6-66-43-66-73-F8-A8-D3-B4-D7-97-2F-FC-A1-D3")]
//[InlineData(5000, "6D-36-F1-F6-37-8F-34-EB-52-F9-2D-11-89-12-B2-2F")]
//[InlineData(6000, "8B-05-78-EB-2E-0E-98-F2-C7-39-76-ED-0F-A9-D2-B8")]
//[InlineData(7000, "01-F8-F2-2A-AB-70-C7-BA-E3-64-19-E7-D2-84-DF-57")]
//[InlineData(8000, "C8-ED-D7-FB-65-66-3A-D9-C6-04-9E-96-E8-CA-4F-2C")]
public static void Bench(int n, string checksum)
{
byte[] bytes = null;
Benchmark.Iterate(() =>
{
bytes = DoBench(n);
});
Assert.True(MatchesChecksum(bytes, checksum));
}
static bool MatchesChecksum(byte[] bytes, string checksum)
{
using (var md5 = MD5.Create())
{
byte[] hash = md5.ComputeHash(bytes);
return (checksum == BitConverter.ToString(hash));
}
}
static byte[] DoBench(int n)
{
double invN = 2.0 / n;
var Crb = new double[n + 7];
for (int i = 0; i < n; i++) { Crb[i] = i * invN - 1.5; }
int lineLen = (n - 1) / 8 + 1;
var data = new byte[n * lineLen];
Parallel.For(0, n, y =>
{
var Ciby = y * invN - 1.0;
var offset = y * lineLen;
for (int x = 0; x < lineLen; x++)
data[offset + x] = getByte(Crb, Ciby, x * 8, y);
});
return data;
}
}
}
|