Class NativeOpenCLDenseLinearAlgebraDoubleBackend<E extends FieldElement<E>>

java.lang.Object

org.episteme.nativ.mathematics.linearalgebra.backends.NativeOpenCLDenseLinearAlgebraDoubleBackend<E>

All Implemented Interfaces:

AutoCloseable, LinearAlgebraProvider<E>, AlgorithmProvider, Backend, ComputeBackend, GPUBackend, NativeBackend

@AutoService({Backend.class,ComputeBackend.class,NativeBackend.class,LinearAlgebraProvider.class,GPUBackend.class})
public class NativeOpenCLDenseLinearAlgebraDoubleBackend<E extends FieldElement<E>>
extends Object
implements LinearAlgebraProvider<E>, NativeBackend, GPUBackend

OpenCL implementation of Dense Linear Algebra Provider for Double precision. Requires fp64 support on the GPU device.

Nested Class Summary

Nested classes/interfaces inherited from interface GPUBackend

GPUBackend.DeviceInfo

Constructor Summary

Constructors

Constructor	Description
NativeOpenCLDenseLinearAlgebraDoubleBackend()

Method Summary

Modifier and Type	Method	Description
Vector<E>	add(Vector<E> a, Vector<E> b)
long	allocateGPUMemory(long sizeBytes)	Allocates GPU memory and returns a handle.
E	angle(Vector<E> a, Vector<E> b)	Returns the angle between two vectors in radians.
CholeskyResult<E>	cholesky(Matrix<E> a)	Computes the Cholesky decomposition of the specified matrix.
void	close()
Matrix<E>	conjugateTranspose(Matrix<E> a)
void	copyFromGPU(long gpuHandle, DoubleBuffer hostBuffer, long sizeBytes)	Copies data from GPU to host.
void	copyToGPU(long gpuHandle, DoubleBuffer hostBuffer, long sizeBytes)	Copies data from host to GPU.
ExecutionContext	createContext()	Creates an execution context for running operations.
Vector<E>	cross(Vector<E> a, Vector<E> b)	Returns the cross product of two 3D vectors.
E	determinant(Matrix<E> a)
E	dot(Vector<E> a, Vector<E> b)
EigenResult<E>	eigen(Matrix<E> a)	Computes the eigenvalue decomposition of the specified matrix.
void	freeGPUMemory(long gpuHandle)	Frees GPU memory.
HardwareAccelerator	getAcceleratorType()	Returns the type of hardware accelerator used by this backend.
GPUBackend.DeviceInfo[]	getDevices()	Returns information about available GPU devices.
String	getId()	Returns the unique identifier for this backend (e.g., "javafx", "jmol", "cuda").
String	getName()	Returns the display name for UI presentation.
String	getNativeLibraryName()	Returns the name of the native library this backend depends on.
int	getPriority()	Priority of this provider (higher means more preferred).
String	getType()	Returns the backend type category (e.g., "plotting", "molecular", "tensor", "math").
Matrix<E>	inverse(Matrix<E> a)
boolean	isAvailable()	Checks if the provider is available in the current environment.
boolean	isCompatible(Ring<?> ring)	Checks if this provider is compatible with the given ring.
boolean	isExplicitlyDisabled()	Checks if this backend is specifically disabled via a system property: "episteme.backend.disable." + getId().
boolean	isLoaded()	Checks if the native library is loaded and available.
LUResult<E>	lu(Matrix<E> a)	Computes the LU decomposition of the specified matrix.
void	matrixMultiply(DoubleBuffer A, DoubleBuffer B, DoubleBuffer C, int m, int n, int k)	Performs matrix multiplication C = A * B on the GPU.
Matrix<E>	multiply(Matrix<E> a, Matrix<E> b)
Vector<E>	multiply(Matrix<E> a, Vector<E> b)
Vector<E>	multiply(Vector<E> v, E s)
E	norm(Vector<E> a)
Vector<E>	normalize(Vector<E> v)	Returns the normalized vector (unit vector).
QRResult<E>	qr(Matrix<E> a)	Computes the QR decomposition of the specified matrix.
Matrix<E>	scale(E scalar, Matrix<E> a)
void	selectDevice(int deviceId)	Selects a GPU device for computation.
void	shutdown()	Called when the provider is no longer needed (e.g., application shutdown).
Vector<E>	solve(Matrix<E> a, Vector<E> b)
Vector<E>	solveTriangular(Matrix<E> a, Vector<E> b, boolean upper, boolean transpose, boolean conjugate, boolean unit)	Solves the triangular system Ax = b.
Matrix<E>	subtract(Matrix<E> a, Matrix<E> b)
Vector<E>	subtract(Vector<E> a, Vector<E> b)
SVDResult<E>	svd(Matrix<E> a)	Computes the Singular Value Decomposition (SVD) of the specified matrix.
void	synchronize()	Synchronizes GPU execution (waits for all kernels to complete).
E	trace(Matrix<E> a)
Matrix<E>	transpose(Matrix<E> a)

