Comments (8)
@DanielLenz. Sorry, I just saw this. I assume you are using nside=2048
, is this right? Are you including foregrounds or just the raw data? Is this for just temperature of temp+polarization? I have to check but if it's temp+polarization auto and cross-correlations for nside=2048
, ~30 minutes sounds about right with ~8 threads?
from namaster.
These times sound about right to me (especially if the 30min include computing the mode-coupling matrix). I just checked that a simple healpy.alm2map
on 2048 maps takes about 2min on my own laptop (which I'd also call high-end).
About lmax
, this is something that is half-implemented (at the C level, but not visible in python), and which I wanted to include as an option once it's better tested.
from namaster.
Thanks for your replies!
Are you including foregrounds or just the raw data? Is this for just temperature of temp+polarization?
I just tested this with the Planck CMB map, nside=2048, temperature only, no foregrounds. I might just be spoiled by the PolSpice implementation, which runs in ~< 1min. It is, however, only approximate and has several other issues that make me really want to switch to nmt, and I'll gladly accept the slightly longer runtime.
Once we have the optional lmax
keyword and can set this to ~2048, the computation of the mode-coupling matrix should go down by a factor of ~(3*2048/2048)^3 = 27 to ~1 minute, which is totally acceptable.
Please feel free to close this if you want, but it would be great if could update this thread once an lmax
branch is available.
from namaster.
@DanielLenz now that I think about it, this may be already possible. If you pass an NmtBin structure that only contains bandpowers up to ell=2048, the computation of the MCM should only go up to that. Let me know if this does or doesn't work!
Note however that you probably want to throw away the last bin or two, since the mode coupling will not be correctly represented for those due to the missing multipoles.
from namaster.
That works perfectly, thank you! It might be worth to add this to the docs, even though it's hinted at in the NmtBin documentation.
from namaster.
OK, great. Will update the docs. Can I check that you've seen that this scales as expected?
from namaster.
Just did a quick analysis, it appears to scale better than lmax^3 for lmax < 2000. I could run larger values of lmax later as well.
from namaster.
Awesome, thanks a lot @DanielLenz!
from namaster.
Related Issues (20)
- nan's with map read from fits file HOT 2
- random seed in synfast_flat
- I am not seeing a real improvement in performance, but it might be my own compiler being a bit outdated. Still worth merging, since libsharp2 is the one being maintained, but will open an issue to investigate further.
- configure fails to find cfitsio HOT 7
- Missing version number in pymaster
- Help installing NaMaster on M2 macOS HOT 17
- cannot find _nmtbin module HOT 4
- Error when setting lmax < 3 * nside - 1 in compute_full_master HOT 2
- Missing different lmax for mask and coupling coefficients for covariances HOT 1
- f_ell does not seem to do anything in NmtBin HOT 1
- Conda installation fails if conda installs gsl=2.7.1 instead of gsl=2.7
- Inconsistency in results between MacOS and Ubuntu HOT 13
- Make exception message more informative
- Weird power in compute_coupled_cell HOT 3
- pymaster installation problem HOT 3
- Pseudo-cl's in the covariance computation - pixel window function HOT 3
- Floating point errors with full sky CAR pixellization HOT 1
- NaMaster installation CFITSIO issue
- Apodisation window functions HOT 3
- Noise bias for flat sky HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from namaster.