Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

The final carsus merge #786

Merged
merged 92 commits into from
Aug 31, 2017
Merged

The final carsus merge #786

merged 92 commits into from
Aug 31, 2017

Conversation

yeganer
Copy link
Contributor

@yeganer yeganer commented Aug 25, 2017

This PR replaces #770

@@ -1,7 +1,8 @@
import os
from matplotlib import pyplot as plt
from matplotlib import colors
from tardis import atomic, plasma_array, util
from tardis import plasma_array, util
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I should fix that, right?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

probably

@@ -1,7 +1,8 @@
import os

from matplotlib import colors
from tardis import atomic, plasma_array, util
from tardis import plasma_array, util
Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Same as above, there is no plasma_array

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

yep

@@ -162,14 +162,14 @@ def _main_nlte_calculation(
(number_of_levels**2, len(t_electrons)))
r_lu_matrix_reshaped[r_lu_index] = B_lus[np.newaxis].T * \
j_blues_filtered * beta_sobolevs_filtered
if atomic_data.has_collision_data:
if atomic_data.collision_data is not None:
if previous_electron_densities is None:
collision_matrix = r_ul_matrix.copy()
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

collision_matrix = np.zeros_like(r_ul_matrix)

initial_t_inner + initial_t_rad + helium_treatment)


def idfn(fixture_value):
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

can you quickly document this

@yeganer
Copy link
Contributor Author

yeganer commented Aug 28, 2017

My verification show that the macroatom scheme does not work right now. These are three plots with current master and this PR for macroatom, downbranch and scatter.
As you can see, for macroatom there are clear differences while downbranch and scatter are the same.
Running my NLTE tests with macroatom results in an infinite loop which further points towards a problem with the macro atom. I'll have a closer look tomorrow.

example
example_downbranch
example_scatter

@yeganer
Copy link
Contributor Author

yeganer commented Aug 28, 2017

Found the bug, nevermind.

@yeganer
Copy link
Contributor Author

yeganer commented Aug 29, 2017

After fixing the bug, all tests involving kurucz_cd23_chianti_h_he.h5 produce near identical results, for example here a setup from @talytha :
tomography

However my NLTE test setup (with no NLTE species) does have a slightly different spectrum:
nlte

@wkerzendorf How should we proceed?

That means we don't need to set indexes in tardis.
Additionally all 'copys' were removed as there does not seem to be a
special need for them.
As a result, there will be SettingWithCopy warnings but my tests showed
they are not problematic.
Additionally: Internally rename 'tardis-refdata' config variable to
'tardis_refdata'.
Also we only close if we opened the HDFStore
This should show up if any property is not found during testing
This is to see the travis error message
@yeganer
Copy link
Contributor Author

yeganer commented Aug 31, 2017

@wkerzendorf Tests are passing. Looks like this might be ready to merge.

@wkerzendorf wkerzendorf merged commit 4dfe08b into tardis-sn:master Aug 31, 2017
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

4 participants