https://kaihsin.github.io/Cytnx/docs/html/index.html
https://hub.docker.com/r/kaihsinwu/cytnx_mkl
* To run:
$docker pull kaihsinwu/cytnx_mkl
$docker run -ti kaihsinwu/cytnx_mkl
* Note:
Once docker image is run, the user code can be compile (for example) with:
$g++-6 -std=c++11 -O3 <your.code.cpp> /opt/cytnx/libcytnx.so
Or, directly using python API, by import cytnx in python script:
from cytnx import ** Boost v1.53+ [check_deleted, atomicadd, intrusive_ptr]
* C++11
* lapack
* blas
* gcc v4.8.5+ (recommand v6+) (required -std=c++11)
[CUDA support]
* CUDA v10+
* cuDNN
[OpenMp support]
* openmp
[Python]
* pybind11 2.2.4
* numpy >= 1.15
[MKL]
* icpc (if using with CUDA, v18- should be used)
sudo apt-get install libboost-all-dev libblas-dev liblapack-dev make curl g++-6 libomp-dev
1.) Set the flags to config the install inside make.inc
There are 4 important flags: ICPC_Enable, OMP_Enable, GPU_Enable, MKL_Enable
* a. The default compiler is g++-6. Change "GCC" for the compiler on your system.
To use intel icpc compiler instead of default compiler, set "ICPC_Enable"=1.
[Note] You can only choose either icpc or gcc.
* In the case where ICPC_Enable=1, GCC will be ignore.
* In the case where ICPC_Enable=1, set "ICPC" to the path of your icpc binary.
* In the case where ICPC_Enable=0, "ICPC" will be ignored.
* b. To enable the GPU support, set "GPU_Enable" =1, otherwise =0.
[Note]
* if GPU_Enable=1, the "CUDA_PATH" should be set to the cuda directory on your system.
* if GPU_Enable=0, the "CUDA_PATH" will be ignored.
* c. To enable the acceleration using OpenMP, set "OMP_Enable" =1, otherwise =0.
* d. The default linalg library are using LAPACK and BLAS. To use intel MKL instead, set "MKL_enable" =1.
2.) compile by running:
$make -Bj
3.) [Option] compile the python API
$make pyobj
* Storage [binded]
* Tensor [binded]
* Bond [binded]
* Accessor [c++ only]
* Symmetry [binded]
* UniTensor [binded]
* Network [binded]
Benefit from both side.
One can do simple prototype on python side
and easy transfer to C++ with small effort!
// c++ version:
#include "cytnx.hpp"
cytnx::Tensor A({3,4,5},cytnx::Type.Double,cytnx::Device.cpu) # python version:
import cytnx
A = cytnx.Tensor((3,4,5),dtype=cytnx.Type.Double,device=cytnx.Device.cpu) The avaliable types are :
| cytnx type | c++ type | Type object
|------------------|----------------------|--------------------
| cytnx_double | double | Type.Double
| cytnx_float | float | Type.Float
| cytnx_uint64 | uint64_t | Type.Uint64
| cytnx_uint32 | uint32_t | Type.Uint32
| cytnx_uint16 | uint16_t | Type.Uint16
| cytnx_int64 | int64_t | Type.Int64
| cytnx_int32 | int32_t | Type.Int32
| cytnx_int16 | int16_t | Type.Int16
| cytnx_complex128 | std::complex<double> | Type.ComplexDouble
| cytnx_complex64 | std::complex<float> | Type.ComplexFloat
| cytnx_bool | bool | Type.Bool
* Memory container with GPU/CPU support.
maintain type conversions (type casting btwn Storages)
and moving btwn devices.
* Generic type object, the behavior is very similar to python.
Storage A(400,Type.Double);
for(int i=0;i<400;i++)
A.at<double>(i) = i;
Storage B = A; // A and B share same memory, this is similar as python
Storage C = A.to(Device.cuda+0); * A tensor, API very similar to numpy and pytorch.
* simple moving btwn CPU and GPU:
Tensor A({3,4},Type.Double,Device.cpu); // create tensor on CPU (default)
Tensor B({3,4},Type.Double,Device.cuda+0); // create tensor on GPU with gpu-id=0
Tensor C = B; // C and B share same memory.
// move A to gpu
Tensor D = A.to(Device.cuda+0);
// inplace move A to gpu
A.to_(Device.cuda+0); * Type conversion in between avaliable:
Tensor A({3,4},Type.Double);
Tensor B = A.astype(Type.Uint64); // cast double to uint64_t * vitual swap and permute. All the permute and swap will not change the underlying memory
* Use Contiguous() when needed to actual moving the memory layout.
Tensor A({3,4,5,2},Type.Double);
A.permute_({0,3,1,2}); // this will not change the memory, only the shape info is changed.
cout << A.is_contiguous() << endl; // this will be false!
A.contiguous_(); // call Configuous() to actually move the memory.
cout << A.is_contiguous() << endl; // this will be true! * access single element using .at
Tensor A({3,4,5},Type.Double);
double val = A.at<double>({0,2,2}); * access elements with python slices similarity:
typedef Accessor ac;
Tensor A({3,4,5},Type.Double);
Tensor out = A.get({ac(0),ac::all(),ac::range(1,4)});
// equivalent to python: out = A[0,:,1:4]
See example/ folder or documentation for how to use API
See test.cpp for using C++ .
See test.py for using python
func | inplace | CPU | GPU | callby tn | Tn | UniTn
--------------|-----------|-----|------|-------------|----|-------
Add | x | Y | Y | Y | Y | Y
Sub | x | Y | Y | Y | Y | Y
Mul | x | Y | Y | Y | Y | Y
Div | x | Y | Y | Y | Y | Y
Cpr | x | Y | Y | Y | Y | x
+,+=[tn] | x | Y | Y | Y (Add_) | Y | Y
-,-=[tn] | x | Y | Y | Y (Sub_) | Y | Y
*,*=[tn] | x | Y | Y | Y (Mul_) | Y | Y
/,/=[tn] | x | Y | Y | Y (Div_) | Y | Y
==[tn] | x | Y | Y | Y (Cpr_) | Y | x
Svd | x | Y | Y | Y | Y | N
*Svd_truncate| x | Y | Y | N | Y | N
Inv | Inv_ | Y | Y | Y | Y | N
Conj | Conj_ | Y | Y | Y | Y | N
Exp | Exp_ | Y | Y | Y | Y | N
Eigh | x | Y | Y | Y | Y | N
Matmul | x | Y | Y | N | Y | N
Diag | x | Y | Y | N | Y | N
*Tensordot | x | Y | Y | N | Y | N
Otimes | x | Y | Y | N | Y | N
* this is a high level linalg
Tensor: zeros(), ones(), arange()
Kai-Hsin Wu kaihsinwu@gmail.com