Methods inherited from class Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

Methods inherited from interface AlgorithmProvider

description, getMetadata

Methods inherited from interface Backend

getAlgorithmProviders, getStatusMessage

Methods inherited from interface ComputeBackend

createBackend, getDescription, supportsComplexNumbers, supportsFloatingPoint, supportsParallelOps

Methods inherited from interface LinearAlgebraProvider

acos, acosh, add, asin, asinh, atan, atanh, cbrt, conditionNumber, configure, cos, cosh, exp, getAlgorithmType, getEnvironmentInfo, log, log10, pow, projection, rank, score, sin, sinh, solve, solve, solve, sqrt, tan, tanh

Constructor Details

NativeOpenCLDenseLinearAlgebraDoubleBackend

public NativeOpenCLDenseLinearAlgebraDoubleBackend()

Method Details

isLoaded

public boolean isLoaded()

Description copied from interface: NativeBackend

Checks if the native library is loaded and available.

Specified by:

isLoaded in interface NativeBackend

Returns:

true if the native library has been successfully loaded

getNativeLibraryName

public String getNativeLibraryName()

Description copied from interface: NativeBackend

Returns the name of the native library this backend depends on.

Specified by:

getNativeLibraryName in interface NativeBackend

Returns:

library name (e.g., "openblas", "cuda", "fftw3")

isAvailable

public boolean isAvailable()

Description copied from interface: AlgorithmProvider

Checks if the provider is available in the current environment.

Specified by:

isAvailable in interface AlgorithmProvider

Specified by:

isAvailable in interface Backend

isExplicitlyDisabled

public boolean isExplicitlyDisabled()

Description copied from interface: Backend

Checks if this backend is specifically disabled via a system property: "episteme.backend.disable." + getId().

It also checks specialized properties like "episteme.linearalgebra.disable." + getId().

Specified by:

isExplicitlyDisabled in interface Backend

Returns:

true if explicitly disabled

isCompatible

public boolean isCompatible(Ring<?> ring)

Description copied from interface: LinearAlgebraProvider

Checks if this provider is compatible with the given ring.

Specified by:

isCompatible in interface LinearAlgebraProvider<E extends FieldElement<E>>

getId

public String getId()

Description copied from interface: Backend

Returns the unique identifier for this backend (e.g., "javafx", "jmol", "cuda").

Specified by:

getId in interface Backend

Specified by:

getId in interface ComputeBackend

getName

public String getName()

Description copied from interface: Backend

Returns the display name for UI presentation.

Specified by:

getName in interface AlgorithmProvider

Specified by:

getName in interface Backend

Specified by:

getName in interface LinearAlgebraProvider<E extends FieldElement<E>>

getPriority

public int getPriority()

Description copied from interface: LinearAlgebraProvider

Priority of this provider (higher means more preferred). Used for automatic backend selection.

Specified by:

getPriority in interface AlgorithmProvider

Specified by:

getPriority in interface Backend

Specified by:

getPriority in interface LinearAlgebraProvider<E extends FieldElement<E>>

shutdown

public void shutdown()

Description copied from interface: AlgorithmProvider

Called when the provider is no longer needed (e.g., application shutdown). Use this to release local resources or close native segments.

Specified by:

shutdown in interface AlgorithmProvider

Specified by:

shutdown in interface Backend

multiply

public Matrix<E> multiply(Matrix<E> a,
 Matrix<E> b)

