multidworlds / eigenspectra

Eigenspectra
MIT License
0 stars 5 forks source link

Clean up the master jupyter notebook so it has clear, worked examples for hot spot and quadrant maps. #43

Open meganmansfield opened 5 years ago

meganmansfield commented 5 years ago

Hey everyone, I've been making progress on this task and have a few remaining questions. First, the in-progress updated jupyter notebook is titled master_cleaned.ipynb. Most of it is complete, but here are my remaining questions:

  1. Arthur - there is a blank section right now where the code to make your fancy looking plots goes. I wasn't sure where to find this code, so I just left it blank for now. Can you fill it in?
  2. Everett - I think you added the function "all_sph_degrees" at one point, but I'm not sure if we need it in the worked example. What do you think?
  3. Everett - When I run run_higher_sph_harm.run_lc_noise_and_fit for a bunch of runs with different degrees of spherical harmonics, I get this error: "OSError: [Errno 24] Too many open files". It works when I try and run only one set of models at a time with one spherical harmonic degree, but then when I try and test lots of different degrees it quits halfway through. Do you know how to fix this?
meganmansfield commented 5 years ago

One more general note, for those of you that have been working on this code: I made some edits to the function find_groups to allow it to work with a user-defined directory instead of just for the one specific hot spot case or whatever. However, now when I run through the code, I can't seem to get it to work for the quadrant model case. In step 1 of the cleaned script (generate lightcurves with starry) it is working and making lightcurves from a planet map with 4 quadrants of different spectra. However, in step 4 (use k means clustering to find groups) it looks like it's still operating on the spectra from the hot spot model? I'm guessing I just misplaced a word somewhere so it's still looking at the old directory for some reason, but for the life of me I can't find it. Can any of you guys figure out what's wrong? Or maybe just try running the quadrant model section (Section II) and see what happens?