-
-
Notifications
You must be signed in to change notification settings - Fork 1
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
fefa834
commit 9342f29
Showing
7 changed files
with
113 additions
and
3,093 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,62 @@ | ||
| module Structs | ||
| export RTPSA, CTPSA, Desc | ||
| const NAMSZ::Int = 16 | ||
|
|
||
| mutable struct RTPSA{T} | ||
| d::Ptr{T} # Ptr to tpsa descriptor | ||
| uid::Cint # Special user field for external use (and padding) | ||
| mo::Cuchar # max ord (allocated) | ||
| lo::Cuchar # lowest used ord | ||
| hi::Cuchar # highest used ord | ||
| nz::Culonglong # zero/nonzero homogenous polynomials. Int64 if 64 bit else 32 bit | ||
| nam::NTuple{NAMSZ,Cuchar} # tpsa name max string length 16 NAMSZ | ||
| coef::Ptr{Cdouble} # warning: must be identical to ctpsa up to coef excluded | ||
| end | ||
|
|
||
| mutable struct CTPSA{T} | ||
| d::Ptr{T} # Ptr to ctpsa descriptor | ||
| uid::Cint # Special user field for external use (and padding) | ||
| mo::Cuchar # max ord (allocated) | ||
| lo::Cuchar # lowest used ord | ||
| hi::Cuchar # highest used ord | ||
| nz::Culonglong # zero/nonzero homogenous polynomials. Int64 if 64 bit else 32 bit | ||
| nam::NTuple{NAMSZ,Cuchar} # tpsa name | ||
| coef::Ptr{ComplexF64} # warning: must be identical to ctpsa up to coef excluded | ||
| end | ||
|
|
||
| struct Desc{T,C} | ||
| id::Cint # index in list of registered descriptors | ||
| nn::Cint # nn = nv+np <= 100000 | ||
| nv::Cint # nv = number of variables | ||
| np::Cint # np = number of parameters | ||
| mo::Cuchar # max orders of vars | ||
| po::Cuchar # max orders of params | ||
| to::Cuchar # global order of truncation. Note ord_t in mad_tpsa is typedef for unsigned char (Cuchar) | ||
| no::Ptr{Cuchar} # orders of each vars and params, no[nn]. In C this is const | ||
|
|
||
| uno::Cint # user provided no | ||
| nth::Cint # max #threads or 1 | ||
| nc::Cuint # number of coefs (max length of TPSA) | ||
|
|
||
| monos::Ptr{Cuchar} # 'matrix' storing the monomials (sorted by var) | ||
| ords::Ptr{Cuchar} # Order of each mono of To | ||
| To::Ptr{Ptr{Cuchar}} # Table by orders -- pointers to monos, sorted by order | ||
| Tv::Ptr{Ptr{Cuchar}} # Table by vars -- pointers to monos, sorted by vars | ||
| ocs::Ptr{Ptr{Cuchar}} # ocs[t,i] -> o; in mul, compute o on thread t; 3 <= o <= mo; terminated with 0 | ||
|
|
||
| ord2idx::Ptr{Cint} # order to polynomial start index in To (i.e. in TPSA coef[]) | ||
| tv2to::Ptr{Cint} # lookup tv->to | ||
| to2tv::Ptr{Cint} # lookup to->tv | ||
| H::Ptr{Cint} # indexing matrix in Tv | ||
| L::Ptr{Ptr{Cint}} # multiplication indexes L[oa,ob]->L_ord; L_ord[ia,ib]->ic | ||
| L_idx::Ptr{Ptr{Ptr{Cint}}} # L_idx[oa,ob]->[start] [split] [end] idxs in L | ||
|
|
||
| size::Culonglong # bytes used by desc. Unsigned Long Int, ikn 32 bit system is int32 but 64 bit int64. Using Culonglong assuming 64 bit | ||
|
|
||
| t::Ptr{Ptr{T}} # tmp for tpsa | ||
| ct::Ptr{Ptr{C}} # tmp for ctpsa | ||
| ti::Ptr{Cint} # idx of tmp ised | ||
| cti::Ptr{Cint} # idx of tmp used | ||
| end | ||
|
|
||
| end |
Oops, something went wrong.