取整截断函数及其在PBC中的使用

李继存 2013-8-17 07:47

取整截断函数及其在PBC中的使用 2013-08-13 14:05:39 取整和截断相关的函数，有着各种不同的名称，不同的定义，很容易让人迷惑。它们之间根据具体情况还可以相互转换成不同的写法，更增添了混乱。下面整理 Fortran 和 C 中的相关函数。 函数 Fortran C 向下取整 ...

个人分类: 数学轮子|13395 次阅读|没有评论

利用PyMOL制作反应势能面示意动画

热度 1 李继存 2013-4-18 23:50

2013-04-18 10:43:35 发布图片 2013-11-14 21:18:18 增加教程 2016-09-07 20:17:18 颜色映射 利用PyMOL内置的OpenGL制作了势能面反应的示意图, 具体制作过程容后再讲. PyMOL是基于Python的分子可视化软件, 在结构生物学中使用十分广泛. 它的主要特色是对生物大分子显示效果好, 并自带一个高效的光 ...

个人分类: 数学轮子|11917 次阅读|2 个评论 热度 1

空间中两随机向量间夹角的概率密度分布

李继存 2013-3-19 09:59

设空间维数为 n ，则此空间中任意两个向量之间的夹角 $\theta$ 符从一定的分布，其概率密度函数如下： $p(\theta)=\frac{\Gamma(\frac{n}{2})}{\Gamma(\frac{n-1}{2})} \frac{\sin^{n-2}(\theta)}{\sqrt\pi} $ 对二维平面， $n=2, p(\theta)=\frac{1}{\pi}$ 对三维空 ...

个人分类: 数学轮子|9764 次阅读|没有评论

一元三次方程求根公式及其Fortran代码

热度 2 李继存 2012-10-31 00:10

一元三次方程的求解比一元二次方程困难，求根公式看起来也很复杂， 中文 维基百科 上原先的三角函数解公式有误，我已经修改过了。但是我觉得还有必要将这些公式简化一下，给出具体的算法，以便编程使用。 首先需要明确，一元三次方程 $f(x)=ax^3+bx^2+cx+d=0,(a\neq0) $ 至少有一个实根 ...

个人分类: 数学轮子|24378 次阅读|5 个评论 热度 2

二维三次卷积插值算法及Fortran代码【更新】

李继存 2012-7-24 13:26

最近有人问起二维三次卷积插值算法（ Cubic Convolution Interpolation ）及其程序的问题，我想到这种插值算法也能用于图像的插值，所以也就有了兴趣，稍微了解了一下，并试着用 Fortran 实现了一下，同时也顺便复习了一下 Fortran90 的一下新特性。如果需要，敬请使用；发现 Bug ，烦请明示。 ...

个人分类: 数学轮子|10400 次阅读|没有评论

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