Commit f3300038 by Martin Reinecke

### enhance demo

parent 0ee1363a
 ... @@ -30,11 +30,11 @@ job = pysharp.sharpjob_d() ... @@ -30,11 +30,11 @@ job = pysharp.sharpjob_d() # m<=l<=lmax. # m<=l<=lmax. # Symmetry: a_l,-m = (-1)**m*conj(a_l,m). # Symmetry: a_l,-m = (-1)**m*conj(a_l,m). # The symmetry implies that all coefficients with m==0 are purely real-valued. # The symmetry implies that all coefficients with m==0 are purely real-valued. # The a_lm are stored in a 1D complex-valued array, in the following order: # a_(0,0), a(1,0), ..., a_(lmax,0), a(1,1), a(2,1), ... a(lmax,1), ..., a(lmax, mmax) # number of required a_lm coefficients # number of required a_lm coefficients nalm = ((mmax+1)*(mmax+2))//2 + (mmax+1)*(lmax-mmax) nalm = ((mmax+1)*(mmax+2))//2 + (mmax+1)*(lmax-mmax) # number of real-valued random numbers to draw nalm_r = nalm*2-lmax-1 # get random a_lm # get random a_lm alm = np.random.uniform(-1., 1., nalm) + 1j*np.random.uniform(-1., 1., nalm) alm = np.random.uniform(-1., 1., nalm) + 1j*np.random.uniform(-1., 1., nalm) # make a_lm with m==0 real-valued # make a_lm with m==0 real-valued ... ...
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!