Giter VIP home page Giter VIP logo

pysdm-examples's People

Contributors

abulenok avatar claresinger avatar edejong-caltech avatar imdula avatar jb-mackay avatar piotrbartman avatar sajjadazimi avatar slayoo avatar tangwhiap avatar trontrytel avatar

Watchers

avatar

Forkers

rxward

pysdm-examples's Issues

Update compressed film parameters

Hi @slayoo,

Big picture: I'm trying to do MCMC to find the optimal parameters values for the compressed film model.

Small picture: I'm trying just to make a plot where I have two curves for the surface tension given different parameter values in the compressed film framework (delta_min and sgm_org). You can see in this notebook that I update the parameters, but the returned surface tension is unchanged. Specifically I have narrowed it down to this:

compressed_film_Ovadnevaite.sgm_org = 40.1 * si.mN / si.m
compressed_film_Ovadnevaite.delta_min = 0.30 * si.nm
formulae_1 = Formulae(surface_tension='CompressedFilm_Ovadnevaite')

compressed_film_Ovadnevaite.sgm_org = 28.5 * si.mN / si.m
compressed_film_Ovadnevaite.delta_min =  0.52 * si.nm
formulae_2 = Formulae(surface_tension='CompressedFilm_Ovadnevaite')

print(formulae_1 == formulae_2)
print(formulae_1.surface_tension == formulae_2.surface_tension)
print(formulae_1.surface_tension.sigma == formulae_2.surface_tension.sigma)

False
True
True

So the formulae are updating, but the formulae.surface_tension are not. How can I get around this? Thanks so much for your help!

Error with Ruehl comp film implementation

Error in PySDM Ruehl surface tension implementation in the step where f_surf is found iteratively.
@slayoo I'm not familiar with this numba error. Any help?

Cannot capture the non-constant value associated with variable 'A_iso' in a function that will escape.

File "../../../PySDM/PySDM/physics/surface_tension/compressed_film_Ruehl.py", line 32:
    def sigma(T, v_wet, v_dry, f_org):
        <source elided>
        # solve implicitly for fraction of organic at surface
        f = lambda f_surf: Cb_iso*(1-f_surf) - C0*np.exp(((A0**2 - (A_iso/f_surf)**2)*m_sigma*sci.N_A)/(2*sci.R*T))
        ^

https://github.com/claresinger/PySDM-examples/blob/cf_ruehl/PySDM_examples/Singer/fig_1.ipynb

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.