// // Math.NET Numerics, part of the Math.NET Project // http://numerics.mathdotnet.com // http://github.com/mathnet/mathnet-numerics // // Copyright (c) 2009-2018 Math.NET // // Permission is hereby granted, free of charge, to any person // obtaining a copy of this software and associated documentation // files (the "Software"), to deal in the Software without // restriction, including without limitation the rights to use, // copy, modify, merge, publish, distribute, sublicense, and/or sell // copies of the Software, and to permit persons to whom the // Software is furnished to do so, subject to the following // conditions: // // The above copyright notice and this permission notice shall be // included in all copies or substantial portions of the Software. // // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, // EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES // OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND // NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT // HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, // WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING // FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR // OTHER DEALINGS IN THE SOFTWARE. // #if NATIVE using System; using IStation.Numerics.Providers.Common.Mkl; namespace IStation.Numerics.Providers.LinearAlgebra.Mkl { /// /// Error codes return from the MKL provider. /// internal enum MklError : int { /// /// Unable to allocate memory. /// MemoryAllocation = -999999 } /// /// Intel's Math Kernel Library (MKL) linear algebra provider. /// internal partial class MklLinearAlgebraProvider : Managed.ManagedLinearAlgebraProvider, IDisposable { const int MinimumCompatibleRevision = 4; readonly string _hintPath; readonly MklConsistency _consistency; readonly MklPrecision _precision; readonly MklAccuracy _accuracy; int _linearAlgebraMajor; int _linearAlgebraMinor; int _vectorFunctionsMajor; int _vectorFunctionsMinor; /// Hint path where to look for the native binaries /// /// Sets the desired bit consistency on repeated identical computations on varying CPU architectures, /// as a trade-off with performance. /// /// VML optimal precision and rounding. /// VML accuracy mode. internal MklLinearAlgebraProvider(string hintPath, MklConsistency consistency, MklPrecision precision, MklAccuracy accuracy) { _hintPath = hintPath; _consistency = consistency; _precision = precision; _accuracy = accuracy; } /// /// Try to find out whether the provider is available, at least in principle. /// Verification may still fail if available, but it will certainly fail if unavailable. /// public override bool IsAvailable() { return MklProvider.IsAvailable(hintPath: _hintPath); } /// /// Initialize and verify that the provided is indeed available. /// If calling this method fails, consider to fall back to alternatives like the managed provider. /// public override void InitializeVerify() { int revision = MklProvider.Load(_hintPath, _consistency, _precision, _accuracy); if (revision < MinimumCompatibleRevision) { throw new NotSupportedException(FormattableString.Invariant($"MKL Native Provider revision r{revision} is too old. Consider upgrading to a newer version. Revision r{MinimumCompatibleRevision} and newer are supported.")); } _linearAlgebraMajor = SafeNativeMethods.query_capability((int)ProviderCapability.LinearAlgebraMajor); _linearAlgebraMinor = SafeNativeMethods.query_capability((int)ProviderCapability.LinearAlgebraMinor); _vectorFunctionsMajor = SafeNativeMethods.query_capability((int)ProviderCapability.VectorFunctionsMajor); _vectorFunctionsMinor = SafeNativeMethods.query_capability((int)ProviderCapability.VectorFunctionsMinor); // we only support exactly one major version, since major version changes imply a breaking change. if (_linearAlgebraMajor != 2) { throw new NotSupportedException(FormattableString.Invariant($"MKL Native Provider not compatible. Expecting linear algebra v2 but provider implements v{_linearAlgebraMajor}.")); } } /// /// Frees memory buffers, caches and handles allocated in or to the provider. /// Does not unload the provider itself, it is still usable afterwards. /// public override void FreeResources() { MklProvider.FreeResources(); } public override string ToString() { return MklProvider.Describe(); } public void Dispose() { FreeResources(); } } } #endif