原始:膨胀方程求解器

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与服务器建立TCP连接后,为什么我看不到服务器发送的消息?   java解析数据到Java8中的2D数组   java Lifecycle OnLifecycleEven已弃用   Java Swing JDBC和MySQL ClientServer LAN设置   java滚动背景边线条的外观   转换日期时发生java异常   对Cucumber 4.7的java IntelliJ支持   mongodb Mongo和java查找查询以及batchsize   java为充气菜单添加自定义图标   java对象到ArrayList   java作为JPackage安装的一部分,我如何在linux中使本机linux应用程序可执行   数字的java JTextfield验证?   java更改Eclipse中默认的Maven存储库,以下载所需的TinyB库   java如何使用Jena从RDF中选择具有特定属性的主题?   java Selenium WD无法在cheaptickets上的弹出窗口中找到webelement。在里面