Astronomical Software Club

看下面这段代码：

#include <stdio.h>
#include <stdlib.h>
#include <omp.h>
 
void para_print()
{
	//automatic vars in a routine called within an omp region is private implicitly
	int i,tid;
	/*if you need shared vars inside omp function call, use static vars;
	* the side effect is static vars keeps its previous value 
        * if not reset a initial value during each call*/
	static int j;
	#pragma omp single
	{
	printf("tid: i,j\n");
	j=10;
	}
	#pragma omp for 
	for(i=0;i<5;i++)
	{
	tid=omp_get_thread_num();
	#pragma omp critical
	{
		j++;

多个并行域之间，或者并行域内部的函数调用之间如果需要在线程内部共享全局变量，必须用threadprivate。如果不将全局变量声明为threadprivate而是声明为private，则在子函数调用时子函数内的全局变量传递不到子函数外面。

比如下面程序：

Matlab 函数进阶：

使用匿名函数(Anonymous Function)和内嵌函数(Nested Function)处理多变量传递问题(Matlab 7.0以上)

问题：有一个多变量函数$f(a,b,c,x)$,现需要分别在$a=a0,b=b0,c=c0$和$a=a1,b=b1,c=c1$的条件下对$f(a,b,c,x)$进行某一操作。

此类问题常在数值积分时出现。

解决方案：

1． 使用全局变量

 打包下载(.zip)：（包括文档，源代码及编译后的程序，以及模拟动画）

程序说明： 

EMRUI V1.0b（2004圣诞预览版）
 ——电磁流变液体固化的二维数值模拟图形界面

一：操作
设定初值后，点“计算”按钮，则开始计算，同时绘图区同步显示运算过程。请耐心等待，待弹出提示"Data ready"后按“放映”即可观看动画。

从这里下载该文件

预览：

matlabcentral为全球matlab用户共享资源的地方，在那里你可以找到各个学科的很多你需要的东西，当然你也可以把自己的程序发布到那里同全球的用户共享交流：http://www.mathworks.com/matlabcentral/

关于天文方面的matlab资源，自己在网上搜罗了下面这些东西，欢迎大家补充

General Websites：