数字旗手

本文参考了百度Mathematica吧的这个帖子。 统计图片中某个区域的面积可以有多种方法和多种软件，比如ImageJ(多谢柏爷推荐~)、ImageMagick(参见这篇教程)、PhotoShop(参见这篇教程)。这里采用Mathematica，仅用少于10行代码即可实现。 Mathematica源码为： 1 2 3 4 img2 = Binarize[img]; totalPixels = Times @@ ImageDimensions[img2]; blackPixels = Count[ImageData[img2], a_ /; a < 0.5, {2}]; blackAre

析出反应 通常来说，过饱和亚稳合金中的析出反应可分为三个不同的阶段： 1. 析出相的形核，这来源于过饱和基体中局部区域的成分起伏，且符合统计规律。 2. 析出相的生长，其受溶质元素扩散控制，且与周围的第二相不产生竞争关系。 3. 析出相的粗化，此过程中小尺寸的颗粒溶解，使得大尺寸的颗粒继续生长。 前两个过程是由相变的自由能差驱动，其在达到平衡相体积分数后就结束，而粗化则只是界面能降低的过程，可以一直持续到仅有一个颗粒为止。另外，生长和粗化这两个阶段较难区分，实际上这里是为了分析方便，对于不同阶段配合真实的边界条件而采用不同的处理方法。 形核 形核模型可采用基于随机性形核的经典形核理

写在前面——如何阅读OpenPhase源码 OpenPhase的源码真的是挺难读的，它是一整套完全的很成熟的求解相场模型的框架，其模块化完成度很高(尤其设计的数据结构面面俱到，号称能处理所有类型的数据)，所以用它来解方程时写的代码很少，但带来一个问题是你必须先读懂它：高度集成的代码让初学者“一进去就出不来”(从而迷失在各个类中)或“压根进不去”(一行代码就完成一个大功能，根本不知道怎么入手)。 我的建议是： 1. 善用Doxygen OpenPhase支持Doxygen，其注释也较全面，使用Doxygen生成类之间的关系图，方便理清思路。 附录有Doxygen生成的dot图

持续更新中~~~ 条件断点 为断点设置条件表达式，只有条件满足时才激活断点 如： break 670 if sum == 3 转自这里 出现”value optimized out”错误 gdb调试程序的时候打印变量值会出现value optimized out情况,可以在gcc编译的时候加上-O0参数项,意思是不进行编译优化,调试的时候就会顺畅了,运行流程不会跳来跳去的,发布项目的时候记得不要在使用-O0参数项,gcc默认编译或加上-O2优化编译会提高程序运行速度. 转自这里。

Update 2016-11-16更新： 在Ubuntu10.10上安装gcc4.8.0时，出现错误： 1 ‘CHAR_BIT’ was not declared in this scope 解决方法： 1 unset LIBRARY_PATH CPATH C_INCLUDE_PATH PKG_CONFIG_PATH CPLUS_INCLUDE_PATH INCLUDE 如果在/etc/profile 或者 ~/.bashrc中设置了ccache的话,也暂时去掉ccache的设置. 参考自这里。 另外： 单纯修改gcc的指向后，会造成其他编译器的不匹配，比如gcc虽然指

写在前面：本文转自原AFEPack论坛afepack.org。 AFEPack教学系列转载自北京大学数学学院李若教授之个人博客。 如果需要转载本教学系列，请保持文章完整性，并注明转载出处。 李老师的话 学校的服务器老是有问题，希望这个地方能够稳定一些。我把AFEPack的使用写一些系列的小例子，从最简单的开始，希望能够为我这个软件包的使用者提供一些方便。 我在大年三十的晚上更新了AFEPack 的主页上的一点内容，现在主要是照抄下来，给想用AFEPack的人一些方向性的指导。在过去的几年中间，我多少花了一些力气来推广这个包，希望能够为计算数学的发展做一点事情。由于一些令人失望的经验，我对

模型概况 该模型是Appolaire等人于2008年在MSE上发表。 B. Appolaire, H. Combeau, G. Lesoult. Modeling of equiaxed growth in multicomponent alloys accounting for convection and for the globular/dendritic morphological transition, Mater. Sci. Eng. A 487 (2008) 33-45. 该模型是一个半解析半数值的模型，一方面对于枝晶的细微结构如一次枝晶尖端半径、二次枝晶间距及直径，对流效应

安装AFEPack软件包，需要提前安好的软件有gcc、g++、doxygen、automake(1.11版本)、dx、emacs、vim。 准备： 我会有一个用户名，比如为qixinbo。在home/qixinbo下建立files，include，lib三个目录。其中include目录是用于存放程序编译时调用的.h头文件，lib目录是用于存放动态链接库.so文件 然后将 boost_1_37_0.tar.bz2, deal.nodoc-6.3.1.tar.gz, AFEPack-snapshot.tar.gz, easymesh.c.gz 四个文件放在files里。 安装boost： 安

纯铁在不同条件下的平衡态计算 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 r t steel1 //read TDB steel的缩写，读取steel1.TDB文件 fe //选择Fe元素 l data //list data的缩写，列出所有的热力学数据 set c t=1000 p=1e5 n=1 //set condition的缩写，设置温度K，压力Pa和摩尔数mol c e //calculate equilibrium的缩写，计算当前条件下的系统平衡态 l,,,, //list的缩写

本文是对徐祖耀老师所著的《相变原理》(科学出版社2000年出版)第一章的读书笔记。 相变的定义 在均匀单相中，或在几个混合相中，出现了不同成分或不同结构(包括原子、离子或电子位置及位向的改变)、不同组织形态或不同性质的相，就称为相变。 相变何时发生 当母相失稳而新相具有较高的稳定性时，就会发生相变。 母相失稳可以需要跨过一定的能垒(需要一定的形核驱动力，即先形核)，也可以连续失稳(没有能垒的连续型相变，不需要进行形核过程，如Spinodal分解)。 在一些固态相变中，除温度因素外，压强和应变也引起母相的失稳。 固态相变热循环时往往发生较大的热滞，即冷却时的相变温度不等于再次加热时的相变温度