uo9xpb74
/
mathnet-numerics
forked from tsai/mathnet-numerics
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity
Browse Source

remove Control.ParallelizeOperation()

provider
Thomas Ibel 13 years ago
parent
8f0a48354c
commit
0da254f1cc
5 changed files with 120 additions and 430 deletions
Whitespace
Show all changes Ignore whitespace when comparing lines Ignore changes in amount of whitespace Ignore changes in whitespace at EOL
Split View
Diff Options
Show Stats Download Patch File Download Diff File
  1. 10
      src/Numerics/Control.cs
  2. 144
      src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Complex.cs
  3. 144
      src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Complex32.cs
  4. 126
      src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Double.cs
  5. 126
      src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Single.cs

10
src/Numerics/Control.cs
View File

@ -183,16 +183,6 @@ namespace MathNet.Numerics
set { _parallelizeElements = Math.Max(3, value); }
}
/// <summary>
/// Given the number elements, should the operation be parallelized.
/// </summary>
/// <param name="elements">The number elements to check.</param>
/// <returns><c>true</c> if the operation should be parallelized; <c>false</c> otherwise.</returns>
public static bool ParallelizeOperation(int elements)
{
return NumberOfParallelWorkerThreads >= 2 && elements >= ParallelizeElements;
}
/// <summary>
/// Maximum number of columns to print in ToString methods by default.
/// </summary>

144
src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Complex.cs
View File

@ -82,43 +82,23 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else if (alpha.IsOne())
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + x[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + x[index];
result[i] = y[i] + x[i];
}
}
});
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + (alpha*x[i]);
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + (alpha*x[index]);
result[i] = y[i] + (alpha*x[i]);
}
}
});
}
}
@ -146,23 +126,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = alpha*x[i];
}
});
}
else
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = alpha*x[index];
result[i] = alpha*x[i];
}
}
});
}
}
@ -178,23 +148,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentNullException("x");
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i].Conjugate();
}
});
}
else
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index].Conjugate();
result[i] = x[i].Conjugate();
}
}
});
}
/// <summary>
@ -262,23 +222,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] + y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] + y[index];
result[i] = x[i] + y[i];
}
}
});
}
/// <summary>
@ -313,23 +263,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] - y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] - y[index];
result[i] = x[i] - y[i];
}
}
});
}
/// <summary>
@ -364,23 +304,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]*y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]*y[index];
result[i] = x[i] * y[i];
}
}
});
}
/// <summary>
@ -415,23 +345,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]/y[i];
}
});
}
else
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]/y[index];
result[i] = x[i] / y[i];
}
}
});
}
/// <summary>

144
src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Complex32.cs
View File

@ -78,43 +78,23 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else if (alpha.IsOne())
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + x[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + x[index];
result[i] = y[i] + x[i];
}
}
});
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + (alpha*x[i]);
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + (alpha*x[index]);
result[i] = y[i] + (alpha * x[i]);
}
}
});
}
}
@ -143,23 +123,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = alpha*x[i];
}
});
}
else
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = alpha*x[index];
result[i] = alpha * x[i];
}
}
});
}
}
@ -175,23 +145,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentNullException("x");
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i].Conjugate();
}
});
}
else
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index].Conjugate();
result[i] = x[i].Conjugate();
}
}
});
}
/// <summary>
@ -260,23 +220,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] + y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] + y[index];
result[i] = x[i] + y[i];
}
}
});
}
/// <summary>
@ -311,23 +261,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] - y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] - y[index];
result[i] = x[i] - y[i];
}
}
});
}
/// <summary>
@ -362,23 +302,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]*y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]*y[index];
result[i] = x[i] * y[i];
}
}
});
}
/// <summary>
@ -413,23 +343,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]/y[i];
}
});
}
else
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]/y[index];
result[i] = x[i] / y[i];
}
}
});
}
/// <summary>

126
src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Double.cs
View File

@ -76,43 +76,23 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else if (alpha == 1.0)
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + x[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + x[index];
result[i] = y[i] + x[i];
}
}
});
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + (alpha*x[i]);
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + (alpha*x[index]);
result[i] = y[i] + (alpha * x[i]);
}
}
});
}
}
@ -140,23 +120,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = alpha*x[i];
}
});
}
else
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = alpha*x[index];
result[i] = alpha * x[i];
}
}
});
}
}
@ -244,23 +214,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] + y[i];
}
});
}
else
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] + y[index];
result[i] = x[i] + y[i];
}
}
});
}
/// <summary>
@ -295,23 +255,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] - y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] - y[index];
result[i] = x[i] - y[i];
}
}
});
}
/// <summary>
@ -346,23 +296,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]*y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]*y[index];
result[i] = x[i] * y[i];
}
}
});
}
/// <summary>
@ -397,23 +337,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]/y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]/y[index];
result[i] = x[i] / y[i];
}
}
});
}
/// <summary>

126
src/Numerics/Providers/LinearAlgebra/ManagedLinearAlgebraProvider.Single.cs
View File

@ -76,43 +76,23 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else if (alpha == 1.0)
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + x[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + x[index];
result[i] = y[i] + x[i];
}
}
});
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = y[i] + (alpha*x[i]);
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = y[index] + (alpha*x[index]);
result[i] = y[i] + (alpha * x[i]);
}
}
});
}
}
@ -140,23 +120,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
}
else
{
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = alpha*x[i];
}
});
}
else
CommonParallel.For(0, x.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = alpha*x[index];
result[i] = alpha * x[i];
}
}
});
}
}
@ -244,23 +214,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] + y[i];
}
});
}
else
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] + y[index];
result[i] = x[i] + y[i];
}
}
});
}
/// <summary>
@ -295,23 +255,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i] - y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index] - y[index];
result[i] = x[i] - y[i];
}
}
});
}
/// <summary>
@ -346,23 +296,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]*y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]*y[index];
result[i] = x[i] * y[i];
}
}
});
}
/// <summary>
@ -397,23 +337,13 @@ namespace MathNet.Numerics.Providers.LinearAlgebra
throw new ArgumentException(Resources.ArgumentVectorsSameLength);
}
if (Control.ParallelizeOperation(x.Length))
{
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (int i = a; i < b; i++)
{
result[i] = x[i]/y[i];
}
});
}
else
CommonParallel.For(0, y.Length, 4096, (a, b) =>
{
for (var index = 0; index < x.Length; index++)
for (int i = a; i < b; i++)
{
result[index] = x[index]/y[index];
result[i] = x[i] / y[i];
}
}
});
}
/// <summary>

Write Preview
Loading…
Cancel
Save