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

Bugfix 1.2.1 #75

Merged
merged 10 commits into from
Apr 14, 2020
Merged

Bugfix 1.2.1 #75

merged 10 commits into from
Apr 14, 2020

Conversation

MuellerSeb
Copy link
Member

@MuellerSeb MuellerSeb commented Mar 22, 2020

  • ModuleNotFoundError not present in py35.
  • Fixing Cressie-Bug Bug: Cressie variogram estimator doubles variance #76
  • Adding analytical formula for integral scales of rational and stable model
  • remove prange from IncomprRandMeth summators to prevent errors on Win and macOS

@MuellerSeb MuellerSeb self-assigned this Mar 22, 2020
@MuellerSeb MuellerSeb added the bug Something isn't working label Mar 22, 2020
@MuellerSeb MuellerSeb requested a review from LSchueler March 22, 2020 19:48
Copy link
Member

@LSchueler LSchueler left a comment

Choose a reason for hiding this comment

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

Your good to go

@@ -325,6 +333,9 @@ def correlation(self, r):
r = np.array(np.abs(r), dtype=np.double)
return np.exp(-np.power(r / self.len_scale, self.alpha))

def calc_integral_scale(self): # noqa: D102
return self.len_scale * sps.gamma(1.0 + 1.0 / self.alpha)
Copy link
Member

Choose a reason for hiding this comment

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

:-) That looks a lot faster

Copy link
Member Author

Choose a reason for hiding this comment

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

It is!!

@@ -77,7 +77,7 @@ cdef inline void normalization_cressie(
if counts[i] == 0:
counts[i] = 1
variogram[i] = (
(1./counts[i] * variogram[i])**4 /
0.5 * (1./counts[i] * variogram[i])**4 /
Copy link
Member

Choose a reason for hiding this comment

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

Thanks for that!

Copy link
Member

@LSchueler LSchueler left a comment

Choose a reason for hiding this comment

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

Approved!

@MuellerSeb
Copy link
Member Author

MuellerSeb commented Apr 8, 2020

@LSchueler should we as well change these lines to use range:

for i in prange(X_len, nogil=True):

for i in prange(X_len, nogil=True):

to get rid of the openmp bug on windows and macOS?

@MuellerSeb MuellerSeb added Documentation enhancement New feature or request labels Apr 14, 2020
@MuellerSeb MuellerSeb merged commit 2c26d74 into master Apr 14, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working Documentation enhancement New feature or request
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants