Hacktoberfest 2026: le issue che i maintainer hanno segnato per ottobre, aperte e adatte ai principianti. Sfoglia le issue Hacktoberfest

Improper ell cuts when loading spectrum from file

Aperta
  • Qualcosa non va?
#119 0 commenti 0 reazioni 0 assegnatari Vedi su GitHub

Nessuno ha ancora preso questa issue.

Valutazione

Difficoltà
3/5
Tempo stimato
1-2 giorni
Idoneità per principianti
45/100
Tipo di issue
Bug
Chiarezza
Abbastanza chiara
Stato di attività
Ferma
Stack tecnologico
python
Ambito
data

Direzione di ricerca

Inizia da pspy_utils.ps_lensed_theory_to_dict e confronta le sue convenzioni per lmax e per l'inizio degli array con quelle delle altre funzioni pspy che accettano lmax. Traccia gli input che iniziano a l=0 o l=2, incluso start_at_zero, e verifica che i multipoli e gli spettri restituiti coprano coerentemente l'intervallo richiesto, senza anteposizioni o spostamenti involontari.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

The function pspy_utils.ps_lensed_theory_to_dict loads a plaintext file and returns a dictionary containing the power spectra, but the way in which the user-requested lmax is handled doesn't work properly. Cf. the code:

def ps_lensed_theory_to_dict(filename, output_type, lmax=None, start_at_zero=False):
    fields = ["TT", "TE", "TB", "ET", "BT", "EE", "EB", "BE", "BB"]
    ps = {}
    l, ps["TT"], ps["EE"], ps["BB"], ps["TE"] = np.loadtxt(filename, unpack=True)
    ps["ET"] = ps["TE"].copy()
    ps["TB"], ps["BT"], ps["EB"], ps["BE"] = np.zeros((4, len(l)))

    if lmax is not None:
        l = l[:lmax]
    scale = l * (l + 1) / (2 * np.pi)
    for f in fields:
        if lmax is not None:
            ps[f] = ps[f][:lmax]
        if output_type == "Cl":
            ps[f] /= scale
        if start_at_zero:
            ps[f] = np.append(np.array([0, 0]), ps[f])
    if start_at_zero:
        l = np.append(np.array([0, 1]), l)
    return l, ps
  1. The lmax cuts are parsed as l = l[:lmax] and ps[f] = ps[f][:lmax], both of which work fine if the array is zero-indexed. However, no check is made to see if the arrays actually start at l = 0.
  2. If the user desires the arrays start at l = 0 instead of l = 2, then two zeros are prepended before the array, but no check is done whether this is needed.

Most other functions in pspy that take lmax as an argument will assume the user wants an array from l = 2 to lmax-1 inclusive, but this function will most likely create arrays from l=2 to lmax+1 inclusive because of the aforementioned design, and more often than not the resulting spectra will create issues when passed on to other functions that expect [2,lmax) arrays.

Lingua principale
Python
Stelle
10
Fork
8
Merge medio
21m
PR unite (30g)
2

Preparare l'ambiente

Questo progetto non fornisce container di sviluppo, Dockerfile né guida per i contributori, quindi l'ambiente è a tuo carico: parti dal suo README e consulta la nostra guida al primo contributo per i passaggi generali.

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di simonsobs/pspy

Tutte le issue di simonsobs/pspy

Issue simili

Altre issue su Python

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.