原始:膨胀方程求解器

primordial的Python项目详细描述


一、二、二、二、二、三、三、三、三、三、三、四、四、四、四、四、四、四、四、四、四、四、四、四、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、https://github.com/williamjameshandley/primordial
:文档:http://primordial.readthedocs.io/

……图片:https://travis-ci.org/williamjameshandley/primordial.svg?branch=master
:目标:https://travis ci.org/williamjameshandley/primordial
:alt:build status
。图像::https://codecov.io/gh/williamjameshandley/primordial/branch/master/graph/badge.svg
:目标:https://codecov.io/gh/williamjameshandley/primordial
:alt:测试覆盖状态
。图片::https://badge.fury.io/py/primordial.svg
:目标:https://badge.fury.io/py/primordial
:alt:pypi location
…图片::https://readthedocs.org/projects/primordial/badge/?version=latest
:目标:https://primordial.readthedocs.io/en/latest/?徽章=最新
:alt:文档状态
…图片::https://zenodo.org/badge/150298714.svg
:target:https://zenodo.org/badge/latestdoi/150298714
:alt:permanent doi for this release





===

``primordial``是一个用于解决宇宙膨胀的python包方程组。

它非常处于beta阶段,目前正在为研究目的而构建。




示例用法
==


代码::python

import numpy
import matplotlib.pyplot as plt
from primordial.solver import solve
from primordial.equations.inflation戡potentials import chaoticpotential
from primordial.equations.t.inflation import equations,kd_初始条件
来自原始方程。事件导入通货膨胀,崩溃


fig,ax=plt。子块(3,sharex=true)
对于k in[-1,0,+1]:
m=1
v=chaoticpotential(m)
方程组(k,v)

事件=[通胀(方程式),记录通胀进出口
通胀(方程式,-1,终端=真),通胀退出时停止
崩溃(方程式,终端机=真)]如果宇宙停止膨胀则停止

phi p=23
t_p=1e-5
ic=kd_初始条件(t_p,n_p,phi_p)
t=numpy.logspace(-5,10,1e6)

sol=solve(方程式,ic,t_eval=t,事件=事件)

ax[0].plot(sol.n(t),sol.phi(t))
ax[0].set戥label(r'$\phi$')

ax[1].set戥yscale('log')
ax[1].set戥label(r'$h$')

ax[2].plot(sol.n(t),1/(sol.h(t)*numpy.exp(sol.n(t)))
ax[2].set_yscale('log')
ax[2].set_ylabel(r'$1/ah$')


ax[-1].set_xlabel('$n$')



plot mode function evolution
--------
。代码::python

import numpy
import matplotlib.pyplot as plt
from primordial.solver import solve
from primordial.equations.inflation\potentials import chaoticpotential
from primordial.equations.t.mukhanov_sasaki import equations,kd_initial_conditions
从primordial.equations.events导入膨胀、折叠、modeexit


fig,axes=plt。子块(3,sharex=true)
对于ax,k in zip(axes,[-1,0,+1]):
ax2=ax.twinx()
m=1
v=chaoticpotential(m)
k=100
方程=方程式(k,v,k)

ps展开
模式退出(方程式,+1,终端=真,值=1e1*k)模式退出时停止
]



n_p=-1.5
phi p=23
t愀p=1e-5
ic=kd愀u初始条件(t愀p,n愀p,phi p)
t=numpy.logspace(-5,10,1E6)

sol=solsolsolsolsolsolsolsolsolsolsolsolsolsolsolsolsolsolsol n(t)
ax.plot(n,sol.r1(t),‘k-’)

ax2.plot(n,—numpy.log(sol.h(t))-n,‘b-’)






















<')
ax.text(0.9,0.9,r'$k=%i$'%k,transform=ax.transaxes)


轴[-1].set_xlabel('$n$')






……|图像0图像::https://raw.githubusercontent.com/williamjameshandley/primordial/master/figures/background.png
…| image1 image::https://raw.githubusercontent.com/williamjameshandley/primordial/master/figures/ms.png


待办事项列表
==
最终希望将此提交给"joss<;https://joss.theoj.org/>;`。在此之前需要做的事情有:

宇宙学
——
*慢滚初始条件
*mukhanov-sazaki进化:数学:`n`
*add:math:`\eta`作为自变量
*add:math:`\phi`作为自变量

code
——
*documentation
*测试
*100%覆盖
*插值
*宇宙学

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java如何向第二个组合框模型项添加数组值从第一个组合框所选项获取数组名称?   使用Java与WebSphere的SSL握手错误   eclipse线程“main”Java中的第一个Java程序异常。lang.NoClassDefFoundError   java将Javafx应用作为Web应用移植的最佳方式   IDE的java右JVM文件夹   java如何在基本适配器中停止文本到语音   java If block使用substring和equals方法以假值执行   在本例中,如何在java中返回多个值?   java第二个主类在maven构建期间覆盖第一个主类   如何在java中设置运行时ArrayList的泛型类型?   java从主机读取文件