| 查看: 664 | 回复: 1 | |||
| 当前主题已经存档。 | |||
ganchunlei至尊木虫 (著名写手)
|
[交流]
元胞自动机模拟
|
||
|
元胞自动机模拟 cellular automaton /*------------------------------------------------------------------------------* * File Name: Cellular Automator example for Origin C * *------------------------------------------------------------------------------*/ #include //********************************** // to implement a command // ca // for a simple celluar automator program // usage: // 1. compile this file // 2. close CodeBuilder for faster execution // 3. create a new matrix // 4. from script window, type ca and press enter //********************************** //////////////////////////////////////////////////////////////////////////////////// //static int temp[8] = {0,0,0,1,1,1,1,0}; // Rule 30 static int my_rules[8] = {0,0,0,1,1,0,1,0}; //////////////////////////////////////////////////////////////////////////////////// // our basic rule function // we take all the 8 possible permutations for the three cells above the // current cell and we need to determine live (1) or die (0) // and the determinations are stored in the my_rules array above // // use static function to avoid being visible from LabTalk static int rule(int left, int top, int right) { if(left && top && right) return my_rules[0]; if(left && top && !right) return my_rules[1]; if(left && !top && right) return my_rules[2]; if(left && !top && !right) return my_rules[3]; if(!left && top && right) return my_rules[4]; if(!left && top && !right) return my_rules[5]; if(!left && !top && right) return my_rules[6]; // for left=top=right=0 return my_rules[7]; } // this is our main function to be called from labtalk void ca() { MatrixLayer ml = Project.ActiveLayer(); if(ml==NULL) { out_str("Must have matrix as active window"  ;return; } ml.SetInternalData(FSI_BYTE); ml.SetNumCols(1000); ml.SetNumRows(500); Matrix int nCols = mat.GetNumCols(); //--------- // init 1st row of matrix to be all zero except 1 in the middle for(int ii = 0; ii < nCols; ii++) { mat[0][ii] = 0; } mat[0][nCols/2] = 1; //--------- int left,top,right; // start from 2nd row progressBox show("calculating...." ;show.SetRange(1,mat.GetNumRows()); for(int nR = 1, nRLast=0; nR < mat.GetNumRows(); nR++,nRLast ++) { if(!show.Set(nR)) // user click cancel break; for(int nC = 0; nC < nCols; nC++) { left = nC > 0? mat[nRLast][nC-1] : 0; top = mat[nRLast][nC]; right = nC < nCols-1? mat[nRLast][nC+1] : 0; mat[nR][nC] = rule(left,top,right); } } ml.SetViewImage(TRUE); ml.LT_execute("Z1=0;Z2=1" ;// this is needed to set Z range to show bitmap as black and white} //----------------------------- EOF |
回复此楼» 猜你喜欢
什么时候开奖？
已经有4人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
2楼2006-10-16 17:47:00
