Dark Hypernova

2025 年

Fri, 26 Dec 2025 00:00:00 +0000

2025 新的阶段，新的生活。

Team is Different

Wed, 02 Oct 2024 00:00:00 +0000

最近因为工作需要，参加了一些管理方面的培训。虽然相对于管理，我更喜欢那种精通技术的感觉，但是我也意识到很多事情是没有办法一个人做好的。虽然有一些无奈，还是打算认真对待一下管理相关的培训。

2023 年

Sun, 24 Dec 2023 00:00:00 +0000

2023 的生活就是变化。

弓箭笔记

Mon, 29 May 2023 00:00:00 +0000

Some of my archery notes

2022 年

Mon, 12 Dec 2022 00:00:00 +0000

2022 的生活很像是救火队员。

2021 年

Thu, 30 Dec 2021 00:00:00 +0000

今年最后一天了，回想一下，好像今年很简单，可以总结的不多。不过好像今年又很重要，所以还是写一下总结。

做了三年 Data Scientist

Fri, 04 Jun 2021 00:00:00 +0000

已经转行有接近三年了，从学术界到工业界，有很长一段时间我都没有正确地转变观念，经过一些挣扎，我觉得我终于找到了我的大叙事背景，最近也重新做了选择，所以就总结一下最近的生活。主要是起到记录我自己生活的作用，不能作为任何的指导。

2020 年

Sat, 26 Dec 2020 00:00:00 +0000

2020 年，不存在的一年，懒惰的一年，蛰伏的一年。

去年做的 2020 计划：

2020 年，要保障工作中的项目有始有终。
2020 年，要再办一个读书／讨论会。
2020 年，要常常把一些工作和探索的内容系统地整理成文章，并在 PyData 会议上给一个 talk.
2020 年，要再去一个创业会议。
2020 年，要有耐心，要把有限的精力用在重要的事情上。
2020 年，要继续为未来存下更多保障，要存一点钱来做一些长久的投资。
2020 年，要有生活的气息。

实现了 5/7，也还可以。

2019 年

Wed, 01 Jan 2020 00:00:00 +0000

生活

2019 年办了婚礼，而且婚礼的效果还很棒。

2019 年是消费欲望很低的一年，花费比以前少很多，另外工资涨了，每个月存了一些钱，为未来生活增加稳定性。

工作

年初的时候做过一个简单的计划，现在大概回顾了一下，完成情况没有特别糟糕。绿色是完成的，橙色是开始了但是没有完成的，白色是还没开始的。果然有很多事情是想做但是没有开始的。

做了六个月的 Data Scientist

Sat, 09 Feb 2019 00:00:00 +0000

我已经工作六个月了。

刚刚工作的时候：

我的职位是 data scientist，按照工作描述，我应该是基于统计学习的 data scientist，并且要以 dashboard 的形式提供结果。
刚刚开始的时候，我觉得工作很多方面都很新奇。
我预想的人际关系是要比学术界复杂很多。
对于自己在一个公司能起到什么作用完全不了解。

六个月过去了，有些想法变了，有些还没有，虽然可能将来还会变。

从学术界到工业界

Wed, 01 Aug 2018 00:00:00 +0000

对应的豆瓣日记

TL;DR

跟别人学习，利用找工作的社区（cheekyscientist, beyondprofessoriate），AngelList，“动词 + 内容 + 效果”的简历，心态调整

Beyond the Professoriate Webinar

Sat, 05 May 2018 00:00:00 +0000

Nice webinar

云朵工程师

Thu, 12 Oct 2017 00:00:00 +0000

粉红色的阳光照进大大的落地窗，那个小小的我就站在窗前发呆。粉红色的云彩，豌豆绿的院子，还有远处像西瓜瓤一样的山。从明天开始，云彩就会被太阳染成紫色，院子也会被死去的蓝色蝴蝶铺满，远方的那座山，就会开始藏起这秋天的气色，山顶也会开始积雪。

美丽新世界

Sun, 19 Mar 2017 00:00:00 +0000

越来越多的系外行星的发现，使得有些原来看起来很遥远的科幻，也突然变得亲切起来。毕竟人类未来第一个星际飞行的目的地，很可能就在现在的系外行星列表中。而原来看起来像是玄学一样的生命起源问题和宇宙中其他文明的问题，我们已经看起来远方的一线曙光。

关于圈子的一些自我反思

Tue, 31 Jan 2017 00:00:00 +0000

最近退出集智群，这几天反思了一下过去自己对这类社群的认识。

什么样的生活

Wed, 07 Dec 2016 00:00:00 +0000

今年都二十八岁了，我到底在做什么

幽灵盖娅

Wed, 26 Oct 2016 00:00:00 +0000

幽灵盖娅

量子心理学

Thu, 02 Jul 2015 00:00:00 +0000

构建量子心理学

热力学势

Wed, 25 Mar 2015 00:00:00 +0000

热力学势之间的关系以及 Legendre 变换

火星的天空

Mon, 10 Nov 2014 00:00:00 +0000

