Skip to content
GitLab
Projects
Groups
Snippets
Help
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
Open sidebar
xQML
xQML
Commits
919aa2e2
Commit
919aa2e2
authored
Feb 11, 2019
by
Syl
Browse files
Merge branch 'master' of gitlab.in2p3.fr:xQML/xQML
i /bin/bash: q: command not found
parents
a213f533
da2d72f3
Changes
2
Hide whitespace changes
Inline
Side-by-side
Showing
2 changed files
with
7 additions
and
5 deletions
+7
-5
src/libcov.c
src/libcov.c
+4
-2
xqml/estimators.py
xqml/estimators.py
+3
-3
No files found.
src/libcov.c
View file @
919aa2e2
...
...
@@ -6,14 +6,16 @@ void build_dSdC( int nside, int nstokes, int npix, int nbin, long *ellbins, long
const
int
lmax1
=
lmax
+
1
;
const
int
ns
=
3
,
nspecall
=
4
;
const
int
nspec
=
nstokes2nspec
(
nstokes
);
const
int
npixtot
=
npix
*
nstokes
;
const
int
64_t
npixtot
=
npix
*
nstokes
;
// fprintf( stdout, "lmax=%d\n", lmax);
// fprintf( stdout, "npix=%d\n", npix);
// fprintf( stdout, "nstokes=%d\n", nstokes);
// fprintf( stdout, "nspec=%d\n", nspec);
memset
(
dSdC
,
0
.,
(
nspec
*
nbin
*
nstokes
*
npix
*
nstokes
*
npix
)
*
sizeof
(
double
));
int64_t
ntot
=
nspec
*
nbin
*
npixtot
*
npixtot
;
// fprintf( stdout, "memset (%d,%d,%d,%d = %lld)...\n", nspec,nbin,nstokes,npix,ntot);
memset
(
dSdC
,
0
.,
ntot
*
sizeof
(
double
));
#pragma omp parallel default(none) shared(stdout,nbin,nside,npix,nstokes,dSdC,ipix,bl,ellbins)
{
...
...
xqml/estimators.py
View file @
919aa2e2
...
...
@@ -146,7 +146,7 @@ def El(invCAA, invCBB, Pl):
def
CrossWindowFunction
(
El
,
Pl
):
"""
Compute mode-mixi
g
n matrix (Tegmark's window matrix)
Compute mode-mixin
g
matrix (Tegmark's window matrix)
Wll = Trace[invCAA.Pl.invCBB.Pl] = Trace[El.Pl]
Parameters
...
...
@@ -181,7 +181,7 @@ def CrossWindowFunction(El, Pl):
def
CrossWindowFunctionLong
(
invCAA
,
invCBB
,
Pl
):
"""
Compute mode-mixi
g
n matrix (Tegmark's window matrix)
Compute mode-mixin
g
matrix (Tegmark's window matrix)
Wll = Trace[invCAA.Pl.invCBB.Pl] = Trace[El.Pl]
Parameters
...
...
@@ -437,4 +437,4 @@ if __name__ == "__main__":
import
doctest
if
np
.
__version__
>=
"1.14.0"
:
np
.
set_printoptions
(
legacy
=
"1.13"
)
doctest
.
testmod
()
\ No newline at end of file
doctest
.
testmod
()
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment