符号计算Sympy

胡友俊 2018-11-28 05:05

用了mathematica 与 maxima 很多年， 这几年开始用sympy, 感觉比前两个都好用。 这些工具对我来说仅是一个积分微分表， 遇到不熟悉的微分积分，直接在里头验证， 不用另查资料了， 方便快捷， 而且结果最可靠。 以下是我常用的几个命令： In : from sympy import * In :x=symbols('x') #defin ...

3887 次阅读|没有评论

托卡马卡磁力线可视化

热度 2 胡友俊 2018-7-27 08:20

下面几副图画了不同安全因子磁面上的磁力线. 画图的源代码见: p.sh and torus.sketch . Fig.1 安全因子q=1/10磁面上的磁力线（托卡马克里没这么小的安全因子）. Fig.2 安全因子q=1/2磁面上的磁力线. Fig.3 安全因子q=1磁面 ...

4354 次阅读|4 个评论 热度 2

PETSc并行计算库

热度 2 胡友俊 2018-7-24 10:38

PETSc 在好几个托卡马克大规模模拟程序（如XGC, GTC)中用到， 变得越来越流行。 PETSc的官方文档里说学习PETSc有很陡的学习曲线，一下子吓走很多人。 PETSc库的一个主要目的是让用户测试各种并行算法变得更容易。 只要懂基本的并行概念， 做过基本的偏微分数值解， 用PETSc一点不难。 我花了一周时间， 读了PET ...

12187 次阅读|4 个评论 热度 2

丘成桐: 花了真功夫，走的冤枉路都不会冤枉

热度 3 胡友俊 2017-3-8 16:25

杨振宁曾经说过 ：“温度场论是很有意思的。这是一个比较深的题目。文章很多，不过我没有去研究它。对于这一类问题，我的一般态度是这样的，如果我决定去研究它，我就一定从头做起，而且不先去看别人的文章。做了一段时间，如果有困难，再去看别人的文章。这样才能很好吸收别人的东西。1959年左右，我和李政道要讨论W ...

3284 次阅读|3 个评论 热度 3

粒子模拟笔记

胡友俊 2016-3-29 11:37

我写的关于粒子模拟的笔记 http://theory.ipp.ac.cn/~yj/research_notes/particle_simulation/index.html 其中有关于Landau damping的模拟结果与解析结果的比较，二者吻合得很好。

4000 次阅读|没有评论

气球模gif动画

胡友俊 2014-11-16 22:54

做了一个气球模的gif动画，如下图： 这是一个m=4的气球模，气球模的特征是低场侧的振幅大于高场侧，动画中低场侧的振幅是高场侧的两倍。 下图是一个有径向相差的气球模结构：

4506 次阅读|没有评论

托卡马克粒子导心轨道动画

热度 1 胡友俊 2014-11-6 16:20

做了个导心运动的三维动画，如下图 图中画的是能量40keV的俘获氘离子在EAST平衡位型(#38300@3.9s)下的运动，可以看到俘获粒子的环向进动

4554 次阅读|4 个评论 热度 1

一篇关于EAST上Alfven波的论文

胡友俊 2014-5-16 15:01

最近写了一篇关于EAST上Alfven本征模的论文， 投稿到PoP，整个过程如下: 这是初稿 595726_0_merged_1394498476-1.pdf 审稿人的意见: ------------------------------------------------------------------------------ Manuscript #POP43137: Referee (Remarks): In general this is a careful, well organized ...

3223 次阅读|没有评论

本页有 2 篇博文因作者的隐私设置或未通过审核而隐藏