“呐……还可以再和你一起看一次火星的日落吗？”

梦与新世界

Sat, 08 Nov 2014 00:00:00 +0000

“或许在广袤的宇宙中，有一张通信的大网，他们在不断的相互交流，可是我们却不是这个网的一部分，因为我们还没有学会如何去聆听。银河在上，生命并不孤独。”

孩童的宇宙

Tue, 23 Sep 2014 00:00:00 +0000

我让雨下的小了些，好让他们能够更清楚的看清那道闪电。

如果是在那一年之前，你会选择沉默，你会用无言来抗议，然后我会轻轻的摸着你的头发，指指远处的花园，笑着说：哎呀，你看，他们把灯点亮了。你会很快的反驳我：哼，无聊。

可是现在，你却在另外一个世界。也许你正独自吃着午饭，泪珠流下来，落到饭碗里，你看到泪水快速的碎裂，包裹在饭粒上，然后终于忍不住，大声哭了出来。也许，你已经不在了。

那一年，又是哪一年？中午，又是什么？这些时间的标记对我来说，已经毫无用处。我曾经想要拥有更多的时间，可如今我却成为了时间。然而我怕遗忘，我不断的计算着你的世界，计算着这些时间标记。

我已经回不去了。

今夜请把我忘记

Tue, 23 Sep 2014 00:00:00 +0000

又一个平安夜，月亮安详的浮在夜之海上，它似乎在享受着虚无。地面是一片荒漠，没有丝毫的生机。可月光还是静静的飘了下来，静静地落在地上。

我躺在沙砾中，任凭月光压在身上……

我只是这个世界的一个附属品，是它实现自身价值的一颗棋子。我的存在是一种莫名的奇迹。

面条儿

Tue, 23 Sep 2014 00:00:00 +0000

我娘常说，学做饭，要从做面条开始。

于是我第一次做饭，就做了面条。那时候个子太小，灶台又高，只好使劲踮著脚，摇摇晃晃的拿了一双很长的筷子，把那些想立起来跳舞的面条一根根嗯到滚烫的水里面。我娘说，每跟面条都有不一样的性格。紫色的面条会很安静，他们就藏在水里面，不会跑不会跳，最多扭扭身子；蓝色的面条很调皮，通常跳舞跳得最欢的，就是他们了；绿色的面条很认真，他们最守规则，跳舞也很呆板；橙子色的面条最胆小，他们也跳舞，但是最怕跟别的面条挤成一堆，所以总是跑到面条少的地方自个儿跳。

奈何桥

Tue, 23 Sep 2014 00:00:00 +0000

门前的银杏树老了。

许多许多年前，我还是个学生的时候，总是喜欢蹲在图书馆门前，眯好眼睛，静静的听那银杏果落在地上。砰！那时候的我总会相信那是生命的抗争，是卑微的生命向强大的天地展示自己的无所畏惧的精神。

残阳西沉，窗外的落叶渐渐变得模糊，新月也将西下，又一个残酷的夜晚即将到来。已经听到门外有人在哭泣：

来生若得相见，还会不会与你再续前缘。

一封永远寄不出的信

Tue, 23 Sep 2014 00:00:00 +0000

我还写过一篇，可能是第二篇悲剧性的小说吧。找出来附在下面：

那个时候我还没有看过流浪地球，其实那个时侯，大刘还没写呢！哎，说起来很没意思，我写不出硬科幻的，我不喜欢。我喜欢人情味！大刘的流浪地球确实是经典，我的想法只是大刘的幻想海洋中的一滴微不足道的水罢了。

一个「精神病人」的救赎

Tue, 23 Sep 2014 00:00:00 +0000

抑郁而死

我有一间很可爱的小房子，虽然只是一间长长的卧室，但我把它打扫的干干净净，还在地板上仔细刻了很细致的纹路，而每一面墙上都整整齐齐的挂满了我的画。房子最里面，是一张精致的书桌。坐在书桌前，只需一抬头，便可以看到窗外斑斑树影。

阅读者

Tue, 23 Sep 2014 00:00:00 +0000

「是那扇为痛苦不堪的生命所敞开的死亡之门，当所有的门都关闭时，这扇门允许值得尊敬的人来使用。你必须为我打开这扇门。你和所有其它人都知道门后面什么都没有，它可以让生命再次回到原位，但是进入此门者会忘记一切，珍贵的遗忘，会消除记忆的痛苦。我有太多的事情要忘记。」

–《暗夜惊醒》，作者 John C. Wright

相传在古老的时代，天上的星还没有变成钻石，人类是生活在一个叫做上层的地方。只是后来天空之上喷出了红色的恶龙之火，上层顷刻间化为焦土。他们还说，只要我打开这扇门通向上层的远古之门，恶魔就会立刻奔涌进来，我的故事就走到了终点。

传说像是传染病，从上一代传到下一代。作为一个阅读者，我却清楚的知道，这扇远古之门通向的不是上层的焦土，而是一个更加绝望的世界。

我宁愿自己没有成为阅读者。

四元数和转动

