मैं C ++ कोड के साथ ट्यूरिंग की प्रतिक्रिया-प्रसार प्रणाली को हल कर रहा हूं। यह बहुत धीमा है: 128x128 पिक्सेल बनावट के लिए, पुनरावृत्तियों की स्वीकार्य संख्या 200 है - जिसके परिणामस्वरूप 2.5 सेकंड की देरी होती है। मुझे दिलचस्प छवि प्राप्त करने के लिए 400 पुनरावृत्तियों की आवश्यकता है - लेकिन प्रतीक्षा के 5 सेकंड बहुत अधिक हैं। इसके अलावा, बनावट का आकार वास्तव में 512x512 होना चाहिए - लेकिन इसके परिणामस्वरूप बहुत इंतजार करना पड़ता है। उपकरण iPad, iPod हैं।
क्या इसे तेजी से करने का कोई मौका है? यूलर विधि धीरे-धीरे (विकिपीडिया) में परिवर्तित हो जाती है - जल्दी विधि होने से पुनरावृत्तियों की संख्या को छोड़ने की अनुमति होगी?
EDIT: जैसा कि थॉमस क्लिंपेल ने बताया, "अगर (m_An [i] [j] <0.0) {...}", "if (m_Bn [i] [j] <0.0) {...}" अभिसरण में देरी कर रहे हैं: हटाने के बाद, 75 पुनरावृत्तियों के बाद सार्थक छवि दिखाई देती है । मैंने नीचे कोड में लाइनें लिखी हैं।
void TuringSystem::solve( int iterations, double CA, double CB ) {
m_iterations = iterations;
m_CA = CA;
m_CB = CB;
solveProcess();
}
void set_torus( int & x_plus1, int & x_minus1, int x, int size ) {
// Wrap "edges"
x_plus1 = x+1;
x_minus1 = x-1;
if( x == size - 1 ) { x_plus1 = 0; }
if( x == 0 ) { x_minus1 = size - 1; }
}
void TuringSystem::solveProcess() {
int n, i, j, i_add1, i_sub1, j_add1, j_sub1;
double DiA, ReA, DiB, ReB;
// uses Euler's method to solve the diff eqns
for( n=0; n < m_iterations; ++n ) {
for( i=0; i < m_height; ++i ) {
set_torus(i_add1, i_sub1, i, m_height);
for( j=0; j < m_width; ++j ) {
set_torus(j_add1, j_sub1, j, m_width);
// Component A
DiA = m_CA * ( m_Ao[i_add1][j] - 2.0 * m_Ao[i][j] + m_Ao[i_sub1][j] + m_Ao[i][j_add1] - 2.0 * m_Ao[i][j] + m_Ao[i][j_sub1] );
ReA = m_Ao[i][j] * m_Bo[i][j] - m_Ao[i][j] - 12.0;
m_An[i][j] = m_Ao[i][j] + 0.01 * (ReA + DiA);
// if( m_An[i][j] < 0.0 ) { m_An[i][j] = 0.0; }
// Component B
DiB = m_CB * ( m_Bo[i_add1][j] - 2.0 * m_Bo[i][j] + m_Bo[i_sub1][j] + m_Bo[i][j_add1] - 2.0 * m_Bo[i][j] + m_Bo[i][j_sub1] );
ReB = 16.0 - m_Ao[i][j] * m_Bo[i][j];
m_Bn[i][j] = m_Bo[i][j] + 0.01 * (ReB + DiB);
// if( m_Bn[i][j] < 0.0 ) { m_Bn[i][j]=0.0; }
}
}
// Swap Ao for An, Bo for Bn
swapBuffers();
}
}