Specified by:

multiply in interface LinearAlgebraProvider<E extends FieldElement<E>>

transpose

public Matrix<E> transpose(Matrix<E> a)

Specified by:

transpose in interface LinearAlgebraProvider<E extends FieldElement<E>>

trace

public E trace(Matrix<E> a)

Specified by:

trace in interface LinearAlgebraProvider<E extends FieldElement<E>>

conjugateTranspose

public Matrix<E> conjugateTranspose(Matrix<E> a)

Specified by:

conjugateTranspose in interface LinearAlgebraProvider<E extends FieldElement<E>>

normalize

public Vector<E> normalize(Vector<E> v)

Description copied from interface: LinearAlgebraProvider

Returns the normalized vector (unit vector).

Specified by:

normalize in interface LinearAlgebraProvider<E extends FieldElement<E>>

cross

public Vector<E> cross(Vector<E> a,
 Vector<E> b)

Description copied from interface: LinearAlgebraProvider

Returns the cross product of two 3D vectors.

Specified by:

cross in interface LinearAlgebraProvider<E extends FieldElement<E>>

angle

public E angle(Vector<E> a,
 Vector<E> b)

Description copied from interface: LinearAlgebraProvider

Returns the angle between two vectors in radians.

Specified by:

angle in interface LinearAlgebraProvider<E extends FieldElement<E>>

subtract

public Matrix<E> subtract(Matrix<E> a,
 Matrix<E> b)

Specified by:

subtract in interface LinearAlgebraProvider<E extends FieldElement<E>>

scale

public Matrix<E> scale(E scalar,
 Matrix<E> a)

Specified by:

scale in interface LinearAlgebraProvider<E extends FieldElement<E>>

solve

public Vector<E> solve(Matrix<E> a,
 Vector<E> b)

Specified by:

solve in interface LinearAlgebraProvider<E extends FieldElement<E>>

inverse

public Matrix<E> inverse(Matrix<E> a)

Specified by:

inverse in interface LinearAlgebraProvider<E extends FieldElement<E>>

determinant

public E determinant(Matrix<E> a)

Specified by:

determinant in interface LinearAlgebraProvider<E extends FieldElement<E>>

lu

public LUResult<E> lu(Matrix<E> a)

Description copied from interface: LinearAlgebraProvider

Computes the LU decomposition of the specified matrix.

Specified by:

lu in interface LinearAlgebraProvider<E extends FieldElement<E>>

qr

public QRResult<E> qr(Matrix<E> a)

Description copied from interface: LinearAlgebraProvider

Computes the QR decomposition of the specified matrix.

Specified by:

qr in interface LinearAlgebraProvider<E extends FieldElement<E>>

svd

public SVDResult<E> svd(Matrix<E> a)

Description copied from interface: LinearAlgebraProvider

Computes the Singular Value Decomposition (SVD) of the specified matrix.

Specified by:

svd in interface LinearAlgebraProvider<E extends FieldElement<E>>

cholesky

public CholeskyResult<E> cholesky(Matrix<E> a)

Description copied from interface: LinearAlgebraProvider

Computes the Cholesky decomposition of the specified matrix.

Specified by:

cholesky in interface LinearAlgebraProvider<E extends FieldElement<E>>

eigen

public EigenResult<E> eigen(Matrix<E> a)

Description copied from interface: LinearAlgebraProvider

Computes the eigenvalue decomposition of the specified matrix.

Specified by:

eigen in interface LinearAlgebraProvider<E extends FieldElement<E>>

solveTriangular

public Vector<E> solveTriangular(Matrix<E> a,
 Vector<E> b,
 boolean upper,
 boolean transpose,
 boolean conjugate,
 boolean unit)

Description copied from interface: LinearAlgebraProvider

Solves the triangular system Ax = b.

Specified by:

solveTriangular in interface LinearAlgebraProvider<E extends FieldElement<E>>

Parameters:

a - the triangular matrix

b - the right-hand side vector

upper - true if A is upper triangular, false if lower

transpose - true if solving A^T x = b

conjugate - true if solving A^H x = b (only if transpose is true and complex)

unit - true if A is unit triangular (diagonal is all ones)

add