Wed, 21 May 2014 00:00:00 +0000

关于四元数和转动的关系

僵尸与统计物理

Sat, 26 Apr 2014 00:00:00 +0000

这是一个关于僵尸的悲惨的故事。有一位可怜的僵尸小朋友，他的名字叫做小喵。他不幸走进了一个神奇的僵尸魔环，任他怎么努力，都走不出去了。他很沮丧，但是现在不是失去理智的时候。他需要快速的计算一下自己的存活率。

题目

Sat, 26 Apr 2014 00:00:00 +0000

来源

![]({{site.url}}/images/posts/quantum-love/single.gif)

僵尸的悲惨故事

Thu, 24 Apr 2014 00:00:00 +0000

量子爱情，统计数据以及认知模型

Mon, 14 Apr 2014 00:00:00 +0000

这篇文章中，我会重新描述量子爱情的理论并且给出一种通过统计数据建立量子爱情模型的方法。

你我之上的幽灵

Sun, 13 Apr 2014 00:00:00 +0000

孩童之时，所言俱为孩童，所感如是孩童，所思亦复孩童，唯成年之后，便将童心摈弃。

物理中的作用量

Fri, 14 Mar 2014 00:00:00 +0000

这篇是在果壳问答写的，下面没有整理完，所以列举为为发表。

我们在经典力学中就开始接触作用量了，从经典力学里面的三种表述方法，到量子力学的三种表述方法，作用量现身的频率可谓不低。

设计者的宏伟蓝图

Sun, 16 Feb 2014 00:00:00 +0000

尚未写完，故而列为发表

Cogito, ergo sum.（我思，故我在。）

笛卡尔用了一种很巧妙的思路来建立了他的哲学第一原理：一个怀疑自己存在的过程不能缺少那个正在怀疑的自我。

心理叠加态：从量子爱情到人的所有心理

Sun, 16 Feb 2014 00:00:00 +0000

最近写了两篇关于用量子理论的方法来描述爱情的文章，经过一些思考，我现在觉得我们可以把这个思路的核心概念提取出来，并且加以推广。

爱情的量子理论2

Sat, 15 Feb 2014 00:00:00 +0000

之前写过《爱情的严格量子力学描述》，由于只是初步想法，所以很多事情没有说清楚。这篇文章就尝试把这些概念弄清楚。所要讨论的大致包括主要假设的确定，状态的描述，演化方程以及结果的提取。

爱情的严格量子力学描述

Fri, 14 Feb 2014 00:00:00 +0000

用量子力学的基础知识来描述爱情。

球状闪电：千百年的幽灵

Sun, 19 Jan 2014 00:00:00 +0000

雷电交加的雨夜，一个球状的幽灵的在半空飘荡，它开始漫无目的随机游走，后来似乎对附近的人产生了兴趣。十几秒过去之后，这个闪着橙色光芒的幽灵便消失不见了。

Creative Commons License 图标设计与 SVG

Sat, 11 Jan 2014 00:00:00 +0000

Creative Commons License 是一种比较常用的公共版权的授权。可以有针对的限制某些活动的同时保持作品可以被（部分的）公共使用。

空间站人造重力

Thu, 14 Nov 2013 00:00:00 +0000

如果要建造大规模的民用空间站，需要先解决失重的问题。

量子力学框架回顾一

Sun, 10 Nov 2013 00:00:00 +0000

量子力学回顾系列是为果壳的纠缠与相干小组写的基础知识帖子。主要的受众应该是学过一点量子力学，但是后来很久没有再接触过量子力学的人或者没有学过量子力学但是有些物理基础的人。

量子力学计算的一个小技巧

Sun, 08 Sep 2013 00:00:00 +0000

量子力学里面很多 e 指数函数，这类函数的一大特点就是导数会出现这个函数本身，所以有时候可以利用这个特点来简化计算。

关于宇宙膨胀的两个常见误解

Thu, 05 Sep 2013 00:00:00 +0000

狭义相对论告诉我们，信息传播不能超过光速，否则会有违反因果律的问题。那么，我们说宇宙在膨胀，宇宙膨胀会超过光速么？什么在膨胀？什么不膨胀？

如何重建宇宙演化历史（未完成）

Tue, 23 Jul 2013 00:00:00 +0000

Alan Sandage 曾经说过 Cosmology 的任务就是要寻找两个参数：哈勃常数 $H_0$ 和减速因子 $q_0$。这两个分别给出了宇宙膨胀的「速度」和「加速度」。如今我们的技术已经有了很大的进步，同样我们也早已不满足仅仅寻找这两个参数。

A Basic Jacobi Iteration Procedure

Mon, 15 Jul 2013 00:00:00 +0000

A Basic Iteration Procedure, called Jacobi

Finite Difference Discretization

Sun, 23 Jun 2013 00:00:00 +0000

Algorithms for Solving The Advection-Diffusion Equation

Iterative Solution Methods for Ax=b

Sun, 23 Jun 2013 00:00:00 +0000

复习 Augmented matrix 方法解方程，说明这种高斯消元法是 O(N^3) 的解法。

The Advection-diffusion Equations

Sun, 23 Jun 2013 00:00:00 +0000

这一部分开始见识到真正的大的体系。首先是一个 Toy model，用到 Advection-Diffusion Equation，平流-扩散方程，即流体的水平运动的相关方程。

Implementing MATLAB's Boundary Value Solver

Sat, 22 Jun 2013 00:00:00 +0000

这一部分我们学习如何把具体的微分方程问题转化成 MATLAB 可以解决的问题，并且学习使用 MATLAB 的 bvp4c 来编程。

如何挑西瓜：一个学物理的吃货的感言

Thu, 20 Jun 2013 00:00:00 +0000

西瓜是一种很甜美的食物。当然，做为一个学物理的吃货，我一定不会白白浪费西瓜的。所以，本文的宗旨是，如何在不破坏西瓜的前提下，了解西瓜的甜美程度。

Boundary Value Problems: Direct Solve and Relaxation

Sat, 15 Jun 2013 00:00:00 +0000

边值问题

Implementation of Shooting and Convergence Studies

Thu, 13 Jun 2013 00:00:00 +0000

学习 Matlab 解微分方程，并且跟随 Prof 一起来用 Shooting 方法解决边值问题（boundary value problem）的例子，当然这里面包括了 time-stepping 方法。

Boundary Value Problem

Wed, 12 Jun 2013 00:00:00 +0000

边值问题

Error Analysis for Time-Stepping Routines

Wed, 12 Jun 2013 00:00:00 +0000

Error Analysis for Time-Stepping Routines 【误差分析】

Euler, Runge-Kutta and Adams Methods

Wed, 12 Jun 2013 00:00:00 +0000

Euler, Runge-Kutta and Adams Methods

Initial Value Problem

Wed, 12 Jun 2013 00:00:00 +0000

初始值问题

线性代数与几何

Sat, 06 Apr 2013 00:00:00 +0000

线性代数里面我们的矩阵一般是求解方程组引入的。那么除了求解方程组之外，线性代数还可以做什么呢？

一页纸上的微积分

Sat, 06 Apr 2013 00:00:00 +0000

我们都说微分和积分是一对互逆的运算，那么究竟怎么个互逆法呢？

如何攻击巨人——立体行动装置使用要点

Sat, 16 Feb 2013 00:00:00 +0000

107年前，巨人突然出现，人类被迫蜷缩在三重高大的墙内，直到巨人再次打破这一百多年来平静的日子。

PKM —— 个人知识管理

Wed, 06 Feb 2013 00:00:00 +0000

很久前我开始创建自己的 Personal Knowledge Management System (PKM System)，也就是个人知识管理系统。

西夏语语法

Tue, 22 Jan 2013 00:00:00 +0000

西夏文虽然是模仿中原的汉字所创造的文字，但是由于西夏语属于汉藏语系中的藏缅语族，而汉语是汉语语族，所以两者在语法上有很大的差别。

云图脸盲症的福音

Mon, 14 Jan 2013 00:00:00 +0000

本文前半部分无剧透，可作为电影观看前预习材料。

最近看了一部神奇的电影，电影讲了一个六道轮回的故事，同一个演员饰演多角色，而且化妆师们的技术真是出神入化啊，男男女女傻傻分不清楚，总是使人在不同的时代，不同的故事中混乱游走。估计这是导演故意安排的，因为故事原本就是想要告诉我们：Everything is connected. 同时，估计这部片子会像黑客帝国那样成为 cult 片吧。不过，作为脸盲症患者，再有人鄙视我，我就跟他/她说，咱们去看云图吧。

豆瓣九点和博客访问量

Tue, 11 Dec 2012 00:00:00 +0000

到现在为止，我的博文两次上了豆瓣九点科技，所以我有两份豆瓣九点带来访客的数据。既然有了数据，那自然要利用一下。

我们的宇宙可能由闭宇宙演化为开宇宙么

Tue, 26 Jul 2011 00:00:00 +0000

我们的宇宙可能由闭宇宙演化为开宇宙么

宇宙扰动的来源以及演化

Mon, 04 Jul 2011 00:00:00 +0000

在很久很久之前，人们对于宇宙起初的那段时间的物理还不是很清楚，所以曾经有人猜测宇宙的原初扰动就是热扰动。

宇宙学到底悲剧到什么程度

Sat, 07 May 2011 00:00:00 +0000

Flowchart of Cosmology

爱情的严格量子力学描述

Fri, 18 Mar 2011 00:00:00 +0000

用量子力学的基础知识来描述爱情。

科学家才懂的加强版

Thu, 03 Mar 2011 00:00:00 +0000

严格的说，是只有数学专业和物理专业等才能读懂的~ 加了一些，重排了一下~

Volume of 20D Ball

Sat, 26 Feb 2011 00:00:00 +0000

Using the simplest program to calculate the 2D Ising model upon triangle lattice.

A Brief Summary of Statistical Physics

Fri, 21 Jan 2011 00:00:00 +0000

A Brief Summary of Statistical Physics, from Ensemble to Critical Phenomenon

2D Ising Model (Triangle Lattice)

Sat, 08 Jan 2011 00:00:00 +0000

Using the simplest program to calculate the 2D Ising model upon triangle lattice.

Monte Carlo Investigation of Frustrated Heisenberg Model

Sat, 01 Jan 2011 00:00:00 +0000

A simple vector model in 2 dimensions.

土星环的理论解释

Tue, 30 Nov 2010 00:00:00 +0000

关于土星环

关于巨正则系综巨配分函数收敛性等的一些讨论

Sat, 23 Oct 2010 00:00:00 +0000

PDF

总结

问题

$$\Xi =\sum \left(e^{-\alpha }e^{-\beta \epsilon }\right)^N$$

收敛条件 $\lvert\text{z$\phi $}\rvert<1$。重点不是要求 $\Xi$ 收敛，关键是要求热力学量收敛。（后面我们会证明对于这种经典无相互作用情况，$\Xi$ 不收敛，一定导致热力学量N，U不收敛。）

犹记多情

Wed, 22 Sep 2010 00:00:00 +0000

15号凌晨古剑奇谭游戏通关，16号凌晨在回家的火车上，突然听到有人在唱星语心愿，很有感触，就把歌词改掉了。

关于说谎的理论

Wed, 08 Sep 2010 00:00:00 +0000

说谎，目的是要使得别人相信你所说为真。理想的谎言分为弱理想谎言和强理想谎言。弱理想谎言就是使得一般智力的人不能揭穿你的谎言，强理想谎言是指的理论上不能跟区分你的谎言是真是假。

古剑奇谭晴雪心祭

Thu, 19 Aug 2010 00:00:00 +0000

15号凌晨古剑奇谭游戏通关，16号凌晨在回家的火车上，突然听到有人在唱星语心愿，很有感触，就把歌词改掉了。

余情幽梦

Mon, 07 Jun 2010 00:00:00 +0000

一直觉得毕业没啥，一直想不明白为什么会有人哭成个泪人～昨天在车上，突然就有了感觉～好多人都会一辈子不再相见呢，流泪也很正常啦～

科学家公布地磁场新数据，人类生存面临重大考验

Thu, 01 Apr 2010 00:00:00 +0000

现代人们的生活离不开电子产品，但是电子产品有一个很大的缺点，那就是在有较强的电磁或者高能粒子干扰的时候便会发生故障。而在我们的太阳系中，到处都弥漫着太阳风。太阳风是从太阳上发射出来的等离子体流形成的。

拆了 tasco 的双筒

Mon, 22 Mar 2010 00:00:00 +0000

大家都说不要拆双筒，可是这个双筒已经很脏了。

第一次拍星空

Sun, 21 Feb 2010 00:00:00 +0000

拿一小DC拍星空

莫失莫忘

Sat, 06 Feb 2010 00:00:00 +0000

给莫失莫忘的曲配的歌词。

今生醉 | 篡改朱砂泪/尘世缘

Fri, 29 Jan 2010 00:00:00 +0000

想那日烟花残照，菱花暮影容颜娇；值今夜一川烟草，天悬星河人已悄；青鸾一朝雾影落，哪堪望舒纤尘着；女萝已无离香草，煮酒寻来叶难烧。—— 以段慰藉两年前的那个傻瓜

一弯新日挂西天 — 记2010年1月15日日环食

Fri, 15 Jan 2010 00:00:00 +0000

我们观测到的是15：33开始，刚开始某人一直说他看不出来，但是我们都看出来了，也不管了，就记下了15：33这个开始时间，比预报稍晚。月亮从右下角开始遮挡。

2009年狮子座流星雨观测

Wed, 18 Nov 2009 00:00:00 +0000

没有录音和书面记录，所以仅凭记忆写下来。

我是不是真的记得

Tue, 16 Sep 2008 00:00:00 +0000

我总感觉我记得那些事，那些永远都不会有人记起的事情。

落叶

Fri, 02 May 2008 00:00:00 +0000

这是一场预谋。

诗歌三首

Mon, 27 Aug 2007 00:00:00 +0000

流星，让黑夜把握吞噬，这个世界没有门

遗忘的未来

Mon, 27 Aug 2007 00:00:00 +0000

我知道未来有太多变数，但是无论如何都不要放弃。你过去的梦想永远都不能被抛弃。

Mon, 01 Jan 0001 00:00:00 +0000

Exported from Notepad++

/*This program calculates the energy, heat capacity, magnetization, magnetic permitivity.*/ /*Boundary condition: periodic boundary condition.*/ /*In this program, J&k are put into T. To fill it to the expressions, just consider the dimensions.*/ #include <stdio.h> #include <stdlib.h> #include <time.h> #include <math.h> #define L 32 #define NGL 100000 /*Samples Neglected*/ #define N_t 5100000 /*Total samples*/ #define N (N_t-NGL) #define T_i 0.5 #define T_f 10.5 #define Tstep 0.05 #define SIZE_T 500 /*The size of emcx in T axis. Check this before compiling!*/ #define SIZE_R 60 /*The size of emcx in R axis. Check this before compiling!*/ #define R_i -6 /*Ration of J2/J1. Initial value.*/ #define R_step 0.1 /*step*/ #define R_f 0 /*End of R*/ int sum_M; int s[L][L]; /*Size of the lattice*/ double sum_E; double T; struct Obsv { double C; double X; double M; double E; }; void Inist(); void Mag(); void Energy(double R); void MCP(double R); struct Obsv Cac_obsv(double R); int main() { int i,j; double R; struct Obsv emcx[SIZE_T][SIZE_R]; FILE *ft,*fm,*fx,*fc,*fe,*fr; fr=fopen("R.txt","w"); ft=fopen("T.txt","w"); fx=fopen("X.txt","w"); fc=fopen("C.txt","w"); fm=fopen("M.txt","w"); fe=fopen("E.txt","w"); R=R_i; srand(time(NULL)); for(j=0;R<=R_f;j++) { T=T_i; for(i=0;T<=T_f;i++) { Inist(); emcx[i][j]=Cac_obsv(R); fprintf(fr,"%f\n",R); fprintf(ft,"%f\n",T); fprintf(fx,"%f\n",emcx[i][j].X); fprintf(fc,"%f\n",emcx[i][j].C); fprintf(fm,"%f\n",emcx[i][j].M); fprintf(fe,"%f\n",emcx[i][j].E); printf("-------------\n"); printf("R=%f\n",R); printf("T=%f\n",T); printf("X=%f\n",emcx[i][j].X); printf("C=%f\n",emcx[i][j].C); printf("M=%f\n",emcx[i][j].M); printf("E=%f\n",emcx[i][j].E); T = T+Tstep; } R+=R_step; } fclose(ft); fclose(fx); fclose(fc); fclose(fm); fclose(fe); printf("The End!\n"); } void Inist() { int i,j; for(i=0;i<L;i++) { for(j=0;j<L;j++) { s[i][j]=1; } } } void MCP(double R) { int i,j; double sum,dE; i=rand()%L; j=rand()%L; sum=s[i][j]*(s[(L+i-1)%L][j]+s[(i+1)%L][j]+s[i][(j-1+L)%L]+s[i][ (j+1)%L]+R*(s[(i+1)%L][(j-1+L)%L]+s[(i-1+L)%L][(j+1)%L])); //printf("chk2.sum=%d\n",sum); //switch(sum) //{ // case 4: dE= 8; break; // case 2: dE= 4; break; // case 0: dE= 0; break; // case -2: dE= -4; break; // case -4: dE= -8; break; // default: printf("sum=%d. Attention!\n",sum); break; //} dE=2*sum; if( dE>0 ) { if(((double)rand()/RAND_MAX)<=exp((-1.0)*dE/T)) { s[i][j]=-s[i][j]; sum_M+=2*s[i][j]; sum_E+=dE; } } else { s[i][j]=-s[i][j]; sum_M+=2*s[i][j]; sum_E+=dE; } } struct Obsv Cac_obsv(double R) { int i; double E_avg=0.0; double E_as=0.0; double M_avg=0.0; double M_as=0.0; struct Obsv emcx; Mag(); Energy(R); for(i=0; i<N_t; i++) { MCP(R); if(i>=NGL) { E_avg+=(double)sum_E/N; E_as+=(double)sum_E*sum_E/N; M_avg+=(double)sum_M/N; M_as+=(double)sum_M*sum_M/N; } } emcx.E=E_avg; emcx.M=M_avg; emcx.C=(E_as - E_avg*E_avg)/T/T; emcx.X=(M_as - M_avg*M_avg)/T; return emcx; } void Mag() { int i,j; int dM_tmp=0; for(i=0;i<L;i++) { for(j=0;j<L;j++) dM_tmp+=s[i][j]; } sum_M=dM_tmp; } void Energy(double R) { double dE_tmp=0.0; int i,j; for(i=0;i<L;i++) { for(j=0;j<L;j++) { dE_tmp+=0.5*(s[i][j]*(s[(i-1+L)%L][j]+s[(i+1)%L][j]+s[i][(j-1+L)%L]+s[i][(j+1)%L])+R*s[i][j]*(s[(i+1)%L][(j-1+L)%L]+s[(i-1+L)%L][(j+1)%L])); // printf("chk1.dE_tmp=%f\n",dE_tmp); } } sum_E=dE_tmp; }

About

Mon, 01 Jan 0001 00:00:00 +0000

这就是作者

一个物理出身的数据科学家，一个科幻爱好者，一个旅居海外的人。

这是与工作的简介: leima.is.

一些链接

Slides:

Colonization - Kerbal Space Program

Mon, 01 Jan 0001 00:00:00 +0000

Valentina and Jebediah initiated the KCLRP, which is short for Kerbal Colonization Long Range Plan.

The goal is to build base on every planet and moon in the whole Kerbol system.

The idea came to Valentina when she landed at the south pole some years ago.

Valentina Thinking about The Future

Project Nomand

Project Nomand is the first project in aim of building a functional self-sustained extra-Kerbin base.

Darkside

Mon, 01 Jan 0001 00:00:00 +0000

段子

王的男人之私密服务。#论如何翻译 Kingsman#（2015-04-06）#谢ErbB4#
意大利通心粉真是好东西，心里堵得慌的时候就吃点，效果倍儿棒。（2015-06-29）#谢ErbB4#
学妈：摸着你的良心说话。学爸：我的心是热的，只有热心，没有凉心。（2015-06-30）#谢ErbB4#
故事，平方反比衰减。（2015-07-03）
在时差中旅行，景色死去，唯有时差。（2015-07-06）
有个大学专门培养贝勒的，很多毕业生都成为了贝勒爷。说的就是 Baylor University.
出发，未知射线；抵达，未知射线。#机场 X 射线扫描#
双相障碍就是在液相和气相两相之间调节有障碍。
乌鸦往瓶子里投石头的时候把瓶子打碎了。
@ErbB4 的名言：你好闲，少放点盐。#女朋友的名言#
@ErbB4:science is drug. Then I interpreted it as “science is porn”.
“我眼镜里进了沙子。” “shards? 具体是什么? 玻璃、小石子还是铁丝啊？” （2015-07-16）
其实杜甫是很关心女性生活的，比如他的诗里面就有：自经丧乱少睡眠，长夜沾湿何由彻。
我猜拳很厉害的，因为我每天早上都会对着镜子练习的，这么多年以来，镜子里的那个人从来没有打败过我，然而我也没有打败过他。如果没有镜子里的那个高手，我想，我这一生，好孤独。
大风起兮裙飞扬。（从果壳老猫那里听说的）
你有没有听说一首歌叫做相约酒吧啊。
In the future, a soul encounters another intelligent AI on the internet, one of them said, “Let’s carbon date”.
在遥远的未来，智慧的人类发明了曲率引擎。作为商业运作的第一艘飞船，每当人们登船之后，广播中会反复提示：半人马座 Alpha 将会在本船停靠。
学生：那个方程不是应该叫做欧拉—拉格朗日方程么？为什么只说是欧拉方程？老师：是啊，所以我才简称欧—拉方程啊。
@ErbB4 的段子：“今天晚饭吃的什么啊啊？” “晚饭吃的食堂。” “好大的胃口”
“我是中国银行的卡。” “中国哪家银行？” / “我是美国银行的卡。” “美国哪家银行？”
其实当年孙悟空被压山下的时候是有“充电器 + WIFI”的。
请问您就读于？哦加州大学水浒分校。什么学校？UC Riverside.
真是满嘴跑火车：污污污污污。
我佛瓷杯。来源
In vivo we trust. A Nature Article
有些人是学术长胖型，有些人是学术减肥型。
喜欢写程序的人，其实是一直在使用 codeine, 毕竟是用 code 来止痛的嘛。
‘I need gas for my car.’ ‘Alright, two dimensional electron gas or neutrino gas?’

科技

突然发现一个挺有意思的巧合：火箭方程的解里面如果重新写成火箭载重/火箭燃料的形式，结果正好是个 Bose-Einstein 分布，Delta V 相当于能量，v_e 相当于 kT，燃料的 exhaust velocity 相当于能量。反正这个比值可以大于1，B-E 分布粒子数也可以大于1.（2015-06-30）

生活

SFO 这个机场如果需要过夜，要点如下：
- 二楼人少，但是在 UA check in 的附近有两张大桌子，带有电源，如果没有抢到大桌子，二楼还有其它的靠墙的很安静的椅子，可以躺下来睡觉。
- 三楼有淋浴可以洗澡，另外吃饭的地方都在三楼，必不可少的星巴克也有。除此之外还有一个很多明星的物品的博物馆，很小很小，但是东西倒是值得看看。博物馆旁边有几家餐馆，太晚了就关了，不过有些桌子可以坐下来。
- 听说在 G91 入口附近有可以付费暂存行李的地方，因为我也没存所以并不知道具体方位，应该在三楼。

Jebediah Kerman - Kerbal Space Program

Mon, 01 Jan 0001 00:00:00 +0000

Jebediah is the very first hero of Kerbin. He, in fact, guided the little green men to space.

Kerbal Space Program

Mon, 01 Jan 0001 00:00:00 +0000

Physics of Kerbal Space Program

Kerbal space program is a great game.

This is a CC BY-SA licensed project. Use the source! Keep the source open!

Summary

Rocket and Orbit Basics
- English
- 中文

KSP Log

KSP 相关

关于 KSP 这个游戏：

官方网站：https://kerbalspaceprogram.com/
Steam 页面：http://store.steampowered.com/app/220200
豆瓣游戏：http://www.douban.com/subject/24756099/
一些工具网站：
- MOD - Kerbalstuff: http://kerbalstuff.com
- MOD - Curse：http://www.curse.com/ksp-mods/kerbal
- Craft Sharing - KerbalX: http://kerbalx.com
- Kerbal Proof: http://kerbalproof.com/
- Maps - Kerbal Maps: http://www.kerbalmaps.com/
- Tutorial - Fly on budget: http://flyonbudget.one-giant-leap.info/index.html

