-
Notifications
You must be signed in to change notification settings - Fork 0
/
.devnote
61 lines (49 loc) · 2.05 KB
/
.devnote
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
----- This file contains personal notes to keep track of ideas, TODO or previous chunks of code which has been removed but could be useful in different context -----
----- Not supposed to be read by anyone other than me -----
cmake -DCMAKE_CXX_COMPILER=clang++-19 -DCMAKE_BUILD_TYPE=Debug -DCMAKE_CXX_FLAGS="-Wfatal-errors" ..
// Extract Elem from a row of YoungTableauSeq at index J
template <std::size_t J, std::size_t Id, class Row>
struct ExtractElem;
template <std::size_t J, std::size_t Id, std::size_t HeadElem, std::size_t... TailElem>
struct ExtractElem<J, Id, std::index_sequence<HeadElem, TailElem...>>
{
using type = std::conditional_t<
J == Id,
std::index_sequence<HeadElem>,
typename ExtractElem<J, Id + 1, std::index_sequence<TailElem...>>::type>;
};
template <std::size_t J, class Row>
using extract_elem_t = ExtractElem<J, 0, Row>::type;
// Type of index used by projectors
template <std::size_t I>
struct ProjX
{
};
template <class Ids>
struct NaturalIndex;
template <std::size_t... Id>
struct NaturalIndex<std::index_sequence<Id...>>
{
template <std::size_t RankId>
struct type : sil::tensor::TensorNaturalIndex<ProjX<Id>...>
{
};
};
template <class NaturalIds, class RankIds>
struct ProjectorIndex;
template <std::size_t... Id, std::size_t... RankId>
struct ProjectorIndex<NaturalIndex<std::index_sequence<Id...>>, std::index_sequence<RankId...>>
{
using type = sil::tensor::FullTensorIndex<
typename NaturalIndex<std::index_sequence<Id...>>::template type<RankId>...>;
};
template <std::size_t Dimension, std::size_t Rank>
using projector_index_t = ProjectorIndex<
NaturalIndex<std::make_index_sequence<Dimension>>,
std::make_index_sequence<2 * Rank>>::type;
----- TODO -----
- apply_lambda probably useless (can be written with ddc::for_each)
- IsTensorIndex::value has to be set to false
- rename mem_id
- sil::tensor::detail::Access implementation is not properly implemented
- specialize tensor_prod using concepts (ie. redirect to natural_tensor_prod if NaturalTensorIndex)