Skip to main content

tensor4all_core/
lib.rs

1#![warn(missing_docs)]
2//! Core tensor operations and types for tensor4all-rs.
3//!
4//! This crate provides the foundational types and operations for tensor networks:
5//!
6//! - **Index types**: [`DynIndex`], [`Index`], [`DynId`] for tensor indices
7//! - **Tag sets**: [`TagSet`], [`TagSetLike`] for metadata tagging
8//! - **Tensors**: [`IdxTensor`] for dynamic-rank dense tensors
9//! - **Operations**: Contraction, SVD, QR decomposition, factorization
10//!
11//! # Example
12//!
13//! ```
14//! use tensor4all_core::{Index, DynIndex, IdxTensor};
15//!
16//! // Create indices with dynamic identity
17//! let i = Index::new_dyn(2);
18//! let j = Index::new_dyn(3);
19//!
20//! // Create a tensor
21//! let data = vec![1.0, 2.0, 3.0, 4.0, 5.0, 6.0];
22//! let t = IdxTensor::from_dense(vec![i.clone(), j.clone()], data).unwrap();
23//! ```
24
25#[cfg(doctest)]
26#[doc = include_str!("../README.md")]
27pub struct ReadmeDoctests;
28
29pub mod prelude;
30
31pub mod col_major_array;
32pub use col_major_array::{ColMajorArray, ColMajorArrayMut, ColMajorArrayRef};
33
34// Common (tags, utilities, scalar)
35/// Dynamic scalar compatibility wrapper built on rank-0 `IdxTensor`.
36pub mod any_scalar;
37pub mod global_default;
38/// Bit-packed integer keys for multi-index maps.
39pub mod index_key;
40pub mod index_like;
41pub mod scalar;
42/// Stack-allocated fixed-capacity string types for ITensors.jl compatibility.
43pub mod smallstring;
44/// Tag set types for tensor metadata.
45pub mod tagset;
46pub mod truncation;
47
48// TCI substrate absorbed from tensor4all-tcicore (dissolved, #639): matrix
49// cross interpolation (LU/LUCI/ACA), common scalar trait, floating-zone
50// pivot search, caching, and index sets.
51pub mod cached_function;
52pub mod error;
53pub mod floating_zone;
54pub mod indexset;
55mod matrix_luci;
56pub mod matrixaca;
57pub mod matrixlu;
58pub mod matrixluci;
59pub mod traits;
60
61pub use self::matrixluci::MatrixLuciScalar;
62pub use cached_function::cache_key::CacheKey;
63pub use cached_function::error::CacheKeyError;
64pub use cached_function::index_int::IndexInt;
65pub use cached_function::CachedFunction;
66pub use error::{MatrixCIError, Result};
67pub use floating_zone::floating_zone_walk;
68pub use indexset::{IndexSet, LocalIndex, MultiIndex};
69pub use matrix_luci::{
70    matrix_luci_factors_from_blocks, matrix_luci_factors_from_matrix,
71    matrix_luci_factors_from_matrix_owned, MatrixLUCI, MatrixLuciFactors,
72};
73pub use matrixaca::MatrixACA;
74pub use matrixlu::{rrlu, rrlu_mut, RrLU, RrLUOptions};
75pub use scalar::Scalar;
76pub use traits::AbstractMatrixCI;
77
78pub use scalar::Scalar as CommonScalar;
79
80// Default concrete type implementations (index, tensor, linalg, etc.)
81pub mod defaults;
82
83// Backwards compatibility: re-export defaults submodules as top-level modules
84// This allows `tensor4all_core::index::...` to work
85pub use defaults::index;
86
87pub use defaults::{DefaultIndex, DefaultTagSet, DynId, DynIndex, Index, TagSet};
88pub use index_like::{sort_indices_deterministic, ConjState, IndexLike};
89
90/// Index operations (replacement, set operations, contraction preparation).
91pub mod index_ops;
92pub use index_ops::{
93    check_unique_indices, common_ind_positions, common_inds, has_common_inds, has_inds, hasind,
94    noncommon_inds, replace_indices, replace_indices_mut, union_inds, unique_inds,
95    ReplaceIndsError,
96};
97pub use smallstring::{SmallChar, SmallString, SmallStringError};
98pub use tagset::{Tag, TagSetError, TagSetLike};
99
100// Tensor (storage, tensor types)
101pub mod tensor_index;
102pub mod tensor_like;
103
104pub use tensor_index::TensorIndex;
105
106// Krylov subspace methods (GMRES, etc.)
107pub mod krylov;
108
109// Block tensor for block matrix GMRES
110pub mod block_tensor;
111
112// Backwards compatibility: re-export defaults::idx_tensor as tensor
113pub use defaults::idx_tensor as tensor;
114
115pub use any_scalar::{AnyScalar, AnyScalarError};
116pub use defaults::idx_tensor::{
117    compute_permutation_from_indices, diag_idx_tensor, unfold_split, IdxTensor, IdxTensorError,
118    StructuredSelectorError, TensorStorageError,
119};
120#[cfg(feature = "tenferro-cuda")]
121pub use defaults::IdxTensorCudaError;
122pub use tensor4all_tensorbackend::TensorElement;
123pub use tensor4all_tensorbackend::{
124    print_and_reset_native_einsum_profile, reset_native_einsum_profile,
125};
126#[cfg(feature = "tenferro-cuda")]
127pub use tensor4all_tensorbackend::{CudaExecutionContext, CudaExecutionContextError, CUDA_ORDINAL};
128pub use tensor_like::{
129    Canonical, DirectSumResult, FactorizeAlg, FactorizeError, FactorizeOptions, FactorizeResult,
130    LinearizationOrder, TensorConstructionLike, TensorContractionLike, TensorFactorizationLike,
131    TensorLike, TensorVectorSpace, TensorVectorSpaceError,
132};
133
134pub use defaults::contract::{
135    contract, contract_owned, contract_owned_with_options, contract_pair,
136    contract_pair_with_operand_options, contract_pair_with_options, contract_with_options,
137    outer_product, print_and_reset_contract_profile, reset_contract_profile, tensordot,
138    ContractionOptions, PairwiseContractionOptions,
139};
140pub use defaults::idx_tensor::{
141    print_and_reset_pairwise_contract_profile, reset_pairwise_contract_profile,
142};
143
144// Re-export linear algebra modules from defaults for backwards compatibility
145// This allows `tensor4all_core::svd::...`, `tensor4all_core::qr::...`, etc.
146pub mod direct_sum {
147    //! Re-export of direct sum operations.
148    pub use crate::defaults::direct_sum::*;
149}
150pub mod factorize {
151    //! Re-export of factorization operations.
152    pub use crate::defaults::factorize::*;
153}
154pub mod qr {
155    //! Re-export of QR decomposition operations.
156    pub use crate::defaults::qr::*;
157}
158pub mod svd {
159    //! Re-export of SVD decomposition operations.
160    pub use crate::defaults::svd::{
161        default_svd_truncation_policy, set_default_svd_truncation_policy, svd, svd_with, SvdError,
162        SvdOptions,
163    };
164}
165
166// Re-export linear algebra items for top-level access
167pub use defaults::direct_sum::direct_sum;
168pub use defaults::factorize::{factorize, factorize_full_rank};
169pub use defaults::qr::{default_qr_rtol, qr, qr_with, set_default_qr_rtol, QrError, QrOptions};
170pub use defaults::svd::{
171    default_svd_truncation_policy, set_default_svd_truncation_policy, svd, svd_with, SvdError,
172    SvdOptions,
173};
174
175// Global default and truncation utilities
176pub use global_default::{GlobalDefault, InvalidRtolError};
177pub use truncation::{
178    validate_svd_truncation_options, DecompositionAlg, InvalidThresholdError, SingularValueMeasure,
179    SvdTruncationOptionsError, SvdTruncationPolicy, ThresholdScale, TruncationRule,
180};