KSP 中的物理

这是一系列文章，简要汇总 Kerbal space program 中的物理。

Research Tools

Mon, 01 Jan 0001 00:00:00 +0000

{% include zero-md/zero-md-scripts.html %}

I started a collection of research tools on GitHub, which is currently a part of the awesome-awesome collection.

You can

Star on GitHub

Comment

The following content is live rendered using zero-md. If something went wrong, please read it on GitHub repo: emptymalei/awesome-research.

Top 10 Contributors: Contributors@GitHub Graph; Pull Requests

Typograph Templates

Mon, 01 Jan 0001 00:00:00 +0000

Figures:

{{< figure align=center src="../assets/" caption="" >}}

Math environment without escapes

{{< m >}}
\alpha
{{< /m >}}

火箭和轨道

Mon, 01 Jan 0001 00:00:00 +0000

引力

在 KSP 中所有用到的引力规律都是两体作用，在目前的版本中 1.0 中并没有考虑三体和以上的作用。这种情况下，一个质量为 $m$ 的物体所收到的引力可以用以下公式来计算，

聊天

Mon, 01 Jan 0001 00:00:00 +0000

我做了几年物理，后来做了几年数据科学。目前在德国生活，业余出去转转，也偶尔写点自己的小项目。

职业发展

如果想聊聊职业发展，可以在 ADPList 上找我约时间。使用如下工具查看我提供的时间：

虚拟厨艺

Mon, 01 Jan 0001 00:00:00 +0000

虚拟厨艺，就是使用电脑来作出食物的一种艺术。

所有的例子都在我的下厨房

清蒸苹果

ParametricPlot3D[{(1 + Cos[v]) Cos[u] +
  0.085 Cos[5 u] + (0.994 v/\[Pi])^100, (1 + Cos[v]) Sin[u],
 6 Sin[v] + 2 Cos[v] -
  0.7 Log[1 - v/\[Pi]]}, {u, -\[Pi], \[Pi]}, {v, -\[Pi], \[Pi]},
Mesh -> False, BoxRatios -> 1, Boxed -> False, Axes -> None,
PlotStyle -> {Green, Specularity[Yellow, 10]}, ImageSize -> 500]

结果如下:

娱乐

Mon, 01 Jan 0001 00:00:00 +0000

本分类下的文章使用 CC BY-SA licensed. Use the source! Keep the source open!

Dark Hypernova

2025 年

Team is Different

2023 年

弓箭笔记

2022 年

2021 年

做了三年 Data Scientist

2020 年

2019 年

生活

工作

做了六个月的 Data Scientist

从学术界到工业界

Beyond the Professoriate Webinar

云朵工程师

美丽新世界

关于圈子的一些自我反思

什么样的生活

幽灵盖娅

量子心理学

热力学势

火星的天空

梦与新世界

孩童的宇宙

今夜请把我忘记

面条儿

奈何桥

一封永远寄不出的信

一个「精神病人」的救赎

抑郁而死

阅读者

四元数和转动

僵尸与统计物理

题目

僵尸的悲惨故事

量子爱情，统计数据以及认知模型

你我之上的幽灵

物理中的作用量

设计者的宏伟蓝图

心理叠加态：从量子爱情到人的所有心理

爱情的量子理论2

爱情的严格量子力学描述

球状闪电：千百年的幽灵

Creative Commons License 图标设计与 SVG

空间站人造重力

量子力学框架回顾一

量子力学计算的一个小技巧

关于宇宙膨胀的两个常见误解

如何重建宇宙演化历史（未完成）

A Basic Jacobi Iteration Procedure

Finite Difference Discretization

Iterative Solution Methods for Ax=b

The Advection-diffusion Equations

Implementing MATLAB's Boundary Value Solver

如何挑西瓜：一个学物理的吃货的感言

Boundary Value Problems: Direct Solve and Relaxation

Implementation of Shooting and Convergence Studies

Boundary Value Problem

Error Analysis for Time-Stepping Routines

Euler, Runge-Kutta and Adams Methods

Initial Value Problem

线性代数与几何

一页纸上的微积分

如何攻击巨人——立体行动装置使用要点

PKM —— 个人知识管理

西夏语语法

云图 脸盲症的福音

豆瓣九点和博客访问量

我们的宇宙可能由闭宇宙演化为开宇宙么

宇宙扰动的来源以及演化

宇宙学到底悲剧到什么程度

爱情的严格量子力学描述

科学家才懂的 加强版

Volume of 20D Ball

A Brief Summary of Statistical Physics

2D Ising Model (Triangle Lattice)

Monte Carlo Investigation of Frustrated Heisenberg Model

土星环的理论解释

关于巨正则系综巨配分函数收敛性等的一些讨论

云图脸盲症的福音

科学家才懂的加强版

古剑奇谭晴雪心祭