public Vector<E> add(Vector<E> a,
 Vector<E> b)

Specified by:

add in interface LinearAlgebraProvider<E extends FieldElement<E>>

subtract

public Vector<E> subtract(Vector<E> a,
 Vector<E> b)

Specified by:

subtract in interface LinearAlgebraProvider<E extends FieldElement<E>>

multiply

public Vector<E> multiply(Vector<E> v,
 E s)

Specified by:

multiply in interface LinearAlgebraProvider<E extends FieldElement<E>>

multiply

public Vector<E> multiply(Matrix<E> a,
 Vector<E> b)

Specified by:

multiply in interface LinearAlgebraProvider<E extends FieldElement<E>>

dot

public E dot(Vector<E> a,
 Vector<E> b)

Specified by:

dot in interface LinearAlgebraProvider<E extends FieldElement<E>>

norm

public E norm(Vector<E> a)

Specified by:

norm in interface LinearAlgebraProvider<E extends FieldElement<E>>

getAcceleratorType

public HardwareAccelerator getAcceleratorType()

Description copied from interface: ComputeBackend

Returns the type of hardware accelerator used by this backend.

Specified by:

getAcceleratorType in interface ComputeBackend

Returns:

the accelerator type

getType

public String getType()

Description copied from interface: Backend

Returns the backend type category (e.g., "plotting", "molecular", "tensor", "math").

Specified by:

getType in interface Backend

Specified by:

getType in interface ComputeBackend

close

public void close()

Specified by:

close in interface AutoCloseable

Specified by:

close in interface LinearAlgebraProvider<E extends FieldElement<E>>

createContext

public ExecutionContext createContext()

Description copied from interface: ComputeBackend

Creates an execution context for running operations.

The context should be used within a try-with-resources block to ensure proper cleanup of resources.

Specified by:

createContext in interface ComputeBackend

Returns:

new execution context

getDevices

public GPUBackend.DeviceInfo[] getDevices()

Description copied from interface: GPUBackend

Returns information about available GPU devices.

Specified by:

getDevices in interface GPUBackend

selectDevice

public void selectDevice(int deviceId)

Description copied from interface: GPUBackend

Selects a GPU device for computation.

Specified by:

selectDevice in interface GPUBackend

Parameters:

deviceId - Device index (0-based)

allocateGPUMemory

public long allocateGPUMemory(long sizeBytes)

Description copied from interface: GPUBackend

Allocates GPU memory and returns a handle.

Specified by:

allocateGPUMemory in interface GPUBackend

Parameters:

sizeBytes - Size in bytes

Returns:

GPU memory handle (opaque pointer)

copyToGPU

public void copyToGPU(long gpuHandle,
 DoubleBuffer hostBuffer,
 long sizeBytes)

Description copied from interface: GPUBackend

Copies data from host to GPU.

Specified by:

copyToGPU in interface GPUBackend

Parameters:

gpuHandle - GPU memory handle

hostBuffer - Host buffer

sizeBytes - Number of bytes to copy

copyFromGPU

public void copyFromGPU(long gpuHandle,
 DoubleBuffer hostBuffer,
 long sizeBytes)

Description copied from interface: GPUBackend

Copies data from GPU to host.

Specified by:

copyFromGPU in interface GPUBackend

Parameters:

gpuHandle - GPU memory handle

hostBuffer - Host buffer

sizeBytes - Number of bytes to copy

freeGPUMemory

public void freeGPUMemory(long gpuHandle)

Description copied from interface: GPUBackend

Frees GPU memory.

Specified by:

freeGPUMemory in interface GPUBackend

Parameters:

gpuHandle - GPU memory handle

synchronize

public void synchronize()

Description copied from interface: GPUBackend

Synchronizes GPU execution (waits for all kernels to complete).

Specified by:

synchronize in interface GPUBackend

matrixMultiply

public void matrixMultiply(DoubleBuffer A,
 DoubleBuffer B,
 DoubleBuffer C,
 int m,
 int n,
 int k)

Description copied from interface: GPUBackend

Performs matrix multiplication C = A * B on the GPU. Dimensions: A(m x k), B(k x n), C(m x n).

Specified by:

matrixMultiply in interface GPUBackend