7
7
8
8
< meta name ="viewport " content ="width=device-width, initial-scale=1.0 " />
9
9
10
- < title > mutis.signal.methods — MUTIS documentation</ title >
10
+ < title > mutis.signal.functions — MUTIS documentation</ title >
11
11
12
12
13
13
37
37
< script src ="../../../_static/jquery.js "> </ script >
38
38
< script src ="../../../_static/underscore.js "> </ script >
39
39
< script src ="../../../_static/doctools.js "> </ script >
40
+ < script src ="../../../_static/copybutton.js "> </ script >
40
41
< script crossorigin ="anonymous " integrity ="sha256-Ae2Vz/4ePdIu6ZyI/5ZGsYnb+m0JlOmKPjt6XZ9JJkA= " src ="https://cdnjs.cloudflare.com/ajax/libs/require.js/2.3.4/require.min.js "> </ script >
41
42
42
43
< script type ="text/javascript " src ="../../../_static/js/theme.js "> </ script >
149
150
150
151
< li > < a href ="../../index.html "> Module code</ a > »</ li >
151
152
152
- < li > mutis.signal.methods </ li >
153
+ < li > mutis.signal.functions </ li >
153
154
154
155
155
156
< li class ="wy-breadcrumbs-aside ">
164
165
< div role ="main " class ="document " itemscope ="itemscope " itemtype ="http://schema.org/Article ">
165
166
< div itemprop ="articleBody ">
166
167
167
- < h1 > Source code for mutis.signal.methods </ h1 > < div class ="highlight "> < pre >
168
+ < h1 > Source code for mutis.signal.functions </ h1 > < div class ="highlight "> < pre >
168
169
< span > </ span > < span class ="c1 "> # Licensed under a 3-clause BSD style license - see LICENSE</ span >
169
170
< span class ="sd "> """Methods for synthetic generation of light curves."""</ span >
170
171
171
172
< span class ="kn "> import</ span > < span class ="nn "> logging</ span >
172
173
174
+ < span class ="kn "> import</ span > < span class ="nn "> nfft</ span >
173
175
< span class ="kn "> import</ span > < span class ="nn "> numpy</ span > < span class ="k "> as</ span > < span class ="nn "> np</ span >
174
176
< span class ="kn "> import</ span > < span class ="nn "> scipy.signal</ span > < span class ="k "> as</ span > < span class ="nn "> scipy_signal</ span >
175
- < span class ="kn "> import</ span > < span class ="nn "> nfft</ span >
176
177
177
178
< span class ="n "> __all__</ span > < span class ="o "> =</ span > < span class ="p "> [</ span >
178
179
< span class ="s2 "> "lc_gen_samp"</ span > < span class ="p "> ,</ span >
@@ -255,6 +256,8 @@ <h1>Source code for mutis.signal.methods</h1><div class="highlight"><pre>
255
256
256
257
257
258
< div class ="viewcode-block " id ="lc_gen_psd_lombscargle "> < a class ="viewcode-back " href ="../../../api/mutis.signal.lc_gen_psd_lombscargle.html#mutis.signal.lc_gen_psd_lombscargle "> [docs]</ a > < span class ="k "> def</ span > < span class ="nf "> lc_gen_psd_lombscargle</ span > < span class ="p "> (</ span > < span class ="n "> times</ span > < span class ="p "> ,</ span > < span class ="n "> signs</ span > < span class ="p "> ):</ span >
259
+ < span class ="sd "> """Description goes here."""</ span >
260
+
258
261
< span class ="k "> if</ span > < span class ="n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> size</ span > < span class ="o "> %</ span > < span class ="mi "> 2</ span > < span class ="o "> !=</ span > < span class ="mi "> 0</ span > < span class ="p "> :</ span >
259
262
< span class ="n "> sigp</ span > < span class ="o "> =</ span > < span class ="n "> signs</ span > < span class ="p "> [:</ span > < span class ="o "> -</ span > < span class ="mi "> 1</ span > < span class ="p "> ]</ span >
260
263
< span class ="n "> tp</ span > < span class ="o "> =</ span > < span class ="n "> times</ span > < span class ="p "> [:</ span > < span class ="o "> -</ span > < span class ="mi "> 1</ span > < span class ="p "> ]</ span >
@@ -276,7 +279,7 @@ <h1>Source code for mutis.signal.methods</h1><div class="highlight"><pre>
276
279
< span class ="n "> s2</ span > < span class ="o "> =</ span > < span class ="n "> nfft</ span > < span class ="o "> .</ span > < span class ="n "> nfft</ span > < span class ="p "> ((</ span > < span class ="n "> times</ span > < span class ="o "> -</ span > < span class ="p "> (</ span > < span class ="n "> times</ span > < span class ="o "> .</ span > < span class ="n "> max</ span > < span class ="p "> ()</ span > < span class ="o "> +</ span > < span class ="n "> times</ span > < span class ="o "> .</ span > < span class ="n "> min</ span > < span class ="p "> ())</ span > < span class ="o "> /</ span > < span class ="mi "> 2</ span > < span class ="p "> )</ span > < span class ="o "> /</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> ptp</ span > < span class ="p "> (</ span > < span class ="n "> times</ span > < span class ="p "> ),</ span > < span class ="n "> fft2</ span > < span class ="p "> ,</ span > < span class ="n "> N</ span > < span class ="p "> ,</ span > < span class ="n "> use_fft</ span > < span class ="o "> =</ span > < span class ="kc "> True</ span > < span class ="p "> )</ span > < span class ="o "> /</ span > < span class ="n "> N</ span >
277
280
278
281
< span class ="c1 "> # fix small deviations</ span >
279
- < span class ="n "> a</ span > < span class ="o "> =</ span > < span class ="p " > ( </ span > < span class =" n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> ()</ span > < span class ="o "> /</ span > < span class ="n "> s2</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> () )</ span >
282
+ < span class ="n "> a</ span > < span class ="o "> =</ span > < span class ="n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> ()</ span > < span class ="o "> /</ span > < span class ="n "> s2</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> ()</ span >
280
283
< span class ="n "> b</ span > < span class ="o "> =</ span > < span class ="n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> mean</ span > < span class ="p "> ()</ span > < span class ="o "> -</ span > < span class ="n "> a</ span > < span class ="o "> *</ span > < span class ="n "> s2</ span > < span class ="o "> .</ span > < span class ="n "> mean</ span > < span class ="p "> ()</ span >
281
284
< span class ="n "> s2</ span > < span class ="o "> =</ span > < span class ="n "> a</ span > < span class ="o "> *</ span > < span class ="n "> s2</ span > < span class ="o "> +</ span > < span class ="n "> b</ span >
282
285
@@ -289,7 +292,7 @@ <h1>Source code for mutis.signal.methods</h1><div class="highlight"><pre>
289
292
290
293
< span class ="n "> nft</ span > < span class ="o "> =</ span > < span class ="n "> nfft</ span > < span class ="o "> .</ span > < span class ="n "> nfft_adjoint</ span > < span class ="p "> ((</ span > < span class ="n "> times</ span > < span class ="o "> -</ span > < span class ="p "> (</ span > < span class ="n "> times</ span > < span class ="o "> .</ span > < span class ="n "> max</ span > < span class ="p "> ()</ span > < span class ="o "> +</ span > < span class ="n "> times</ span > < span class ="o "> .</ span > < span class ="n "> min</ span > < span class ="p "> ())</ span > < span class ="o "> /</ span > < span class ="mi "> 2</ span > < span class ="p "> )</ span > < span class ="o "> /</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> ptp</ span > < span class ="p "> (</ span > < span class ="n "> times</ span > < span class ="p "> ),</ span > < span class ="n "> signs</ span > < span class ="p "> ,</ span > < span class ="n "> N</ span > < span class ="p "> ,</ span > < span class ="n "> use_fft</ span > < span class ="o "> =</ span > < span class ="kc "> True</ span > < span class ="p "> )</ span >
291
294
292
- < span class ="c1 "> # build random phase to get real signal: </ span >
295
+ < span class ="c1 "> # build random phase to get real signal</ span >
293
296
< span class ="n "> phase</ span > < span class ="o "> =</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> random</ span > < span class ="o "> .</ span > < span class ="n "> random</ span > < span class ="p "> (</ span > < span class ="n "> N</ span > < span class ="o "> //</ span > < span class ="mi "> 2</ span > < span class ="p "> )</ span >
294
297
< span class ="n "> phase</ span > < span class ="o "> =</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> concatenate</ span > < span class ="p "> ((</ span > < span class ="o "> -</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> flip</ span > < span class ="p "> (</ span > < span class ="n "> phase</ span > < span class ="p "> ),</ span > < span class ="p "> [</ span > < span class ="mi "> 0</ span > < span class ="p "> ],</ span > < span class ="n "> phase</ span > < span class ="p "> [:</ span > < span class ="o "> -</ span > < span class ="mi "> 1</ span > < span class ="p "> ]))</ span >
295
298
@@ -298,7 +301,7 @@ <h1>Source code for mutis.signal.methods</h1><div class="highlight"><pre>
298
301
< span class ="n "> s2</ span > < span class ="o "> =</ span > < span class ="n "> np</ span > < span class ="o "> .</ span > < span class ="n "> real</ span > < span class ="p "> (</ span > < span class ="n "> s2</ span > < span class ="p "> )</ span > < span class ="c1 "> # np.real to fix small imaginary part from numerical error</ span >
299
302
300
303
< span class ="c1 "> # fix small mean, std difference from numerical error</ span >
301
- < span class ="n "> a</ span > < span class ="o "> =</ span > < span class ="p " > ( </ span > < span class =" n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> ()</ span > < span class ="o "> /</ span > < span class ="n "> s2</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> () )</ span >
304
+ < span class ="n "> a</ span > < span class ="o "> =</ span > < span class ="n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> ()</ span > < span class ="o "> /</ span > < span class ="n "> s2</ span > < span class ="o "> .</ span > < span class ="n "> std</ span > < span class ="p "> ()</ span >
302
305
< span class ="n "> b</ span > < span class ="o "> =</ span > < span class ="n "> signs</ span > < span class ="o "> .</ span > < span class ="n "> mean</ span > < span class ="p "> ()</ span > < span class ="o "> -</ span > < span class ="n "> a</ span > < span class ="o "> *</ span > < span class ="n "> s2</ span > < span class ="o "> .</ span > < span class ="n "> mean</ span > < span class ="p "> ()</ span >
303
306
< span class ="n "> s2</ span > < span class ="o "> =</ span > < span class ="n "> a</ span > < span class ="o "> *</ span > < span class ="n "> s2</ span > < span class ="o "> +</ span > < span class ="n "> b</ span >
304
307
0 commit comments