Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed
Identifieur interne : 000004 ( Pmc/Corpus ); précédent : 000003; suivant : 000005Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed
Auteurs : Eric J. Nichols ; Axel HuttSource :
- Frontiers in Neuroinformatics [ 1662-5196 ] ; 2015.
Abstract
Neural Field models (NFM) play an important role in the understanding of neural population dynamics on a mesoscopic spatial and temporal scale. Their numerical simulation is an essential element in the analysis of their spatio-temporal dynamics. The simulation tool described in this work considers scalar spatially homogeneous neural fields taking into account a finite axonal transmission speed and synaptic temporal derivatives of first and second order. A text-based interface offers complete control of field parameters and several approaches are used to accelerate simulations. A graphical output utilizes video hardware acceleration to display running output with reduced computational hindrance compared to simulators that are exclusively software-based. Diverse applications of the tool demonstrate breather oscillations, static and dynamic Turing patterns and activity spreading with finite propagation speed. The simulator is open source to allow tailoring of code and this is presented with an extension use case.
Url:
DOI: 10.3389/fninf.2015.00025
PubMed: 26539105
PubMed Central: 4611063
Links to Exploration step
PMC:4611063Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed</title>
<author><name sortKey="Nichols, Eric J" sort="Nichols, Eric J" uniqKey="Nichols E" first="Eric J." last="Nichols">Eric J. Nichols</name>
</author>
<author><name sortKey="Hutt, Axel" sort="Hutt, Axel" uniqKey="Hutt A" first="Axel" last="Hutt">Axel Hutt</name>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">26539105</idno>
<idno type="pmc">4611063</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4611063</idno>
<idno type="RBID">PMC:4611063</idno>
<idno type="doi">10.3389/fninf.2015.00025</idno>
<date when="2015">2015</date>
<idno type="wicri:Area/Pmc/Corpus">000004</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000004</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed</title>
<author><name sortKey="Nichols, Eric J" sort="Nichols, Eric J" uniqKey="Nichols E" first="Eric J." last="Nichols">Eric J. Nichols</name>
</author>
<author><name sortKey="Hutt, Axel" sort="Hutt, Axel" uniqKey="Hutt A" first="Axel" last="Hutt">Axel Hutt</name>
</author>
</analytic>
<series><title level="j">Frontiers in Neuroinformatics</title>
<idno type="eISSN">1662-5196</idno>
<imprint><date when="2015">2015</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Neural Field models (NFM) play an important role in the understanding of neural population dynamics on a mesoscopic spatial and temporal scale. Their numerical simulation is an essential element in the analysis of their spatio-temporal dynamics. The simulation tool described in this work considers scalar spatially homogeneous neural fields taking into account a finite axonal transmission speed and synaptic temporal derivatives of first and second order. A text-based interface offers complete control of field parameters and several approaches are used to accelerate simulations. A graphical output utilizes video hardware acceleration to display running output with reduced computational hindrance compared to simulators that are exclusively software-based. Diverse applications of the tool demonstrate breather oscillations, static and dynamic Turing patterns and activity spreading with finite propagation speed. The simulator is open source to allow tailoring of code and this is presented with an extension use case.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Atay, F M" uniqKey="Atay F">F. M. Atay</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Beim Graben, P" uniqKey="Beim Graben P">P. beim Graben</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brackley, C A" uniqKey="Brackley C">C. A. Brackley</name>
</author>
<author><name sortKey="Turner, M" uniqKey="Turner M">M. Turner</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bressloff, P" uniqKey="Bressloff P">P. Bressloff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bressloff, P C" uniqKey="Bressloff P">P. C. Bressloff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buckwar, E" uniqKey="Buckwar E">E. Buckwar</name>
</author>
<author><name sortKey="Kuske, R" uniqKey="Kuske R">R. Kuske</name>
</author>
<author><name sortKey="Mohammed, S" uniqKey="Mohammed S">S. Mohammed</name>
</author>
<author><name sortKey="Shardlow, T" uniqKey="Shardlow T">T. Shardlow</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buckwar, E" uniqKey="Buckwar E">E. Buckwar</name>
</author>
<author><name sortKey="Winkler, R" uniqKey="Winkler R">R. Winkler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Buckwar, E" uniqKey="Buckwar E">E. Buckwar</name>
</author>
<author><name sortKey="Winkler, R" uniqKey="Winkler R">R. Winkler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Carletti, M" uniqKey="Carletti M">M. Carletti</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Carnevale, N" uniqKey="Carnevale N">N. Carnevale</name>
</author>
<author><name sortKey="Hines, M" uniqKey="Hines M">M. Hines</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S. Coombes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S. Coombes</name>
</author>
<author><name sortKey="Beim Graben, P" uniqKey="Beim Graben P">P. beim Graben</name>
</author>
<author><name sortKey="Potthast, R" uniqKey="Potthast R">R. Potthast</name>
</author>
<author><name sortKey="Wright, J" uniqKey="Wright J">J. Wright</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S. Coombes</name>
</author>
<author><name sortKey="Laing, C" uniqKey="Laing C">C. Laing</name>
</author>
<author><name sortKey="Schmidt, H" uniqKey="Schmidt H">H. Schmidt</name>
</author>
<author><name sortKey="Svanstedt, N" uniqKey="Svanstedt N">N. Svanstedt</name>
</author>
<author><name sortKey="Wyller, J" uniqKey="Wyller J">J. Wyller</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S. Coombes</name>
</author>
<author><name sortKey="Owen, M" uniqKey="Owen M">M. Owen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coombes, S" uniqKey="Coombes S">S. Coombes</name>
</author>
<author><name sortKey="Venkov, N" uniqKey="Venkov N">N. Venkov</name>
</author>
<author><name sortKey="Shiau, L" uniqKey="Shiau L">L. Shiau</name>
</author>
<author><name sortKey="Bojak, I" uniqKey="Bojak I">I. Bojak</name>
</author>
<author><name sortKey="Liley, D" uniqKey="Liley D">D. Liley</name>
</author>
<author><name sortKey="Laing, C" uniqKey="Laing C">C. Laing</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Deco, G" uniqKey="Deco G">G. Deco</name>
</author>
<author><name sortKey="Jirsa, V K" uniqKey="Jirsa V">V. K. Jirsa</name>
</author>
<author><name sortKey="Robinson, P" uniqKey="Robinson P">P. Robinson</name>
</author>
<author><name sortKey="Breakspear, M" uniqKey="Breakspear M">M. Breakspear</name>
</author>
<author><name sortKey="Friston, K" uniqKey="Friston K">K. Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Demeulemeester, H" uniqKey="Demeulemeester H">H. Demeulemeester</name>
</author>
<author><name sortKey="Vandesande, F" uniqKey="Vandesande F">F. Vandesande</name>
</author>
<author><name sortKey="Orban, G" uniqKey="Orban G">G. Orban</name>
</author>
<author><name sortKey="Brandon, C" uniqKey="Brandon C">C. Brandon</name>
</author>
<author><name sortKey="Vanderhaeghen, J J" uniqKey="Vanderhaeghen J">J. J. Vanderhaeghen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Elvin, A J" uniqKey="Elvin A">A. J. Elvin</name>
</author>
<author><name sortKey="Laing, C R" uniqKey="Laing C">C. R. Laing</name>
</author>
<author><name sortKey="Roberts, M G" uniqKey="Roberts M">M. G. Roberts</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Faye, G" uniqKey="Faye G">G. Faye</name>
</author>
<author><name sortKey="Faugeras, O" uniqKey="Faugeras O">O. Faugeras</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Folias, S E" uniqKey="Folias S">S. E. Folias</name>
</author>
<author><name sortKey="Bressloff, P C" uniqKey="Bressloff P">P. C. Bressloff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Friston, K J" uniqKey="Friston K">K. J. Friston</name>
</author>
<author><name sortKey="Kahan, J" uniqKey="Kahan J">J. Kahan</name>
</author>
<author><name sortKey="Biswal, B" uniqKey="Biswal B">B. Biswal</name>
</author>
<author><name sortKey="Razi, A" uniqKey="Razi A">A. Razi</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hashemi, M" uniqKey="Hashemi M">M. Hashemi</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Sleigh, J" uniqKey="Sleigh J">J. Sleigh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Atay, F" uniqKey="Atay F">F. Atay</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Bestehorn, M" uniqKey="Bestehorn M">M. Bestehorn</name>
</author>
<author><name sortKey="Wennekers, T" uniqKey="Wennekers T">T. Wennekers</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Buhry, L" uniqKey="Buhry L">L. Buhry</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Lefebvre, J" uniqKey="Lefebvre J">J. Lefebvre</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Rougier, N" uniqKey="Rougier N">N. Rougier</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Rougier, N" uniqKey="Rougier N">N. Rougier</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Sutherland, C" uniqKey="Sutherland C">C. Sutherland</name>
</author>
<author><name sortKey="Longtin, A" uniqKey="Longtin A">A. Longtin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Idiart, M A P" uniqKey="Idiart M">M. A. P. Idiart</name>
</author>
<author><name sortKey="Abbott, L F" uniqKey="Abbott L">L. F. Abbott</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jirsa, V K" uniqKey="Jirsa V">V. K. Jirsa</name>
</author>
<author><name sortKey="Haken, H" uniqKey="Haken H">H. Haken</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jirsa, V" uniqKey="Jirsa V">V. Jirsa</name>
</author>
<author><name sortKey="Jantzen, K" uniqKey="Jantzen K">K. Jantzen</name>
</author>
<author><name sortKey="Fuchs, A" uniqKey="Fuchs A">A. Fuchs</name>
</author>
<author><name sortKey="Kelso, J" uniqKey="Kelso J">J. Kelso</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Kilpatrick, Z P" uniqKey="Kilpatrick Z">Z. P. Kilpatrick</name>
</author>
<author><name sortKey="Bressloff, P C" uniqKey="Bressloff P">P. C. Bressloff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Laing, C" uniqKey="Laing C">C. Laing</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Langtangen, H P" uniqKey="Langtangen H">H. P. Langtangen</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Molaee Ardekani, B" uniqKey="Molaee Ardekani B">B. Molaee-Ardekani</name>
</author>
<author><name sortKey="Senhadji, L" uniqKey="Senhadji L">L. Senhadji</name>
</author>
<author><name sortKey="Shamsollahi, M B" uniqKey="Shamsollahi M">M. B. Shamsollahi</name>
</author>
<author><name sortKey="Vosoughi Vahdat, B" uniqKey="Vosoughi Vahdat B">B. Vosoughi-Vahdat</name>
</author>
<author><name sortKey="Wodey, E" uniqKey="Wodey E">E. Wodey</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nunez, P" uniqKey="Nunez P">P. Nunez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nunez, P" uniqKey="Nunez P">P. Nunez</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nunez, P" uniqKey="Nunez P">P. Nunez</name>
</author>
<author><name sortKey="Srinivasan, R" uniqKey="Srinivasan R">R. Srinivasan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Owen, M R" uniqKey="Owen M">M. R. Owen</name>
</author>
<author><name sortKey="Laing, C R" uniqKey="Laing C">C. R. Laing</name>
</author>
<author><name sortKey="Coombes, S" uniqKey="Coombes S">S. Coombes</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pinotsis, D" uniqKey="Pinotsis D">D. Pinotsis</name>
</author>
<author><name sortKey="Friston, K" uniqKey="Friston K">K. Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pinotsis, D A" uniqKey="Pinotsis D">D. A. Pinotsis</name>
</author>
<author><name sortKey="Moran, R J Friston K J" uniqKey="Moran R">R. J, Friston, K. J. Moran</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pinotsis, D A" uniqKey="Pinotsis D">D. A. Pinotsis</name>
</author>
<author><name sortKey="Schwarzkopf, D S" uniqKey="Schwarzkopf D">D. S. Schwarzkopf</name>
</author>
<author><name sortKey="Litvak, V" uniqKey="Litvak V">V. Litvak</name>
</author>
<author><name sortKey="Rees, G" uniqKey="Rees G">G. Rees</name>
</author>
<author><name sortKey="Barnes, G" uniqKey="Barnes G">G. Barnes</name>
</author>
<author><name sortKey="Friston, K J" uniqKey="Friston K">K. J. Friston</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pinto, D" uniqKey="Pinto D">D. Pinto</name>
</author>
<author><name sortKey="Ermentrout, G" uniqKey="Ermentrout G">G. Ermentrout</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Qubbaj, M R" uniqKey="Qubbaj M">M. R. Qubbaj</name>
</author>
<author><name sortKey="Jirsa, V K" uniqKey="Jirsa V">V. K. Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rosenkilde, C" uniqKey="Rosenkilde C">C. Rosenkilde</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rossant, C" uniqKey="Rossant C">C. Rossant</name>
</author>
<author><name sortKey="Harris, K D" uniqKey="Harris K">K. D. Harris</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rougier, N P" uniqKey="Rougier N">N. P. Rougier</name>
</author>
<author><name sortKey="Fix, J" uniqKey="Fix J">J. Fix</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rougier, N" uniqKey="Rougier N">N. Rougier</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sanz Leon, P" uniqKey="Sanz Leon P">P. Sanz Leon</name>
</author>
<author><name sortKey="Knock, S A" uniqKey="Knock S">S. A. Knock</name>
</author>
<author><name sortKey="Spiegler, A" uniqKey="Spiegler A">A. Spiegler</name>
</author>
<author><name sortKey="Jirsa, V K" uniqKey="Jirsa V">V. K. Jirsa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sanz Leon, P" uniqKey="Sanz Leon P">P. Sanz Leon</name>
</author>
<author><name sortKey="Knock, S A" uniqKey="Knock S">S. A. Knock</name>
</author>
<author><name sortKey="Woodman, M M" uniqKey="Woodman M">M. M. Woodman</name>
</author>
<author><name sortKey="Domide, L" uniqKey="Domide L">L. Domide</name>
</author>
<author><name sortKey="Mersmann, J" uniqKey="Mersmann J">J. Mersmann</name>
</author>
<author><name sortKey="Mcintosh, A" uniqKey="Mcintosh A">A. McIntosh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Schmidt, H" uniqKey="Schmidt H">H. Schmidt</name>
</author>
<author><name sortKey="Hutt, A" uniqKey="Hutt A">A. Hutt</name>
</author>
<author><name sortKey="Schimansky Geier, L" uniqKey="Schimansky Geier L">L. Schimansky-Geier</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Steyn Ross, M L" uniqKey="Steyn Ross M">M. L. Steyn-Ross</name>
</author>
<author><name sortKey="Steyn Ross, D A" uniqKey="Steyn Ross D">D. A. Steyn-Ross</name>
</author>
<author><name sortKey="Wilson, M T" uniqKey="Wilson M">M. T. Wilson</name>
</author>
<author><name sortKey="Sleigh, J W" uniqKey="Sleigh J">J. W. Sleigh</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stimberg, M" uniqKey="Stimberg M">M. Stimberg</name>
</author>
<author><name sortKey="Goodman, D F" uniqKey="Goodman D">D. F. Goodman</name>
</author>
<author><name sortKey="Benichoux, V" uniqKey="Benichoux V">V. Benichoux</name>
</author>
<author><name sortKey="Brette, R" uniqKey="Brette R">R. Brette</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Turing, A" uniqKey="Turing A">A. Turing</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Van Loan, C" uniqKey="Van Loan C">C. Van Loan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Veltz, R" uniqKey="Veltz R">R. Veltz</name>
</author>
<author><name sortKey="Faugeras, O" uniqKey="Faugeras O">O. Faugeras</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Veltz, R" uniqKey="Veltz R">R. Veltz</name>
</author>
<author><name sortKey="Faugeras, O" uniqKey="Faugeras O">O. Faugeras</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, X J" uniqKey="Wang X">X.-J. Wang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wang, Y" uniqKey="Wang Y">Y. Wang</name>
</author>
<author><name sortKey="Markram, H" uniqKey="Markram H">H. Markram</name>
</author>
<author><name sortKey="Goodman, P H" uniqKey="Goodman P">P. H. Goodman</name>
</author>
<author><name sortKey="Berger, T K" uniqKey="Berger T">T. K. Berger</name>
</author>
<author><name sortKey="Ma, J" uniqKey="Ma J">J. Ma</name>
</author>
<author><name sortKey="Goldman Rakic, P S" uniqKey="Goldman Rakic P">P. S. Goldman-Rakic</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wertheimer, M" uniqKey="Wertheimer M">M. Wertheimer</name>
</author>
<author><name sortKey="Spillmann, L" uniqKey="Spillmann L">L. Spillmann</name>
</author>
<author><name sortKey="Sarris, V" uniqKey="Sarris V">V. Sarris</name>
</author>
<author><name sortKey="Sekuler, R" uniqKey="Sekuler R">R. Sekuler</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wilson, H R" uniqKey="Wilson H">H. R. Wilson</name>
</author>
<author><name sortKey="Cowan, J" uniqKey="Cowan J">J. Cowan</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wright, J J" uniqKey="Wright J">J. J. Wright</name>
</author>
<author><name sortKey="Kydd, R" uniqKey="Kydd R">R. Kydd</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wright, J J Liley D T J" uniqKey="Wright J">J. J. Liley, D. T. J. Wright</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Front Neuroinform</journal-id>
<journal-id journal-id-type="iso-abbrev">Front Neuroinform</journal-id>
<journal-id journal-id-type="publisher-id">Front. Neuroinform.</journal-id>
<journal-title-group><journal-title>Frontiers in Neuroinformatics</journal-title>
</journal-title-group>
<issn pub-type="epub">1662-5196</issn>
<publisher><publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">26539105</article-id>
<article-id pub-id-type="pmc">4611063</article-id>
<article-id pub-id-type="doi">10.3389/fninf.2015.00025</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Neuroscience</subject>
<subj-group><subject>Methods</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed</article-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Nichols</surname>
<given-names>Eric J.</given-names>
</name>
<xref ref-type="author-notes" rid="fn001"><sup>*</sup>
</xref>
<uri xlink:type="simple" xlink:href="http://loop.frontiersin.org/people/147472/overview"></uri>
</contrib>
<contrib contrib-type="author"><name><surname>Hutt</surname>
<given-names>Axel</given-names>
</name>
<uri xlink:type="simple" xlink:href="http://loop.frontiersin.org/people/5234/overview"></uri>
</contrib>
</contrib-group>
<aff><institution>Team Neurosys, Loria, Centre National de la Recherche Scientifique, INRIA, UMR no. 7503, Université de Lorraine</institution>
<country>Nancy, France</country>
</aff>
<author-notes><fn fn-type="edited-by"><p>Edited by: Gaute T. Einevoll, Norwegian University of Life Sciences, Norway</p>
</fn>
<fn fn-type="edited-by"><p>Reviewed by: Viktor Jirsa, Aix-Marseille University, France; Sid Visser, University of Exeter, UK</p>
</fn>
<corresp id="fn001">*Correspondence: Eric J. Nichols <email xlink:type="simple">ericjnichols@gmail.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub"><day>20</day>
<month>10</month>
<year>2015</year>
</pub-date>
<pub-date pub-type="collection"><year>2015</year>
</pub-date>
<volume>9</volume>
<elocation-id>25</elocation-id>
<history><date date-type="received"><day>20</day>
<month>5</month>
<year>2015</year>
</date>
<date date-type="accepted"><day>02</day>
<month>10</month>
<year>2015</year>
</date>
</history>
<permissions><copyright-statement>Copyright © 2015 Nichols and Hutt.</copyright-statement>
<copyright-year>2015</copyright-year>
<copyright-holder>Nichols and Hutt</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract><p>Neural Field models (NFM) play an important role in the understanding of neural population dynamics on a mesoscopic spatial and temporal scale. Their numerical simulation is an essential element in the analysis of their spatio-temporal dynamics. The simulation tool described in this work considers scalar spatially homogeneous neural fields taking into account a finite axonal transmission speed and synaptic temporal derivatives of first and second order. A text-based interface offers complete control of field parameters and several approaches are used to accelerate simulations. A graphical output utilizes video hardware acceleration to display running output with reduced computational hindrance compared to simulators that are exclusively software-based. Diverse applications of the tool demonstrate breather oscillations, static and dynamic Turing patterns and activity spreading with finite propagation speed. The simulator is open source to allow tailoring of code and this is presented with an extension use case.</p>
</abstract>
<kwd-group><kwd>neural field</kwd>
<kwd>numerical simulation</kwd>
<kwd>delay</kwd>
<kwd>pattern formation</kwd>
</kwd-group>
<counts><fig-count count="8"></fig-count>
<table-count count="1"></table-count>
<equation-count count="15"></equation-count>
<ref-count count="66"></ref-count>
<page-count count="11"></page-count>
<word-count count="6324"></word-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>1. Introduction</title>
<p>The understanding of spatio-temporal electric activity in neural tissue is essential in the study of neurobiological phenomena. To achieve this, mesoscopic-scale models such as neural mass and neural fields (NFM) which describe the dynamics of a large population of neurons reflecting coarse-grained properties of single neurons (Wilson and Cowan, <xref rid="B64" ref-type="bibr">1973</xref>
; Deco et al., <xref rid="B16" ref-type="bibr">2008</xref>
; Bressloff, <xref rid="B4" ref-type="bibr">2012</xref>
; Hutt and Buhry, <xref rid="B27" ref-type="bibr">2014</xref>
) play an important role. NFMs serve as a good description of the dynamic source of Local Field Potentials and encephalographic data (Nunez, <xref rid="B39" ref-type="bibr">1974</xref>
, <xref rid="B40" ref-type="bibr">2000</xref>
; Wright and Kydd, <xref rid="B65" ref-type="bibr">1992</xref>
; Wright and Liley, <xref rid="B66" ref-type="bibr">1994</xref>
; Jirsa et al., <xref rid="B34" ref-type="bibr">2002</xref>
; Nunez and Srinivasan, <xref rid="B41" ref-type="bibr">2006</xref>
; Coombes et al., <xref rid="B12" ref-type="bibr">2014</xref>
). They allow to consider diverse single neuron features that may tune neural population dynamics, such as somatic (Molaee-Ardekani et al., <xref rid="B38" ref-type="bibr">2007</xref>
) and synaptic adaptation (Coombes and Owen, <xref rid="B14" ref-type="bibr">2005</xref>
; Kilpatrick and Bressloff, <xref rid="B35" ref-type="bibr">2010</xref>
), extra-synaptic receptor dynamics (Hashemi et al., <xref rid="B22" ref-type="bibr">2014</xref>
; Hutt and Buhry, <xref rid="B27" ref-type="bibr">2014</xref>
) or finite axonal transmission speed (Jirsa and Haken, <xref rid="B33" ref-type="bibr">1996</xref>
; Pinto and Ermentrout, <xref rid="B46" ref-type="bibr">2001</xref>
; Hutt et al., <xref rid="B26" ref-type="bibr">2003</xref>
, <xref rid="B31" ref-type="bibr">2008</xref>
; Coombes, <xref rid="B11" ref-type="bibr">2005</xref>
; Faye and Faugeras, <xref rid="B19" ref-type="bibr">2010</xref>
; Veltz and Faugeras, <xref rid="B59" ref-type="bibr">2011</xref>
, <xref rid="B60" ref-type="bibr">2013</xref>
). All these applications make NFMs valuable in order to understand spatio-temporal dynamics of neural population activity.</p>
<p>Mathematical analysis and the numerical integration of NFMs are complementary. The recent years have shown strong attention of research on the mathematical properties of NFMs, whereas the numerical simulation of NFM solutions has been less considered in research. Since NFMs generalize partial differential equations (Coombes et al., <xref rid="B15" ref-type="bibr">2007</xref>
; Hutt, <xref rid="B23" ref-type="bibr">2007</xref>
) while involving finite transmission delay interactions, they allow to study a large class of pattern forming systems, cf. Hutt (<xref rid="B23" ref-type="bibr">2007</xref>
). In recent years, several software tools have been developed to simulate neural network dynamics. Examples for simulators for networks of spiking neurons are <italic>BRIAN</italic>
(Stimberg et al., <xref rid="B56" ref-type="bibr">2014</xref>
) and <italic>Neuron</italic>
(Carnevale and Hines, <xref rid="B10" ref-type="bibr">2006</xref>
). <italic>The Virtual Brain</italic>
(Sanz Leon et al., <xref rid="B53" ref-type="bibr">2013</xref>
) allows to simulate networks of neural mass models to reproduce global brain activity. The simulation platform <italic>DANA</italic>
(Rougier and Fix, <xref rid="B50" ref-type="bibr">2012</xref>
) simulates a hierarchy of coupled Dynamic Neural Fields which are decentralized, i.e., are updated numerically in time asynchronously (Rougier and Hutt, <xref rid="B51" ref-type="bibr">2011</xref>
). These latter software tools are powerful, general and highly adaptive to the framework of their neural network types. However, they do not provide the effective computation for the specific NFM given in Equation (1) which is a stochastic delayed integral-differential equation in two spatial dimensions. The tool presented here fills a gap in the landscape of neural simulator tools which are typically very general and adaptive and, hence, not efficient for NFM. A simulation tool for NFM allows to explore rapidly and in a user-friendly way the solution space of Equation (1) in order to reproduce numerically experimental spatio-temporal dynamics, e.g., to understand neuroimaging data (Friston et al., <xref rid="B21" ref-type="bibr">2014</xref>
; Pinotsis and Friston, <xref rid="B43" ref-type="bibr">2014</xref>
), retrieve neural sources and lateral connections (Pinotsis et al., <xref rid="B45" ref-type="bibr">2013</xref>
), and understand power spectra of electroencephalographic activity (Pinotsis et al., <xref rid="B44" ref-type="bibr">2012</xref>
). In addition, the tool promises to allow detection of new numerical solutions, cf. Section 3. The numerical analysis is non-trivial and challenging if NFMs become more complex, e.g., by involving complex dynamical features rendering the model high-dimensional or by considering delayed interactions. The present work considers a two-dimensional spatial embedding of neural populations similar to several previous studies (Laing, <xref rid="B36" ref-type="bibr">2005</xref>
; Owen et al., <xref rid="B42" ref-type="bibr">2007</xref>
) while taking into account finite axonal transmission speed (Hutt and Rougier, <xref rid="B30" ref-type="bibr">2010</xref>
, <xref rid="B29" ref-type="bibr">2014</xref>
). By virtue of its modularity, the simulator allows subsequent extensions with additional features, such as extra-synaptic receptor effects or several interacting populations.</p>
<p>The combination of finite axonal transmission speed and two-dimensional spatial embedding is challenging from a numerical simulation perspective due to the missing convolution structure (Hutt and Rougier, <xref rid="B30" ref-type="bibr">2010</xref>
, <xref rid="B29" ref-type="bibr">2014</xref>
) leading to long simulation durations. To overcome this problem, a numerical technique has been developed in recent years (Hutt and Rougier, <xref rid="B30" ref-type="bibr">2010</xref>
, <xref rid="B29" ref-type="bibr">2014</xref>
). Since future research in neural fields will investigate spatio-temporal dynamics involving finite axonal transmission speed, we have developed an open-source simulation toolbox that allows to gain spatio-temporal solutions of NFM models in two spatial dimensions, visualize them and save them, if necessary, as movies. We hope that the tool will provide an essential tool for the computational neuroscience community to advance the research field and the insight into the brain.</p>
<p>The simulator in this work obeys integral-differential equations of the type</p>
<disp-formula id="E1"><label>(1)</label>
<mml:math id="M1"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mi>η</mml:mi>
<mml:mfrac><mml:mrow><mml:msup><mml:mo>∂</mml:mo>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow><mml:mo>∂</mml:mo>
<mml:msup><mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mi>γ</mml:mi>
<mml:mfrac><mml:mo>∂</mml:mo>
<mml:mrow><mml:mo>∂</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mi>V</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>=</mml:mo>
<mml:mi>I</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>+</mml:mo>
<mml:mstyle displaystyle="true"><mml:mrow><mml:msub><mml:mo>∫</mml:mo>
<mml:mi>Ω</mml:mi>
</mml:msub>
<mml:mtext> </mml:mtext>
</mml:mrow>
</mml:mstyle>
<mml:mi>K</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo stretchy="false">)</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mi>S</mml:mi>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mi>V</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>−</mml:mo>
<mml:mfrac><mml:mrow><mml:mo>∥</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>∥</mml:mo>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:msup><mml:mtext>d</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>y</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>with a two-dimensional square spatial domain Ω and periodic boundary condition. The mean neuron potential <italic>V</italic>
∈ <inline-formula><mml:math id="M2"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
at location <bold>x</bold>
∈ Ω is evolved by the external stimulus <italic>I</italic>
(<italic>x, t</italic>
) ∈ <inline-formula><mml:math id="M3"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
and the integral of the synaptic connectivity kernel <italic>K</italic>
: <inline-formula><mml:math id="M4"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
<sup>2</sup>
→ <inline-formula><mml:math id="M5"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
and population firing rate <italic>S</italic>
∈ <inline-formula><mml:math id="M6"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
which depend on the distance between spatial locations <bold>x</bold>
and <bold>y</bold>
with a finite axonal transmission speed <italic>c</italic>
. Equation (1) represents the core of most NFM in the sense that most NFMs consider extensions of this equation.</p>
<p>Motivation for the work arises from a need for a visualization tool that is useful to the largest number of NFM researchers, allows for the tailoring of code and has fast while visually appealing output. The simulator can operate on all major operating systems. Output of data in three dimensions is provided by PyOpenGL which brings the speed and graphical detail of low-level OpenGL to the agile Python language. It is open source, enabling modification of the simulator in any beneficial way.</p>
</sec>
<sec id="s2"><title>2. Materials and methods</title>
<p>The cross-platform simulator is written in Python (version 2.7) and uses the NumPy library in consideration of its speed being close to the computational rate of the platform-dependent C language (Langtangen, <xref rid="B37" ref-type="bibr">2006</xref>
). The simulator can be downloaded<xref ref-type="fn" rid="fn0001"><sup>1</sup>
</xref>
in a package along with documentation<xref ref-type="fn" rid="fn0002"><sup>2</sup>
</xref>
describing its installation, running, features, and examples and the code is registered in ModelDB<xref ref-type="fn" rid="fn0003"><sup>3</sup>
</xref>
.</p>
<p>The following Section 2.1 describes the comprehensive access to field parameters, the subsequent section details the techniques applied to accelerate the simulation and Section 2.3 discusses the 3D visualization.</p>
<sec><title>2.1. Field parameters</title>
<p>A textual interface named values.py is provided in the root directory of the simulator code. It allows field values to be changed without knowledge of the inner workings of the simulator. For example, if η in Equation (1) is initialized as a non-zero number, a second order derivative is calculated to solve <italic>V</italic>
. Conversely, the interface eliminates the knowledge requirement of the numerical implementation of the derivatives and all other underlying code implementations. The interface has additional benefits of easily modifying variables in one place without searching through the code. This implementation permits changing parameters easily and sharing code amongst others working with similar simulations by the exchange of a single file.</p>
<p>The most important aspect of a text-based interface from its user experience is its facilitation of novelty by allowing absolute control of all terms of Equation (1). For instance, matrix <italic>I</italic>
can be defined in the interface with as many lines of Python code as necessary given the definition ends with an assignment (i.e., <italic>I = …</italic>
). Assigning the first parameter in the values.py file, named showData, a value of 3 displays <italic>I</italic>
in the simulator, which can be useful when refining its values. Time-varying spatio-temporal input is available in the interface by uncommenting and modifying the body of a function named <italic>updateI</italic>
in any manner while maintaining that <italic>I</italic>
is returned. Neural field investigations are thereby efficiently implemented with free choice over all the variables accessible through the interface while retaining the full performance.</p>
</sec>
<sec><title>2.2. Accelerated simulation</title>
<p>The simulator is advantageous in its acceleration of spatial and temporal integration. Multiple approaches are used to increase the simulation speed.</p>
<sec><title>2.2.1. Spatial and temporal integration</title>
<p>Equation (1) includes a spatial integral with a homogeneous kernel <italic>K</italic>
. Please recall that homogeneous kernels just depend on the difference vector <bold>x</bold>
− <bold>y</bold>
between two spatial locations <bold>x, y</bold>
including isotropic kernels, i.e., <italic>K</italic>
= <italic>K</italic>
(∥<bold>x</bold>
− <bold>y</bold>
∥), as a specific case. In the absence of the finite transmission delay term, this integral would represent a spatial convolution and would be solvable numerically efficiently by a Fast Fourier Transform (FFT) (Van Loan, <xref rid="B58" ref-type="bibr">1991</xref>
). For non-vanishing transmission delay, the convolution structure is less obvious and the FFT is not applicable directly. Nevertheless, it is possible to re-write the spatial integration to utilize a FFT in space (Owen et al., <xref rid="B42" ref-type="bibr">2007</xref>
; Hutt and Rougier, <xref rid="B30" ref-type="bibr">2010</xref>
, <xref rid="B29" ref-type="bibr">2014</xref>
) as</p>
<disp-formula id="E2"><mml:math id="M7"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:mrow><mml:msub><mml:mo>∫</mml:mo>
<mml:mi>Ω</mml:mi>
</mml:msub>
<mml:mrow><mml:msup><mml:mtext>d</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi>y</mml:mi>
<mml:mtext> </mml:mtext>
<mml:mi>K</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mi>S</mml:mi>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mi>V</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>−</mml:mo>
<mml:mfrac><mml:mrow><mml:mo>∥</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>∥</mml:mo>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:mrow><mml:msub><mml:mo>∫</mml:mo>
<mml:mi>Ω</mml:mi>
</mml:msub>
<mml:mrow><mml:msup><mml:mtext>d</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi>y</mml:mi>
<mml:msubsup><mml:mstyle displaystyle="true"><mml:mo>∫</mml:mo>
</mml:mstyle>
<mml:mn>0</mml:mn>
<mml:mrow><mml:msub><mml:mi>τ</mml:mi>
<mml:mi>m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mi>d</mml:mi>
<mml:mi>τ</mml:mi>
<mml:mi>L</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>τ</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mi>S</mml:mi>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mi>V</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>τ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>with the maximum delay time τ<sub><italic>m</italic>
</sub>
and the spatio-temporal kernel function <italic>L</italic>
(<bold>x</bold>
, <italic>t</italic>
) = <italic>K</italic>
(<bold>x</bold>
)δ(∥ <bold>x</bold>
∥ ∕<italic>c</italic>
− <italic>t</italic>
). We observe that the spatial summation represents an integration over delayed spatial rings, which are convolved spatially with the transfer function <italic>S</italic>
in Equation (1). Introducing a regular rectangular spatial grid for spatial discretization, finite axonal speed <italic>c</italic>
yields rings of width</p>
<disp-formula id="E3"><label>(2)</label>
<mml:math id="M8"><mml:mrow><mml:mi>w</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo> </mml:mo>
<mml:mi>c</mml:mi>
<mml:mo>·</mml:mo>
<mml:mi>Δ</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>·</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>l</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>delineated within the field, where <italic>n</italic>
and <italic>l</italic>
are the number of discretized spatial units and the length of the field, respectively, and Δ<italic>t</italic>
is the finite integration time step. The Pythagorean theorem gives the maximum radius of the rings in the field <inline-formula><mml:math id="M9"><mml:mi>r</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>∕</mml:mo>
<mml:msqrt><mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msqrt>
</mml:math>
</inline-formula>
over which the spatial integration is performed, which is applied to obtain the number of rings in a field as</p>
<disp-formula id="E4"><label>(3)</label>
<mml:math id="M10"><mml:msub><mml:mi>n</mml:mi>
<mml:mrow><mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mrow><mml:mo>⌊</mml:mo>
<mml:mrow><mml:mi>r</mml:mi>
<mml:mo>/</mml:mo>
<mml:mi>w</mml:mi>
</mml:mrow>
<mml:mo>⌋</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mrow><mml:mo>⌊</mml:mo>
<mml:mrow><mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:msqrt><mml:mn>2</mml:mn>
</mml:msqrt>
<mml:mi>c</mml:mi>
<mml:mi>Δ</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo>⌋</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>defining the maximum delay to τ<sub><italic>m</italic>
</sub>
= <italic>n</italic>
<sub><italic>rings</italic>
</sub>
Δ<italic>t</italic>
. The spatio-temporal kernel <italic>L</italic>
is determined by the spatial kernel <italic>K</italic>
and the axonal speed <italic>c</italic>
(Hutt and Rougier, <xref rid="B30" ref-type="bibr">2010</xref>
, <xref rid="B29" ref-type="bibr">2014</xref>
).</p>
<p>Equation (1) involves distance-dependent delays which represent a specific type of distributed delays (Hutt and Lefebvre, <xref rid="B28" ref-type="bibr">in press</xref>
). To this end, it is necessary to initialize the field variable <italic>V</italic>
in an initial time interval and the toolbox allows the user to set the initial values arbitrarily. The external input <italic>I</italic>
may be deterministic or stochastic and the user may choose it according to her needs, e.g., implementing spatial correlations in stochastic inputs. To integrate the evolution equation in time the user may choose between different integration methods for delay differential equations (Buckwar and Winkler, <xref rid="B7" ref-type="bibr">2006</xref>
, <xref rid="B8" ref-type="bibr">2007</xref>
). Standard methods discretize time regularly in steps of duration Δ<italic>t</italic>
yielding results (Equations 2, 3). In the case of stochastic input, the toolbox includes numerical implementations of the delayed Euler-Maruyama method (Buckwar et al., <xref rid="B6" ref-type="bibr">2008</xref>
) and the stochastic version of the Runge-Kutta method for delayed differential equation (Carletti, <xref rid="B9" ref-type="bibr">2006</xref>
). For deterministic inputs, the equivalent deterministic methods are available.</p>
<p>If there is no modification to <italic>K</italic>
and <italic>c</italic>
during the simulation, then <italic>L</italic>
is calculated once only before the start of the simulation while <italic>S</italic>
(·) changes over time. The convolution of <italic>L</italic>
and <italic>S</italic>
is performed using a FFT what greatly increases the speed of the integral convolution compared to conventional integration. This can be understood easily recalling that the two-dimensional FFT needs to sum up <inline-formula><mml:math id="M11"><mml:msup><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msubsup><mml:mrow><mml:mo class="qopname">log</mml:mo>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>
terms leading to summands of the total number of <inline-formula><mml:math id="M12"><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>F</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msup><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msubsup><mml:mrow><mml:mo class="qopname">log</mml:mo>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mrow><mml:mo stretchy="false">(</mml:mo>
<mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
. In contrast, conventional integration sums up terms of number <italic>n</italic>
<sup>4</sup>
for each delay time and hence the total number of computation <inline-formula><mml:math id="M13"><mml:msub><mml:mrow><mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msup><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>4</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msub><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
, cf. Appendix I. Hence the FFT implementation speeds up the integration by a factor of</p>
<disp-formula id="E5"><label>(4)</label>
<mml:math id="M14"><mml:mrow><mml:msub><mml:mi>f</mml:mi>
<mml:mrow><mml:mi>s</mml:mi>
<mml:mi>p</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>u</mml:mi>
<mml:mi>p</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:mi>F</mml:mi>
<mml:mi>F</mml:mi>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow><mml:msub><mml:mi>N</mml:mi>
<mml:mrow><mml:mi>c</mml:mi>
<mml:mi>o</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>=</mml:mo>
<mml:mfrac><mml:mrow><mml:msup><mml:mi>n</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow><mml:msubsup><mml:mrow><mml:mi>log</mml:mi>
</mml:mrow>
<mml:mn>2</mml:mn>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<p>The axonal speed's implementation is described in detail in Hutt and Rougier (<xref rid="B30" ref-type="bibr">2010</xref>
). It is important to note that other (conventional) numerical software tools not taking into account the convolution structure have to sum up <italic>N</italic>
<sub><italic>conv</italic>
</sub>
terms in case of fully connected networks. For instance, this holds true for the simulation tools <italic>BRIAN, Neuron</italic>
and <italic>The Virtual Brain</italic>
(Sanz Leon et al., <xref rid="B52" ref-type="bibr">2015</xref>
) which have to memorize the history and advance the stored activity field <italic>n</italic>
<sub><italic>rings</italic>
</sub>
times. We also note that the method proposed may be implemented in these simulation tools in the future since they also may consider spatially homogeneous neural fields as specific cases. The FFT-based method presented here computes the network interactions faster than these tools by <italic>f</italic>
<sub><italic>speedup</italic>
</sub>
given in Equation (4). For instance, for a typical number of spatial grid intervals of <italic>n</italic>
= 512 as used in the application showed in Section 3, we obtain the huge speed up factor of <italic>f</italic>
<sub><italic>speedup</italic>
</sub>
= ≈3236.</p>
</sec>
<sec><title>2.2.2. Self-writing code</title>
<p>The second approach to increase the simulation speed employs self-writing code to reduce the simulator's instruction set. The simulator writes and executes its own code to increase the efficiency of simulations and display only the user defined features. The simulation code is based on interface selections and is self-written by an initialization module at the onset of the program. The interface offers features, such as a second derivative calculation, <italic>I</italic>
and <italic>K</italic>
updates and added noise, that conditionally run during the simulation and are not performed over time if the user selects to view <italic>V</italic>
, <italic>I</italic>
, or <italic>K</italic>
at <italic>t</italic>
= 0. For example, the visual interface offers the choice of viewing the spatial kernel <italic>K</italic>
for its design and visualization. Only the code that initializes and displays <italic>K</italic>
is written to the executing module if this choice is selected. The self-writing code is also favorable when the full simulation is run with calculations executed unconditionally. The result is very efficient code that is changed with every modification to the interface.</p>
</sec>
<sec><title>2.2.3. Implementation on GPUs</title>
<p>The third approach parallelizes the output calculations on the graphics processing unit (GPU). The displayed matrix is put onto the running system's GPU for hardware acceleration of the visualization. Vertex buffer objects improve visualization throughput by uploading vertices to video device memory where vertex and fragment shaders transform and write neural field data in parallel to the framebuffer for display. The simulator also avoids the CPU to GPU information transfer bottleneck by its storage of data on the video device memory. This is accomplished with the <italic>OpenGL Shading Language</italic>
that is used through PyOpenGL to achieve a better visual description of information than is provided in other tools. A background on PyOpenGL and its comparison to other visualization libraries can be found in Rossant and Harris (<xref rid="B49" ref-type="bibr">2013</xref>
).</p>
</sec>
<sec><title>2.2.4. Optimal visualization rate</title>
<p>The fourth approach to accelerate simulations is to display field matrices at a rate optimized for continuous visualization perception. Two images are perceived simultaneously when there is an interval of less than 30 ms between them (Wertheimer et al., <xref rid="B63" ref-type="bibr">2012</xref>
). The simulator takes advantage of the temporal lag in biological visual perception by stopping the numerical calculations to submit the field data to the GPU once within every 30 ms. This allows for the numerical part of the simulations to continue with fewer stoppages, resulting in faster simulations.</p>
</sec>
</sec>
<sec><title>2.3. 3D visualization</title>
<p>The open source and cross-platform <italic>show3D</italic>
library was written for the Neural Field Simulator to display field information. The library's visualization of neural fields expands two dimensional neural field data into a third dimension to better observe the differences in field locations. This is achieved by raising every value in the 2D spatial plane to a third dimension position [<italic>x</italic>
, <italic>y</italic>
] ↦ [<italic>z</italic>
] relative to other 2D field values.</p>
<p>Color values are efficiently manipulated with the keyboard keys shown in Appendix II (Table <xref ref-type="table" rid="TA1">A1</xref>
). There is a selection of 8 colors, cf. Figure <xref ref-type="fig" rid="F1">1</xref>
, available for the background, minimum, middle, and maximum graph values.</p>
<fig id="F1" position="float"><label>Figure 1</label>
<caption><p><bold>Selection of colors that can be applied to the background and ranges of the plotted matrix</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0001"></graphic>
</fig>
<p>Intermediate color transformations are encoded in a dictionary containing 8<sup>2</sup>
unique 3 element vectors, each representing red, green, and blue colors. The appropriate color transformation vector is uploaded to the GPU where the vector elements represent one mutually inclusive index of [0, 1, <italic>Z</italic>
, 1-<italic>Z</italic>
] with Z axis locations ∈ <bold>R</bold>
∣ 0 < <italic>Z</italic>
< 1. Each location on the Z axis is subsequently colored in parallel by the GPU with the appropriate shade. Graph value colors are interpolated with two choices of ranges: [minimum, maximum] and [minimum, middle], [middle, maximum] graph value colors. Different depths of the graph can be highlighted by raising or lowering the ranges of colors.</p>
<p>Scrolling the mouse rotates neural fields in the direction of the mouse and the keyboard is used to move fields in various ways, cf. Figure <xref ref-type="fig" rid="F2">2</xref>
.</p>
<fig id="F2" position="float"><label>Figure 2</label>
<caption><p><bold>Keyboard keys and their corresponding movements</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0002"></graphic>
</fig>
<p>Images and videos of simulations are saved, respectively in .png and .mp4 formats by using the keyboard keys in Appendix II. Visualization parameters are saved by the library after every simulation to reduce graphical adjustments during subsequent simulations of neural fields.</p>
<p>The show3D graphical visualization library is not limited to neural field data. Every two dimensional NumPy matrix can be displayed in 3 dimensions using the show3D library. Documentation for the show3D library's use, including a tutorial and code API, is online<xref ref-type="fn" rid="fn0004"><sup>4</sup>
</xref>
and packaged with example code along with the library<xref ref-type="fn" rid="fn0005"><sup>5</sup>
</xref>
. However, there is no requirement for the library's separate download for use with the simulator because it is integrated into the Neural Field Simulator.</p>
</sec>
</sec>
<sec id="s3"><title>3. Applications</title>
<p>The simulator can be used to analyze spatio-temporal neural field dynamics. The simulator's open source code allows modifications and extensions to be added to the code. The subsequent sections describe few of these possible applications.</p>
<p>Introducing finite axonal transmission speed in neural fields substantially slows numerical computation. However, to omit finite transmission speed is to neglect biological physiology (Idiart and Abbott, <xref rid="B32" ref-type="bibr">1993</xref>
). Hutt and Rougier (<xref rid="B30" ref-type="bibr">2010</xref>
) have suggested to implement finite axonal transmission speed in a computationally efficient manner that is utilized by the simulator. Numerically, the speed is infinite if <inline-formula><mml:math id="M15"><mml:mi>c</mml:mi>
<mml:mo>≥</mml:mo>
<mml:mi>l</mml:mi>
<mml:mo>∕</mml:mo>
<mml:msqrt><mml:mrow><mml:mn>2</mml:mn>
</mml:mrow>
</mml:msqrt>
<mml:mi>Δ</mml:mi>
<mml:mi>t</mml:mi>
</mml:math>
</inline-formula>
and there is increasing delay as <italic>c</italic>
decreases.</p>
<sec><title>3.1. Breather</title>
<p>Breather oscillations have been reported in theory (Folias and Bressloff, <xref rid="B20" ref-type="bibr">2005</xref>
; Hutt and Rougier, <xref rid="B30" ref-type="bibr">2010</xref>
) and experiments (Wang, <xref rid="B61" ref-type="bibr">2010</xref>
). As shown in Hutt and Rougier (<xref rid="B30" ref-type="bibr">2010</xref>
), breathers are solutions of Equation (1) for finite axonal transmission speeds. They can be obtained and visualized in the simulator by assuming a temporally constant external input <italic>I</italic>
in Equation (1). For a Gaussian-shape input with its apex at the center of the field, one overwrites the <italic>I</italic>
variable section in the values.py file as</p>
<disp-formula id="E6"><mml:math id="M16"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sigma</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">5</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">65685425</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext>I </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">20</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mtext>x</mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sigma</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sigma</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>and change the showData variable assignment near the beginning of the values.py file to</p>
<disp-formula id="E7"><mml:math id="M17"><mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>showData</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">3</mml:mn>
</mml:mrow>
</mml:math>
</disp-formula>
<p>to show the input <italic>I</italic>
in the simulation. In the definition of <italic>I</italic>
, the space variable <italic>x</italic>
∈ <inline-formula><mml:math id="M18"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
<sup>2</sup>
is defined to cover the spatial domain Ω (not shown in the code snippet). A field input similar to Figure <xref ref-type="fig" rid="F3">3A</xref>
can be seen when the simulator is run.</p>
<fig id="F3" position="float"><label>Figure 3</label>
<caption><p><bold>Breather parameters plotted in the simulator for (A) <italic>I</italic>
and (B) <italic>K</italic>
in Equation (1)</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0003"></graphic>
</fig>
<p>An inhibitory synaptic connectivity kernel, <italic>K</italic>
in Equation (1), can be implemented for the breather and viewed by changing the <italic>showData</italic>
and <italic>K</italic>
variables in the values.py file to</p>
<disp-formula id="E8"><mml:math id="M19"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>showData</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">4</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">4</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mtext>x</mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mn mathvariant="monospace">3</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">18</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>and running the simulator. Here, <italic>x</italic>
∈ <inline-formula><mml:math id="M20"><mml:mi mathvariant="-tex-caligraphic">R</mml:mi>
</mml:math>
</inline-formula>
<sup>2</sup>
. An inhibitory synaptic kernel similar to the one in Figure <xref ref-type="fig" rid="F3">3B</xref>
can be subsequently viewed.</p>
<p>After overwriting <italic>I</italic>
and <italic>K</italic>
as noted above, replace the following variables and function in the values.py file with the values below:</p>
<disp-formula id="E9"><mml:math id="M21"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>showData</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>endTime</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">1</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>dt</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">002</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>gamma</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1.0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>eta</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">500</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>l</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">30</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">512</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>zeros</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>noiseVcont</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>a</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">32</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>b</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mn mathvariant="monospace">32</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">32</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">0.1</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sqrt</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>dt</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mrow></mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>def</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>updateS</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>return</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">10000</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">005</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>Spatially localized breather oscillations are replicated by running the program. Figure <xref ref-type="fig" rid="F4">4</xref>
shows two cycles of the oscillations after setting the minimum and maximum z axis values by typing</p>
<disp-formula id="E10"><mml:math id="M22"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace"><mml:mtext>n </mml:mtext>
</mml:mstyle>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">00</mml:mn>
<mml:mstyle mathvariant="monospace"><mml:mtext>48 </mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace"><mml:mtext>Enter key</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace"><mml:mtext>y </mml:mtext>
</mml:mstyle>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">00</mml:mn>
<mml:mstyle mathvariant="monospace"><mml:mtext>58 </mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace"><mml:mtext>Enter key</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>after running the program.</p>
<fig id="F4" position="float"><label>Figure 4</label>
<caption><p><bold>Two cycles of the breather oscillations</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0004"></graphic>
</fig>
</sec>
<sec><title>3.2. Turing patterns</title>
<p>Turing patterns (Turing, <xref rid="B57" ref-type="bibr">1952</xref>
) have been reported in neural field models (Atay and Hutt, <xref rid="B1" ref-type="bibr">2006</xref>
; Elvin et al., <xref rid="B18" ref-type="bibr">2009</xref>
; Steyn-Ross et al., <xref rid="B55" ref-type="bibr">2010</xref>
). The Neural Field Simulator is able to compute and display noisy neural field activity evolving into Turing patterns.</p>
<p>Static Turing patterns emerge from noisy initial conditions with the following interface properties:</p>
<disp-formula id="E11"><mml:math id="M23"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>showData</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>endTime</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">10</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>dt</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">01</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>gamma</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>eta</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">6364</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>l</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">90</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">512</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>ones</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">5</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mtext>4 </mml:mtext>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>random</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>normal</mml:mtext>
</mml:mstyle>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mn mathvariant="monospace">0.1</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>noiseVcont</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>None</mml:mtext>
</mml:mstyle>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>I</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>zeros</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>lins</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>linspace</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">9</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">1</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>zeros</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>for</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>in</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>range</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">:</mml:mo>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>lins</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>I</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mn mathvariant="monospace">150</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>for</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>in</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>range</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>lins</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mn mathvariant="monospace">200</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mrow></mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>def</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>updateS</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>return</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">24</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">3</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>With <italic>c</italic>
= 6364 the effective speed is infinite for the given <italic>l</italic>
and <italic>dt</italic>
properties. Figure <xref ref-type="fig" rid="F5">5</xref>
shows the simulation starting with random field potential noise. A pattern begins to emerge at 1 s and evolves into a temporally constant Turing pattern after approximately 5 s.</p>
<fig id="F5" position="float"><label>Figure 5</label>
<caption><p><bold>Static Turing pattern emerging over 5 s</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0005"></graphic>
</fig>
<p>Dynamic Turing patterns emerging over time in the simulator with interface values:</p>
<disp-formula id="E12"><mml:math id="M24"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>showData</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> 1 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>endTime</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">40</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>dt</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">00</mml:mn>
<mml:mtext>5 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>gamma</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0.</mml:mn>
<mml:mtext>82 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>eta</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">10</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>l</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">10</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> 256 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>ones</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">4</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mtext>1 </mml:mtext>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>random</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>normal</mml:mtext>
</mml:mstyle>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mn mathvariant="monospace">0.1</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>noiseVcont</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>None</mml:mtext>
</mml:mstyle>
</mml:mstyle>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>Uexcite</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>zeros</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>I</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>ones</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>lins</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>linspace</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">7</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mtext>1 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>localStrong</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>linspace</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>zeros</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>for</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>in</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>range</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">:</mml:mo>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>lins</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>localStrong</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>for</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>in</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>range</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>lins</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>localStrong</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>i</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mrow></mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>def</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>updateS</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>return</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">5</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">5</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">3</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>Figure <xref ref-type="fig" rid="F6">6</xref>
shows a typical simulation, given random starting field potential noise, with different Turing patterns materializing. Activity patterns form at varying intervals, generally every few seconds, throughout the simulation. The times in Figure <xref ref-type="fig" rid="F6">6</xref>
were chosen for clear displays of different (vertical, cone, and horizontal) patterns.</p>
<fig id="F6" position="float"><label>Figure 6</label>
<caption><p><bold>Turing patterns in neural field activity forming over time during the same simulation</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0006"></graphic>
</fig>
</sec>
<sec><title>3.3. Finite spreading speed</title>
<p>Stimulating a neural population at a single location, as is done in typical physiological experiments applying external stimuli, the neural activity spreads in the population. Since finite transmission speed represents delayed spatial interaction in the population under study, it affects the spreading speed of activity (Hutt, <xref rid="B23" ref-type="bibr">2007</xref>
). If the transmission speed is infinite, the activity spreads diffusively involving the instantaneous activation at all spatial locations. Conversely, finite transmission speed delays the activity spread leading to a slowly-moving spreading front (Hutt and Atay, <xref rid="B25" ref-type="bibr">2006</xref>
; Hutt, <xref rid="B24" ref-type="bibr">2009</xref>
). Figure <xref ref-type="fig" rid="F7">7</xref>
shows numerical simulations for large (top row) and small speeds (bottom row), other parameters are identical.</p>
<fig id="F7" position="float"><label>Figure 7</label>
<caption><p><bold>Activity spread with large speed <italic>c</italic>
(top) and small speed <italic>c</italic>
(bottom)</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0007"></graphic>
</fig>
<p>The simulator allows the transmission speed to be examined closely in the visualization window by decreasing the maximum z axis value to be close to the original field value. This was done in Figure <xref ref-type="fig" rid="F7">7</xref>
by typing</p>
<disp-formula id="E13"><mml:math id="M25"><mml:mrow><mml:mstyle mathvariant="monospace"><mml:mtext>y </mml:mtext>
</mml:mstyle>
<mml:mn mathvariant="monospace">2.0000</mml:mn>
<mml:mstyle mathvariant="monospace"><mml:mtext>2 </mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">[</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace"><mml:mtext>Enter key</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>after starting the simulator and before beginning the simulation.</p>
<p>The parameters to achieve the results in Figure <xref ref-type="fig" rid="F7">7</xref>
are</p>
<disp-formula id="E14"><mml:math id="M26"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>showData</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> 1 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>endTime</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> 1 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>dt</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">00</mml:mn>
<mml:mtext>4 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>gamma</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>eta</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0.</mml:mn>
<mml:mtext>35 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">10</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>l</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">10</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> 256 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>ones</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>noiseVcont</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>None</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>Uexcite</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>zeros</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">,</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>I</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>x</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">04</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mn mathvariant="monospace">0</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">04</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>phi</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext>3 </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">_</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">10</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>pi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>l</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">=</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">0.1</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>cos</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">_</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>a</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">\</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>cos</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>k</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">_</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>a</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>cos</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>phi</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>b</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>phi</mml:mtext>
</mml:mstyle>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace" stretchy="false">)</mml:mo>
<mml:mtext>)</mml:mtext>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">\</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>cos</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>K</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">_</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>c</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>a</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>cos</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>phi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>b</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>sin</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>phi</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace" stretchy="false">)</mml:mo>
<mml:mo mathvariant="monospace" stretchy="false">)</mml:mo>
<mml:mo mathvariant="monospace" stretchy="false">)</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">\</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">−</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>x</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mn mathvariant="monospace">10</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>l</mml:mtext>
</mml:mstyle>
</mml:mstyle>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>float</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>n</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mn mathvariant="monospace">2</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mrow></mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>def</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>updateS</mml:mtext>
</mml:mstyle>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace">:</mml:mo>
<mml:mtext> </mml:mtext>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mtext> </mml:mtext>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="blue"><mml:mtext>return</mml:mtext>
</mml:mstyle>
<mml:mtext> </mml:mtext>
<mml:mn mathvariant="monospace">2</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace">/</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mn mathvariant="monospace">1</mml:mn>
<mml:mo mathvariant="monospace">+</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>np</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>exp</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace" stretchy="false">(</mml:mo>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">5</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">5</mml:mn>
<mml:mo mathvariant="monospace">*</mml:mo>
<mml:mrow><mml:mo mathvariant="monospace">(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="monospace" class="text" mathcolor="red"><mml:mtext>V</mml:mtext>
</mml:mstyle>
<mml:mo mathvariant="monospace">−</mml:mo>
<mml:mn mathvariant="monospace">3</mml:mn>
<mml:mo mathvariant="monospace">.</mml:mo>
<mml:mn mathvariant="monospace">0</mml:mn>
</mml:mrow>
<mml:mo mathvariant="monospace">)</mml:mo>
</mml:mrow>
<mml:mo mathvariant="monospace" stretchy="false">)</mml:mo>
<mml:mo mathvariant="monospace" stretchy="false">)</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>where <italic>c</italic>
is chosen according to the values given in Figure <xref ref-type="fig" rid="F7">7</xref>
.</p>
</sec>
<sec><title>3.4. Extensions</title>
<p>The simulator, being open source, allows the tailoring of code to provide modifications and extensions. An example extension is the addition of a graphical interface to modify parameters and simulate neural fields. Figure <xref ref-type="fig" rid="F8">8</xref>
shows an example interface coded with the wxPython<xref ref-type="fn" rid="fn0006"><sup>6</sup>
</xref>
toolkit. Simulations are started, paused, continued, and started anew by clicking a button.</p>
<fig id="F8" position="float"><label>Figure 8</label>
<caption><p><bold>A graphical interface for the simulator</bold>
.</p>
</caption>
<graphic xlink:href="fninf-09-00025-g0008"></graphic>
</fig>
<p>Neural field parameters can be modified prior to and during simulations by clicking on the appropriate area of the interface and completing a pop-up dialogue. Running simulations are automatically paused when a mouse hovers above parameter areas of the interface. There is a symbiosis among the show3D library discussed in Section 2.3 and the parameter selection interface. It is possible to view the external input stimulus, kernel, and firing rate in the GLFW window by adding a mouse event and hovering over these sections to automatically view the respective matrices. Viewing these elements while changing their parameters can help to fine-tune field parameters. Moving the mouse away from these areas unpauses paused simulations and the field potential matrix is shown in the GLFW window. Further synergy between the interface and show3D library is implemented with the option to alter graph values from the interface where z axis limits and colors can be selected.</p>
</sec>
</sec>
<sec id="s4"><title>4. Discussion</title>
<p>The Neural Field Simulator and its dependencies are cross-platform. However, the simulator interacts with graphics hardware using system-specific drivers which can result in problems on some operating systems. The graphical user interface in Section 3.4 is an example of this where the cross platform wxPython toolkit uses OpenGL to draw to the screen. The show3D library likewise uses OpenGL to interact with GPU. The graphical user interface and show3D library function symbiotically on Mac systems via the Apple Graphics Library. Conversely, on other operating systems such as Linux and Windows, unfortunately the separate utilization of OpenGL causes the simulator to crash. To this end, the current version of the simulator is released without the addition of extensions in order to operate properly on every major operating system. Nevertheless, a graphical interface can be a good choice with an appropriate single operating system.</p>
<p>Apart from the software implementation, in future work some model assumptions can be released. The square geometry can be recast easily to a rectangular geometry, whereas more general geometries (e.g., the impressive implementation work in <italic>The Virtual Brain</italic>
Sanz Leon et al., <xref rid="B52" ref-type="bibr">2015</xref>
) will take more numerical effort. Periodic boundary conditions guarantee the simple application of the FFT, effective implementations of other boundary conditions like Dirichlet conditions [<italic>V</italic>
(<bold>z</bold>
, <italic>t</italic>
) = const, <bold>z</bold>
∈ ∂Ω] will demand some implementation changes in the spatial integral computation. Such modifications may still retain the fundamental implementation of the FFT. In contrast, rejecting the homogeneity hypothesis of spatial interactions, i.e., <italic>K</italic>
= <italic>K</italic>
(<bold>x, y</bold>
) ≠ <italic>K</italic>
(<bold>x</bold>
− <bold>y</bold>
), abolishes the convolution structure and slows down the numerical simulation, cf. Appendix I.</p>
<p>Future work will extend the NFM model to multiple equations to render the model even more biologically plausible. In addition, an extension of the implementation to a mixture of constant and space-dependent delays as considered by Veltz and Faugeras (<xref rid="B59" ref-type="bibr">2011</xref>
) will be interesting.</p>
<sec><title>Conflict of interest statement</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
</sec>
</body>
<back><ack><p>The authors thank Kevin Green for valuable testing of the software under MacOS and Linux and for providing system parameters for Figure <xref ref-type="fig" rid="F6">6</xref>
. This work is funded by the European Research Council for support under the European Union's Seventh Framework Programme (FP7/2007-2013), ERC grant agreement No. 257253 (MATHANA project).</p>
</ack>
<fn-group><fn id="fn0001"><p><sup>1</sup>
<ext-link ext-link-type="uri" xlink:href="https://gforge.inria.fr/projects/nfsimulator/">https://gforge.inria.fr/projects/nfsimulator/</ext-link>
.</p>
</fn>
<fn id="fn0002"><p><sup>2</sup>
<ext-link ext-link-type="uri" xlink:href="http://nfsimulator.gforge.inria.fr">http://nfsimulator.gforge.inria.fr</ext-link>
.</p>
</fn>
<fn id="fn0003"><p><sup>3</sup>
<ext-link ext-link-type="uri" xlink:href="https://senselab.med.yale.edu/ModelDB/showModel.cshtml?model=184479">https://senselab.med.yale.edu/ModelDB/showModel.cshtml?model=184479</ext-link>
.</p>
</fn>
<fn id="fn0004"><p><sup>4</sup>
<ext-link ext-link-type="uri" xlink:href="http://show3d.gforge.inria.fr/index.html">http://show3d.gforge.inria.fr/index.html</ext-link>
.</p>
</fn>
<fn id="fn0005"><p><sup>5</sup>
<ext-link ext-link-type="uri" xlink:href="https://gforge.inria.fr/projects/show3d/">https://gforge.inria.fr/projects/show3d/</ext-link>
.</p>
</fn>
<fn id="fn0006"><p><sup>6</sup>
<ext-link ext-link-type="uri" xlink:href="http://www.wxpython.org">http://www.wxpython.org</ext-link>
.</p>
</fn>
</fn-group>
<ref-list><title>References</title>
<ref id="B1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Atay</surname>
<given-names>F. M.</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2006</year>
). <article-title>Neural fields with distributed transmission speeds and long-range feedback delays</article-title>
. <source>SIAM J. Appl. Dyn. Syst.</source>
<volume>5</volume>
, <fpage>670</fpage>
–<lpage>698</lpage>
. <pub-id pub-id-type="doi">10.1137/050629367</pub-id>
</mixed-citation>
</ref>
<ref id="B2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>beim Graben</surname>
<given-names>P.</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>Attractor and saddle node dynamics in heterogeneous neural fields</article-title>
. <source>EPJ Nonlin. Biomed. Phys.</source>
<volume>2</volume>
, <fpage>4</fpage>
<pub-id pub-id-type="doi">10.1140/epjnbp17</pub-id>
</mixed-citation>
</ref>
<ref id="B3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brackley</surname>
<given-names>C. A.</given-names>
</name>
<name><surname>Turner</surname>
<given-names>M.</given-names>
</name>
</person-group>
(<year>2009</year>
). <article-title>Persistent fluctuations of activity in undriven continuum neural field models with power-law connections</article-title>
. <source>Phys. Rev. E</source>
<volume>79</volume>
:<fpage>011918</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.79.011918</pub-id>
<pub-id pub-id-type="pmid">19257080</pub-id>
</mixed-citation>
</ref>
<ref id="B4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bressloff</surname>
<given-names>P.</given-names>
</name>
</person-group>
(<year>2012</year>
). <article-title>Spatiotemporal dynamics of continuum neural fields</article-title>
. <source>J. Phys. A</source>
<volume>45</volume>
:<fpage>033001</fpage>
<pub-id pub-id-type="doi">10.1088/1751-8113/45/3/033001</pub-id>
</mixed-citation>
</ref>
<ref id="B5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bressloff</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
(<year>2001</year>
). <article-title>Traveling fronts and wave propagation failure in an inhomogeneous neural network</article-title>
. <source>Phys. D</source>
<volume>155</volume>
, <fpage>83</fpage>
–<lpage>100</lpage>
. <pub-id pub-id-type="doi">10.1016/S0167-2789(01)00266-4</pub-id>
</mixed-citation>
</ref>
<ref id="B6"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buckwar</surname>
<given-names>E.</given-names>
</name>
<name><surname>Kuske</surname>
<given-names>R.</given-names>
</name>
<name><surname>Mohammed</surname>
<given-names>S.</given-names>
</name>
<name><surname>Shardlow</surname>
<given-names>T.</given-names>
</name>
</person-group>
(<year>2008</year>
). <article-title>Weak convergence of the euler scheme for stochastic differential delay equations</article-title>
. <source>LMS J. Comput. Math.</source>
<volume>11</volume>
, <fpage>60</fpage>
–<lpage>99</lpage>
. <pub-id pub-id-type="doi">10.1112/S146115700000053X</pub-id>
</mixed-citation>
</ref>
<ref id="B7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buckwar</surname>
<given-names>E.</given-names>
</name>
<name><surname>Winkler</surname>
<given-names>R.</given-names>
</name>
</person-group>
(<year>2006</year>
). <article-title>Multi-step methods for sdes and their application to problems with small noise</article-title>
. <source>SIAM J. Num. Anal.</source>
<volume>44</volume>
, <fpage>779</fpage>
–<lpage>803</lpage>
. <pub-id pub-id-type="doi">10.1137/040602857</pub-id>
</mixed-citation>
</ref>
<ref id="B8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Buckwar</surname>
<given-names>E.</given-names>
</name>
<name><surname>Winkler</surname>
<given-names>R.</given-names>
</name>
</person-group>
(<year>2007</year>
). <article-title>Multi-step maruyama methods for stochastic delay differential equations</article-title>
. <source>Stoch. Anal. Appl.</source>
<volume>25</volume>
, <fpage>933</fpage>
–<lpage>959</lpage>
. <pub-id pub-id-type="doi">10.1080/07362990701540311</pub-id>
</mixed-citation>
</ref>
<ref id="B9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Carletti</surname>
<given-names>M.</given-names>
</name>
</person-group>
(<year>2006</year>
). <article-title>Numerical solution of stochastic differential problems in the biosciences</article-title>
. <source>J. Comp. Appl. Math.</source>
<volume>185</volume>
, <fpage>422</fpage>
–<lpage>440</lpage>
. <pub-id pub-id-type="doi">10.1016/j.cam.2005.03.020</pub-id>
</mixed-citation>
</ref>
<ref id="B10"><mixed-citation publication-type="book"><person-group person-group-type="editor"><name><surname>Carnevale</surname>
<given-names>N.</given-names>
</name>
<name><surname>Hines</surname>
<given-names>M.</given-names>
</name>
</person-group>
(eds.). (<year>2006</year>
). <source>The NEURON Book</source>
. <publisher-loc>Cambridge, UK</publisher-loc>
: <publisher-name>Cambridge University Press</publisher-name>
.</mixed-citation>
</ref>
<ref id="B11"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
</person-group>
(<year>2005</year>
). <article-title>Waves, bumps and patterns in neural field theories</article-title>
. <source>Biol. Cybern.</source>
<volume>93</volume>
, <fpage>91</fpage>
–<lpage>108</lpage>
. <pub-id pub-id-type="doi">10.1007/s00422-005-0574-y</pub-id>
<pub-id pub-id-type="pmid">16059785</pub-id>
</mixed-citation>
</ref>
<ref id="B12"><mixed-citation publication-type="book"><person-group person-group-type="editor"><name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
<name><surname>beim Graben</surname>
<given-names>P.</given-names>
</name>
<name><surname>Potthast</surname>
<given-names>R.</given-names>
</name>
<name><surname>Wright</surname>
<given-names>J.</given-names>
</name>
</person-group>
(eds.). (<year>2014</year>
). <source>Neural Fields: Theory and Applications</source>
. <publisher-loc>New York, NY</publisher-loc>
: <publisher-name>Springer</publisher-name>
.</mixed-citation>
</ref>
<ref id="B13"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
<name><surname>Laing</surname>
<given-names>C.</given-names>
</name>
<name><surname>Schmidt</surname>
<given-names>H.</given-names>
</name>
<name><surname>Svanstedt</surname>
<given-names>N.</given-names>
</name>
<name><surname>Wyller</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<year>2012</year>
). <article-title>Waves in random neural media</article-title>
. <source>Disc. Cont. Dyn. Syst. A</source>
<volume>32</volume>
, <fpage>2951</fpage>
–<lpage>2970</lpage>
. <pub-id pub-id-type="doi">10.3934/dcds.2012.32.2951</pub-id>
</mixed-citation>
</ref>
<ref id="B14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
<name><surname>Owen</surname>
<given-names>M.</given-names>
</name>
</person-group>
(<year>2005</year>
). <article-title>Bumps, breathers, and waves in a neural network with spike frequency adaptation</article-title>
. <source>Phys. Rev. Lett.</source>
<volume>94</volume>
:<fpage>148102</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevLett.94.148102</pub-id>
<pub-id pub-id-type="pmid">15904116</pub-id>
</mixed-citation>
</ref>
<ref id="B15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
<name><surname>Venkov</surname>
<given-names>N.</given-names>
</name>
<name><surname>Shiau</surname>
<given-names>L.</given-names>
</name>
<name><surname>Bojak</surname>
<given-names>I.</given-names>
</name>
<name><surname>Liley</surname>
<given-names>D.</given-names>
</name>
<name><surname>Laing</surname>
<given-names>C.</given-names>
</name>
</person-group>
(<year>2007</year>
). <article-title>Modeling electrocortical activity through improved local approximations of integral neural field equations</article-title>
. <source>Phys. Rev. E</source>
<volume>76</volume>
:<fpage>051901</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.76.051901</pub-id>
<pub-id pub-id-type="pmid">18233681</pub-id>
</mixed-citation>
</ref>
<ref id="B16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deco</surname>
<given-names>G.</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>V. K.</given-names>
</name>
<name><surname>Robinson</surname>
<given-names>P.</given-names>
</name>
<name><surname>Breakspear</surname>
<given-names>M.</given-names>
</name>
<name><surname>Friston</surname>
<given-names>K.</given-names>
</name>
</person-group>
(<year>2008</year>
). <article-title>The dynamic brain: from spiking neurons to neural masses and cortical fields</article-title>
. <source>PLoS Comput. Biol.</source>
<volume>4</volume>
:<fpage>e1000092</fpage>
. <pub-id pub-id-type="doi">10.1371/journal.pcbi.1000092</pub-id>
<pub-id pub-id-type="pmid">18769680</pub-id>
</mixed-citation>
</ref>
<ref id="B17"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Demeulemeester</surname>
<given-names>H.</given-names>
</name>
<name><surname>Vandesande</surname>
<given-names>F.</given-names>
</name>
<name><surname>Orban</surname>
<given-names>G.</given-names>
</name>
<name><surname>Brandon</surname>
<given-names>C.</given-names>
</name>
<name><surname>Vanderhaeghen</surname>
<given-names>J. J.</given-names>
</name>
</person-group>
(<year>1988</year>
). <article-title>Heterogeneity of gabaergic cells in cat visual cortex</article-title>
. <source>J. Neurosci.</source>
<volume>8</volume>
, <fpage>988</fpage>
–<lpage>1000</lpage>
. <pub-id pub-id-type="pmid">2894415</pub-id>
</mixed-citation>
</ref>
<ref id="B18"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Elvin</surname>
<given-names>A. J.</given-names>
</name>
<name><surname>Laing</surname>
<given-names>C. R.</given-names>
</name>
<name><surname>Roberts</surname>
<given-names>M. G.</given-names>
</name>
</person-group>
(<year>2009</year>
). <article-title>Transient turing patterns in a neural field model</article-title>
. <source>Phys. Rev. E</source>
<volume>79</volume>
:<fpage>011911</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.79.011911</pub-id>
<pub-id pub-id-type="pmid">19257073</pub-id>
</mixed-citation>
</ref>
<ref id="B19"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Faye</surname>
<given-names>G.</given-names>
</name>
<name><surname>Faugeras</surname>
<given-names>O.</given-names>
</name>
</person-group>
(<year>2010</year>
). <article-title>Some theoretical and numerical results for delayed neural field equations</article-title>
. <source>Phys. D</source>
<volume>239</volume>
, <fpage>561</fpage>
–<lpage>578</lpage>
. <pub-id pub-id-type="doi">10.1016/j.physd.2010.01.010</pub-id>
</mixed-citation>
</ref>
<ref id="B20"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Folias</surname>
<given-names>S. E.</given-names>
</name>
<name><surname>Bressloff</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
(<year>2005</year>
). <article-title>Breathers in two-dimensional neural media</article-title>
. <source>Phys. Rev. Lett.</source>
<volume>95</volume>
:<fpage>208107</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevLett.95.208107</pub-id>
<pub-id pub-id-type="pmid">16384107</pub-id>
</mixed-citation>
</ref>
<ref id="B21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Friston</surname>
<given-names>K. J.</given-names>
</name>
<name><surname>Kahan</surname>
<given-names>J.</given-names>
</name>
<name><surname>Biswal</surname>
<given-names>B.</given-names>
</name>
<name><surname>Razi</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>A DCM for resting state fMRI</article-title>
. <source>Neuroimage</source>
<volume>94</volume>
, <fpage>396</fpage>
–<lpage>407</lpage>
. <pub-id pub-id-type="doi">10.1016/j.neuroimage.2013.12.009</pub-id>
<pub-id pub-id-type="pmid">24345387</pub-id>
</mixed-citation>
</ref>
<ref id="B22"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hashemi</surname>
<given-names>M.</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Sleigh</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>Anesthetic action on extra-synaptic receptors: effects in neural population models of EEG activity</article-title>
. <source>Front. Syst. Neurosci.</source>
<volume>8</volume>
:<issue>232</issue>
. <pub-id pub-id-type="doi">10.3389/fnsys.2014.00232</pub-id>
<pub-id pub-id-type="pmid">25540612</pub-id>
</mixed-citation>
</ref>
<ref id="B23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2007</year>
). <article-title>Generalization of the reaction-diffusion, Swift-Hohenberg, and Kuramoto-Sivashinsky equations and effects of finite propagation speeds</article-title>
. <source>Phys. Rev. E</source>
<volume>75</volume>
:<fpage>026214</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.75.026214</pub-id>
<pub-id pub-id-type="pmid">17358412</pub-id>
</mixed-citation>
</ref>
<ref id="B24"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2009</year>
). <article-title>Oscillatory activity in excitable neural systems</article-title>
. <source>Contemp. Phys.</source>
<volume>51</volume>
, <fpage>3</fpage>
–<lpage>16</lpage>
. <pub-id pub-id-type="doi">10.1080/00107510903293710</pub-id>
</mixed-citation>
</ref>
<ref id="B25"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Atay</surname>
<given-names>F.</given-names>
</name>
</person-group>
(<year>2006</year>
). <article-title>Effects of distributed transmission speeds on propagating activity in neural populations</article-title>
. <source>Phys. Rev. E</source>
<volume>73</volume>
:<fpage>021906</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.73.021906</pub-id>
<pub-id pub-id-type="pmid">16605361</pub-id>
</mixed-citation>
</ref>
<ref id="B26"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Bestehorn</surname>
<given-names>M.</given-names>
</name>
<name><surname>Wennekers</surname>
<given-names>T.</given-names>
</name>
</person-group>
(<year>2003</year>
). <article-title>Pattern formation in intracortical neuronal fields</article-title>
. <source>Netw. Comput. Neural Syst.</source>
<volume>14</volume>
, <fpage>351</fpage>
–<lpage>368</lpage>
. <pub-id pub-id-type="doi">10.1088/0954-898X_14_2_310</pub-id>
<pub-id pub-id-type="pmid">12790189</pub-id>
</mixed-citation>
</ref>
<ref id="B27"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Buhry</surname>
<given-names>L.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>Study of GABAergic extra-synaptic tonic inhibition in single neurons and neural populations by traversing neural scales: application to propofol-induced anaesthesia</article-title>
. <source>J. Comput. Neurosci.</source>
<volume>37</volume>
, <fpage>417</fpage>
–<lpage>437</lpage>
. <pub-id pub-id-type="doi">10.1007/s10827-014-0512-x</pub-id>
<pub-id pub-id-type="pmid">24976146</pub-id>
</mixed-citation>
</ref>
<ref id="B28"><mixed-citation publication-type="webpage"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Lefebvre</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<year>in press</year>
). <article-title>Stochastic center manifold analysis in scalar nonlinear systems involving distributed delays and additive noise</article-title>
. <source>Markov Process. Rel. Fields</source>
. Available online at: <ext-link ext-link-type="uri" xlink:href="https://hal.inria.fr/hal-01006388">https://hal.inria.fr/hal-01006388</ext-link>
</mixed-citation>
</ref>
<ref id="B29"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Rougier</surname>
<given-names>N.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>Numerical simulation scheme of one- and two-dimensional neural fields involving space-dependent delays</article-title>
, in <source>Neural Field Theory</source>
, eds <person-group person-group-type="editor"><name><surname>beim Graben</surname>
<given-names>P.</given-names>
</name>
<name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
<name><surname>Potthast</surname>
<given-names>R.</given-names>
</name>
<name><surname>Wright</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<publisher-loc>Berlin</publisher-loc>
: <publisher-name>Springer-Verlag</publisher-name>
), <fpage>175</fpage>
–<lpage>183</lpage>
.</mixed-citation>
</ref>
<ref id="B30"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Rougier</surname>
<given-names>N.</given-names>
</name>
</person-group>
(<year>2010</year>
). <article-title>Activity spread and breathers induced by finite transmission speeds in two-dimensional neural fields</article-title>
. <source>Phys. Rev. E</source>
<volume>82</volume>
:<fpage>R055701</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.82.055701</pub-id>
<pub-id pub-id-type="pmid">21230538</pub-id>
</mixed-citation>
</ref>
<ref id="B31"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Sutherland</surname>
<given-names>C.</given-names>
</name>
<name><surname>Longtin</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2008</year>
). <article-title>Driving neural oscillations with correlated spatial input and topographic feedback</article-title>
. <source>Phys. Rev. E</source>
<volume>78</volume>
:<fpage>021911</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.78.021911</pub-id>
<pub-id pub-id-type="pmid">18850869</pub-id>
</mixed-citation>
</ref>
<ref id="B32"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Idiart</surname>
<given-names>M. A. P.</given-names>
</name>
<name><surname>Abbott</surname>
<given-names>L. F.</given-names>
</name>
</person-group>
(<year>1993</year>
). <article-title>Propagation of excitation in neural network models</article-title>
. <source>Netw. Comput. Neural Syst.</source>
<volume>4</volume>
, <fpage>285</fpage>
–<lpage>294</lpage>
. <pub-id pub-id-type="doi">10.1088/0954-898X_4_3_003</pub-id>
</mixed-citation>
</ref>
<ref id="B33"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jirsa</surname>
<given-names>V. K.</given-names>
</name>
<name><surname>Haken</surname>
<given-names>H.</given-names>
</name>
</person-group>
(<year>1996</year>
). <article-title>Field theory of electromagnetic brain activity</article-title>
. <source>Phys. Rev. Lett.</source>
<volume>77</volume>
, <fpage>960</fpage>
–<lpage>963</lpage>
. <pub-id pub-id-type="pmid">10062950</pub-id>
</mixed-citation>
</ref>
<ref id="B34"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jirsa</surname>
<given-names>V.</given-names>
</name>
<name><surname>Jantzen</surname>
<given-names>K.</given-names>
</name>
<name><surname>Fuchs</surname>
<given-names>A.</given-names>
</name>
<name><surname>Kelso</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<year>2002</year>
). <article-title>Spatiotemporal forward solution of the EEG and MEG using network modelling</article-title>
. <source>IEEE Trans. Med. Imag.</source>
<volume>21</volume>
, <fpage>493</fpage>
–<lpage>504</lpage>
. <pub-id pub-id-type="doi">10.1109/TMI.2002.1009385</pub-id>
</mixed-citation>
</ref>
<ref id="B35"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kilpatrick</surname>
<given-names>Z. P.</given-names>
</name>
<name><surname>Bressloff</surname>
<given-names>P. C.</given-names>
</name>
</person-group>
(<year>2010</year>
). <article-title>Binocular rivalry in a competitive neural network with synaptic depression</article-title>
. <source>SIAM J. Appl. Dyn. Syst.</source>
<volume>9</volume>
, <fpage>1303</fpage>
–<lpage>1347</lpage>
. <pub-id pub-id-type="doi">10.1137/100788872</pub-id>
</mixed-citation>
</ref>
<ref id="B36"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Laing</surname>
<given-names>C.</given-names>
</name>
</person-group>
(<year>2005</year>
). <article-title>Spiral waves in nonlocal equations</article-title>
. <source>SIAM J. Appl. Dyn. Syst.</source>
<volume>4</volume>
, <fpage>588</fpage>
–<lpage>606</lpage>
. <pub-id pub-id-type="doi">10.1137/040612890</pub-id>
</mixed-citation>
</ref>
<ref id="B37"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Langtangen</surname>
<given-names>H. P.</given-names>
</name>
</person-group>
(<year>2006</year>
). <article-title>Numerical computing in python</article-title>
, in <source>Python Scripting for Computational Science, Volume 3 of Texts in Computational Science and Engineering</source>
, eds <person-group person-group-type="editor"><name><surname>Barth</surname>
<given-names>T. J.</given-names>
</name>
<name><surname>Griebel</surname>
<given-names>M.</given-names>
</name>
<name><surname>Keyes</surname>
<given-names>D. E.</given-names>
</name>
<name><surname>Nieminen</surname>
<given-names>R. M.</given-names>
</name>
<name><surname>Roose</surname>
<given-names>D.</given-names>
</name>
<name><surname>Schlick</surname>
<given-names>T.</given-names>
</name>
</person-group>
(<publisher-loc>Berlin; Heidelberg</publisher-loc>
: <publisher-name>Springer</publisher-name>
), <fpage>131</fpage>
–<lpage>181</lpage>
.</mixed-citation>
</ref>
<ref id="B38"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Molaee-Ardekani</surname>
<given-names>B.</given-names>
</name>
<name><surname>Senhadji</surname>
<given-names>L.</given-names>
</name>
<name><surname>Shamsollahi</surname>
<given-names>M. B.</given-names>
</name>
<name><surname>Vosoughi-Vahdat</surname>
<given-names>B.</given-names>
</name>
<name><surname>Wodey</surname>
<given-names>E.</given-names>
</name>
</person-group>
(<year>2007</year>
). <article-title>Brain activity modeling in general anesthesia: enhancing local mean-field models using a slow adaptive firing rate</article-title>
. <source>Phys. Rev. E</source>
<volume>76</volume>
:<fpage>041911</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevE.76.041911</pub-id>
<pub-id pub-id-type="pmid">17995030</pub-id>
</mixed-citation>
</ref>
<ref id="B39"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nunez</surname>
<given-names>P.</given-names>
</name>
</person-group>
(<year>1974</year>
). <article-title>The brain wave equation: a model for the EEG</article-title>
. <source>Math. Biosci.</source>
<volume>21</volume>
, <fpage>279</fpage>
–<lpage>291</lpage>
.</mixed-citation>
</ref>
<ref id="B40"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nunez</surname>
<given-names>P.</given-names>
</name>
</person-group>
(<year>2000</year>
). <article-title>Toward a quantitative description of large-scale neocortical dynamic function and EEG</article-title>
. <source>Behav. Brain Sci.</source>
<volume>23</volume>
, <fpage>371</fpage>
–<lpage>437</lpage>
. <pub-id pub-id-type="doi">10.1017/S0140525X00003253</pub-id>
<pub-id pub-id-type="pmid">11301576</pub-id>
</mixed-citation>
</ref>
<ref id="B41"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Nunez</surname>
<given-names>P.</given-names>
</name>
<name><surname>Srinivasan</surname>
<given-names>R.</given-names>
</name>
</person-group>
(<year>2006</year>
). <source>Electric Fields of the Brain: The Neurophysics of EEG</source>
. <publisher-loc>New York, NY; Oxford</publisher-loc>
: <publisher-name>Oxford University Press</publisher-name>
.</mixed-citation>
</ref>
<ref id="B42"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Owen</surname>
<given-names>M. R.</given-names>
</name>
<name><surname>Laing</surname>
<given-names>C. R.</given-names>
</name>
<name><surname>Coombes</surname>
<given-names>S.</given-names>
</name>
</person-group>
(<year>2007</year>
). <article-title>Bumps and rings in a two-dimensional neural field: splitting and rotational instabilities</article-title>
. <source>New J. Phys.</source>
<volume>9</volume>
:<fpage>378</fpage>
<pub-id pub-id-type="doi">10.1088/1367-2630/9/10/378</pub-id>
</mixed-citation>
</ref>
<ref id="B43"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pinotsis</surname>
<given-names>D.</given-names>
</name>
<name><surname>Friston</surname>
<given-names>K.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>Extracting novel information from neuroimaging data using neural fields</article-title>
. <source>EPJ Nonlin. Biomed. Phys.</source>
<volume>2</volume>
, <fpage>5</fpage>
<pub-id pub-id-type="doi">10.1140/epjnbp18</pub-id>
</mixed-citation>
</ref>
<ref id="B44"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pinotsis</surname>
<given-names>D. A.</given-names>
</name>
<name><surname>Moran</surname>
<given-names>R. J, Friston, K. J.</given-names>
</name>
</person-group>
(<year>2012</year>
). <article-title>Dynamic causal modeling with neural fields</article-title>
. <source>Neuroimage</source>
<volume>59</volume>
, <fpage>1261</fpage>
–<lpage>1274</lpage>
. <pub-id pub-id-type="doi">10.1016/j.neuroimage.2011.08.020</pub-id>
<pub-id pub-id-type="pmid">21924363</pub-id>
</mixed-citation>
</ref>
<ref id="B45"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pinotsis</surname>
<given-names>D. A.</given-names>
</name>
<name><surname>Schwarzkopf</surname>
<given-names>D. S.</given-names>
</name>
<name><surname>Litvak</surname>
<given-names>V.</given-names>
</name>
<name><surname>Rees</surname>
<given-names>G.</given-names>
</name>
<name><surname>Barnes</surname>
<given-names>G.</given-names>
</name>
<name><surname>Friston</surname>
<given-names>K. J.</given-names>
</name>
</person-group>
(<year>2013</year>
). <article-title>Dynamic causal modelling of lateral interactions in the visual cortex</article-title>
. <source>Neuroimage</source>
<volume>66</volume>
, <fpage>563</fpage>
–<lpage>576</lpage>
. <pub-id pub-id-type="doi">10.1016/j.neuroimage.2012.10.078</pub-id>
<pub-id pub-id-type="pmid">23128079</pub-id>
</mixed-citation>
</ref>
<ref id="B46"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pinto</surname>
<given-names>D.</given-names>
</name>
<name><surname>Ermentrout</surname>
<given-names>G.</given-names>
</name>
</person-group>
(<year>2001</year>
). <article-title>Spatially structured activity in synaptically coupled neuronal networks: I. travelling fronts and pulses</article-title>
. <source>SIAM J. Appl. Math.</source>
<volume>62</volume>
, <fpage>206</fpage>
–<lpage>225</lpage>
. <pub-id pub-id-type="doi">10.1137/S0036139900346453</pub-id>
</mixed-citation>
</ref>
<ref id="B47"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Qubbaj</surname>
<given-names>M. R.</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>V. K.</given-names>
</name>
</person-group>
(<year>2007</year>
). <article-title>Neural field dynamics with heterogeneous connection topology</article-title>
. <source>Phys. Rev. Lett.</source>
<volume>98</volume>
:<fpage>238102</fpage>
. <pub-id pub-id-type="doi">10.1103/PhysRevLett.98.238102</pub-id>
<pub-id pub-id-type="pmid">17677938</pub-id>
</mixed-citation>
</ref>
<ref id="B48"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rosenkilde</surname>
<given-names>C.</given-names>
</name>
</person-group>
(<year>1979</year>
). <article-title>Functional heterogeneity of the prefrontal cortex in the monkey: a review</article-title>
. <source>Behav. Neural Biol.</source>
<volume>25</volume>
, <fpage>301</fpage>
–<lpage>345</lpage>
. <pub-id pub-id-type="pmid">111666</pub-id>
</mixed-citation>
</ref>
<ref id="B49"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rossant</surname>
<given-names>C.</given-names>
</name>
<name><surname>Harris</surname>
<given-names>K. D.</given-names>
</name>
</person-group>
(<year>2013</year>
). <article-title>Hardware-accelerated interactive data visualization for neuroscience in python</article-title>
. <source>Front. Neuroinform.</source>
<volume>7</volume>
:<issue>36</issue>
. <pub-id pub-id-type="doi">10.3389/fninf.2013.00036</pub-id>
<pub-id pub-id-type="pmid">24391582</pub-id>
</mixed-citation>
</ref>
<ref id="B50"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rougier</surname>
<given-names>N. P.</given-names>
</name>
<name><surname>Fix</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<year>2012</year>
). <article-title>DANA: distributed numerical and adaptive modelling framework</article-title>
. <source>Network</source>
<volume>23</volume>
, <fpage>237</fpage>
–<lpage>253</lpage>
. <pub-id pub-id-type="doi">10.3109/0954898X.2012.721573</pub-id>
<pub-id pub-id-type="pmid">22994650</pub-id>
</mixed-citation>
</ref>
<ref id="B51"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rougier</surname>
<given-names>N.</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>2011</year>
). <article-title>Synchronous and asynchronous evaluation of dynamic neural fields</article-title>
. <source>J. Diff. Equ. Appl.</source>
<volume>17</volume>
, <fpage>1119</fpage>
–<lpage>1133</lpage>
. <pub-id pub-id-type="doi">10.1080/10236190903051575</pub-id>
</mixed-citation>
</ref>
<ref id="B52"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sanz Leon</surname>
<given-names>P.</given-names>
</name>
<name><surname>Knock</surname>
<given-names>S. A.</given-names>
</name>
<name><surname>Spiegler</surname>
<given-names>A.</given-names>
</name>
<name><surname>Jirsa</surname>
<given-names>V. K.</given-names>
</name>
</person-group>
(<year>2015</year>
). <article-title>Mathematical framework for large-scale brain network modeling in the virtual brain</article-title>
. <source>Neuroimage</source>
<volume>111</volume>
, <fpage>385</fpage>
–<lpage>430</lpage>
. <pub-id pub-id-type="doi">10.1016/j.neuroimage.2015.01.002</pub-id>
<pub-id pub-id-type="pmid">25592995</pub-id>
</mixed-citation>
</ref>
<ref id="B53"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sanz Leon</surname>
<given-names>P.</given-names>
</name>
<name><surname>Knock</surname>
<given-names>S. A.</given-names>
</name>
<name><surname>Woodman</surname>
<given-names>M. M.</given-names>
</name>
<name><surname>Domide</surname>
<given-names>L.</given-names>
</name>
<name><surname>Mersmann</surname>
<given-names>J.</given-names>
</name>
<name><surname>McIntosh</surname>
<given-names>A.</given-names>
</name>
<etal></etal>
</person-group>
. (<year>2013</year>
). <article-title>The Virtual Brain: a simulator of primate brain network dynamics</article-title>
. <source>Front. Neuroinform.</source>
<volume>7</volume>
:<issue>10</issue>
. <pub-id pub-id-type="doi">10.3389/fninf.2013.00010</pub-id>
<pub-id pub-id-type="pmid">23781198</pub-id>
</mixed-citation>
</ref>
<ref id="B54"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schmidt</surname>
<given-names>H.</given-names>
</name>
<name><surname>Hutt</surname>
<given-names>A.</given-names>
</name>
<name><surname>Schimansky-Geier</surname>
<given-names>L.</given-names>
</name>
</person-group>
(<year>2009</year>
). <article-title>Wave fronts in inhomogeneous neural field models</article-title>
. <source>Phys. D</source>
<volume>238</volume>
, <fpage>1101</fpage>
–<lpage>1112</lpage>
. <pub-id pub-id-type="doi">10.1016/j.physd.2009.02.017</pub-id>
</mixed-citation>
</ref>
<ref id="B55"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Steyn-Ross</surname>
<given-names>M. L.</given-names>
</name>
<name><surname>Steyn-Ross</surname>
<given-names>D. A.</given-names>
</name>
<name><surname>Wilson</surname>
<given-names>M. T.</given-names>
</name>
<name><surname>Sleigh</surname>
<given-names>J. W.</given-names>
</name>
</person-group>
(<year>2010</year>
). <article-title>Cortical patterns and gamma genesis are modulated by reversal potentials and gap-junction diffusion</article-title>
, in <source>Modeling Phase Transitions in the Brain, Vol. 4 of Springer Series in Computational Neuroscience</source>
, eds <person-group person-group-type="editor"><name><surname>Steyn-Ross</surname>
<given-names>D. A.</given-names>
</name>
<name><surname>Steyn-Ross</surname>
<given-names>M.</given-names>
</name>
</person-group>
(<publisher-loc>New York, NY</publisher-loc>
: <publisher-name>Springer</publisher-name>
), <fpage>271</fpage>
–<lpage>299</lpage>
.</mixed-citation>
</ref>
<ref id="B56"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stimberg</surname>
<given-names>M.</given-names>
</name>
<name><surname>Goodman</surname>
<given-names>D. F.</given-names>
</name>
<name><surname>Benichoux</surname>
<given-names>V.</given-names>
</name>
<name><surname>Brette</surname>
<given-names>R.</given-names>
</name>
</person-group>
(<year>2014</year>
). <article-title>Equation-oriented specification of neural models for simulations</article-title>
. <source>Front. Neuroinform.</source>
<volume>8</volume>
:<issue>6</issue>
. <pub-id pub-id-type="doi">10.3389/fninf.2014.00006</pub-id>
<pub-id pub-id-type="pmid">24550820</pub-id>
</mixed-citation>
</ref>
<ref id="B57"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Turing</surname>
<given-names>A.</given-names>
</name>
</person-group>
(<year>1952</year>
). <article-title>The chemical basis of morphogenesis</article-title>
. <source>Philos. Trans. R. Soc. Lond. B Biol. Sci.</source>
<volume>237</volume>
, <fpage>37</fpage>
–<lpage>72</lpage>
.</mixed-citation>
</ref>
<ref id="B58"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Van Loan</surname>
<given-names>C.</given-names>
</name>
</person-group>
(<year>1991</year>
). <source>Computational Frameworks for the Fast Fourier Transform</source>
. <publisher-loc>Philadelphia, PA</publisher-loc>
: <publisher-name>SIAM</publisher-name>
.</mixed-citation>
</ref>
<ref id="B59"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Veltz</surname>
<given-names>R.</given-names>
</name>
<name><surname>Faugeras</surname>
<given-names>O.</given-names>
</name>
</person-group>
(<year>2011</year>
). <article-title>Stability of the stationary solutions of neural field equations with propagation delays</article-title>
. <source>J. Math. Neurosci.</source>
<volume>1</volume>
:<fpage>1</fpage>
. <pub-id pub-id-type="doi">10.1186/2190-8567-1-1</pub-id>
<pub-id pub-id-type="pmid">22655751</pub-id>
</mixed-citation>
</ref>
<ref id="B60"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Veltz</surname>
<given-names>R.</given-names>
</name>
<name><surname>Faugeras</surname>
<given-names>O.</given-names>
</name>
</person-group>
(<year>2013</year>
). <article-title>A center manifold result for delayed neural fields equations</article-title>
. <source>SIAM J. Math. Anal.</source>
<volume>45</volume>
, <fpage>1527</fpage>
–<lpage>1562</lpage>
. <pub-id pub-id-type="doi">10.1137/110856162</pub-id>
</mixed-citation>
</ref>
<ref id="B61"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>X.-J.</given-names>
</name>
</person-group>
(<year>2010</year>
). <article-title>Neurophysiological and computational principles of cortical rhythms in cognition</article-title>
. <source>Physiol. Rev.</source>
<volume>90</volume>
, <fpage>1195</fpage>
–<lpage>1268</lpage>
. <pub-id pub-id-type="doi">10.1152/physrev.00035.2008</pub-id>
<pub-id pub-id-type="pmid">20664082</pub-id>
</mixed-citation>
</ref>
<ref id="B62"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name><surname>Markram</surname>
<given-names>H.</given-names>
</name>
<name><surname>Goodman</surname>
<given-names>P. H.</given-names>
</name>
<name><surname>Berger</surname>
<given-names>T. K.</given-names>
</name>
<name><surname>Ma</surname>
<given-names>J.</given-names>
</name>
<name><surname>Goldman-Rakic</surname>
<given-names>P. S.</given-names>
</name>
</person-group>
(<year>2006</year>
). <article-title>Heterogeneity in the pyramidal network of the medial prefrontal cortex</article-title>
. <source>Nat. Neurosci.</source>
<volume>9</volume>
, <fpage>534</fpage>
–<lpage>542</lpage>
. <pub-id pub-id-type="doi">10.1038/nn1670</pub-id>
<pub-id pub-id-type="pmid">16547512</pub-id>
</mixed-citation>
</ref>
<ref id="B63"><mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Wertheimer</surname>
<given-names>M.</given-names>
</name>
<name><surname>Spillmann</surname>
<given-names>L.</given-names>
</name>
<name><surname>Sarris</surname>
<given-names>V.</given-names>
</name>
<name><surname>Sekuler</surname>
<given-names>R.</given-names>
</name>
</person-group>
(<year>2012</year>
). <source>On Perceived Motion and Figural Organization</source>
. <publisher-loc>Cambridge, MA</publisher-loc>
: <publisher-name>MITPress</publisher-name>
.</mixed-citation>
</ref>
<ref id="B64"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wilson</surname>
<given-names>H. R.</given-names>
</name>
<name><surname>Cowan</surname>
<given-names>J.</given-names>
</name>
</person-group>
(<year>1973</year>
). <article-title>A mathematical theory of the functional dynamics of cortical and thalamic nervous tissue</article-title>
. <source>Kybernetik</source>
<volume>13</volume>
, <fpage>55</fpage>
–<lpage>80</lpage>
. <pub-id pub-id-type="pmid">4767470</pub-id>
</mixed-citation>
</ref>
<ref id="B65"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wright</surname>
<given-names>J. J.</given-names>
</name>
<name><surname>Kydd</surname>
<given-names>R.</given-names>
</name>
</person-group>
(<year>1992</year>
). <article-title>The electroencephalogram and cortical neural networks</article-title>
. <source>Network</source>
<volume>3</volume>
, <fpage>341</fpage>
–<lpage>362</lpage>
. <pub-id pub-id-type="doi">10.1088/0954-898X_3_3_006</pub-id>
</mixed-citation>
</ref>
<ref id="B66"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wright</surname>
<given-names>J. J. Liley, D. T. J.</given-names>
</name>
</person-group>
(<year>1994</year>
). <article-title>A millimetric-scale simulation of electrocortical wave dynamics based on anatomical estimates of cortical synaptic density</article-title>
. <source>Netw. Comput. Neural Syst.</source>
<volume>5</volume>
, <fpage>191</fpage>
–<lpage>202</lpage>
. <pub-id pub-id-type="doi">10.1088/0954-898X_5_2_005</pub-id>
</mixed-citation>
</ref>
</ref-list>
<app-group><app id="A1"><title>Appendix i</title>
<sec><title>Heterogeneous neural fields</title>
<p>Heterogeneous neural fields have attracted increased attention in recent years (Bressloff, <xref rid="B5" ref-type="bibr">2001</xref>
; Qubbaj and Jirsa, <xref rid="B47" ref-type="bibr">2007</xref>
; Brackley and Turner, <xref rid="B3" ref-type="bibr">2009</xref>
; Schmidt et al., <xref rid="B54" ref-type="bibr">2009</xref>
; Coombes et al., <xref rid="B13" ref-type="bibr">2012</xref>
; beim Graben and Hutt, <xref rid="B2" ref-type="bibr">2014</xref>
) since they have been found in biological neural networks such as the prefrontal cortex (Rosenkilde, <xref rid="B48" ref-type="bibr">1979</xref>
; Wang et al., <xref rid="B62" ref-type="bibr">2006</xref>
) and visual cortex (Demeulemeester et al., <xref rid="B17" ref-type="bibr">1988</xref>
). In order to study the neural population activity in such systems, the present implementation could be extended along the following mathematical reasoning. The integral in Equation (1) may be re-written as</p>
<disp-formula id="E15"><mml:math id="M27"><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:mrow><mml:msub><mml:mo>∫</mml:mo>
<mml:mi>Ω</mml:mi>
</mml:msub>
<mml:mrow><mml:msup><mml:mtext>d</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi>y</mml:mi>
<mml:mtext> </mml:mtext>
<mml:mi>K</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mi>S</mml:mi>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mi>V</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>−</mml:mo>
<mml:mfrac><mml:mrow><mml:mo>∥</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>∥</mml:mo>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:mrow><mml:msub><mml:mo>∫</mml:mo>
<mml:mi>Ω</mml:mi>
</mml:msub>
<mml:mrow><mml:msup><mml:mtext>d</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi>y</mml:mi>
<mml:mstyle displaystyle="true"><mml:mrow><mml:msubsup><mml:mo>∫</mml:mo>
<mml:mrow><mml:mo>−</mml:mo>
<mml:mi>∞</mml:mi>
</mml:mrow>
<mml:mi>∞</mml:mi>
</mml:msubsup>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:mstyle>
<mml:mi>τ</mml:mi>
<mml:mi>K</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo stretchy="false">)</mml:mo>
<mml:mi>δ</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>τ</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac><mml:mrow><mml:mo>∥</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>∥</mml:mo>
</mml:mrow>
<mml:mi>c</mml:mi>
</mml:mfrac>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mi>S</mml:mi>
<mml:mrow><mml:mo>[</mml:mo>
<mml:mrow><mml:mi>V</mml:mi>
<mml:mrow><mml:mo>(</mml:mo>
<mml:mrow><mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>τ</mml:mi>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr><mml:mtd><mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true"><mml:mrow><mml:msub><mml:mo>∫</mml:mo>
<mml:mi>Ω</mml:mi>
</mml:msub>
<mml:mrow><mml:msup><mml:mtext>d</mml:mtext>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mstyle>
<mml:mi>y</mml:mi>
<mml:mstyle displaystyle="true"><mml:mrow><mml:msubsup><mml:mo>∫</mml:mo>
<mml:mn>0</mml:mn>
<mml:mrow><mml:msub><mml:mi>τ</mml:mi>
<mml:mrow><mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:mstyle>
<mml:mi>T</mml:mi>
<mml:mi>D</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>−</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
<mml:mtext> </mml:mtext>
<mml:mi>R</mml:mi>
<mml:mo stretchy="false">(</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>x</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mstyle mathvariant="bold" mathsize="normal"><mml:mi>y</mml:mi>
</mml:mstyle>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>−</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>with <italic>D</italic>
(<bold>x</bold>
, <italic>t</italic>
) = δ(∥ <bold>x</bold>
∥ ∕<italic>c</italic>
− <italic>t</italic>
), <italic>R</italic>
(<bold>x, y</bold>
, <italic>t</italic>
) = <italic>K</italic>
(<bold>x, y</bold>
)<italic>S</italic>
[<italic>V</italic>
(<bold>y</bold>
, <italic>t</italic>
)]. This integral still has a spatial convolution structure; however, it is not perfect since <italic>R</italic>
includes the spatial location <bold>x</bold>
. The numerical simulation of the integral consequently involves <inline-formula><mml:math id="M28"><mml:msup><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mn>4</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msub><mml:mrow><mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow><mml:mi>r</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
summands and the numerical implementation is slower than in the homogeneous case. The formulation of heterogeneous neural fields is not implemented yet, but it will be part of a future update.</p>
</sec>
</app>
</app-group>
<app-group><app id="A2"><title>Appendix ii</title>
<table-wrap id="TA1" position="anchor"><label>Table A1</label>
<caption><p><bold>Keyboard keys and their actions</bold>
.</p>
</caption>
<table frame="hsides" rules="groups"><thead><tr><th valign="top" align="left" rowspan="1" colspan="1"><bold>Key</bold>
</th>
<th valign="top" align="left" rowspan="1" colspan="1"><bold>Action</bold>
</th>
</tr>
</thead>
<tbody><tr style="border-top: thin solid #000000;"><td valign="top" align="left" rowspan="1" colspan="1">2</td>
<td valign="top" align="left" rowspan="1" colspan="1">Interpolate between min and max graph value colors</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">3</td>
<td valign="top" align="left" rowspan="1" colspan="1">Interpolate between min-mid and mid-max graph value colors</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">↑</td>
<td valign="top" align="left" rowspan="1" colspan="1">Rotate the field up</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">↓</td>
<td valign="top" align="left" rowspan="1" colspan="1">Rotate the field down</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">←</td>
<td valign="top" align="left" rowspan="1" colspan="1">Rotate the field left</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">→</td>
<td valign="top" align="left" rowspan="1" colspan="1">Rotate the field right</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">a</td>
<td valign="top" align="left" rowspan="1" colspan="1">Modulate minimum graph value color</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">b</td>
<td valign="top" align="left" rowspan="1" colspan="1">Modulate background graph color</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">d</td>
<td valign="top" align="left" rowspan="1" colspan="1">Move the field down</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">e</td>
<td valign="top" align="left" rowspan="1" colspan="1">Move the field up</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">Esc</td>
<td valign="top" align="left" rowspan="1" colspan="1">Exit simulation</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">f</td>
<td valign="top" align="left" rowspan="1" colspan="1">Move the field right</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">g</td>
<td valign="top" align="left" rowspan="1" colspan="1">Change number of axes lines</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">i</td>
<td valign="top" align="left" rowspan="1" colspan="1">Save an image</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">j</td>
<td valign="top" align="left" rowspan="1" colspan="1">Set min and max z axis limits to min and max field values</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">k</td>
<td valign="top" align="left" rowspan="1" colspan="1">Change color distribution to a higher range</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">l</td>
<td valign="top" align="left" rowspan="1" colspan="1">Change color distribution to a reduced range</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">m</td>
<td valign="top" align="left" rowspan="1" colspan="1">Set minimum z axis limit to minimum field value</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">n</td>
<td valign="top" align="left" rowspan="1" colspan="1">Change minimum z axis limit</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">o</td>
<td valign="top" align="left" rowspan="1" colspan="1">Equally distribute color range</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">p</td>
<td valign="top" align="left" rowspan="1" colspan="1">Pause and resume simulation</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1"><italic>pg up</italic>
</td>
<td valign="top" align="left" rowspan="1" colspan="1">Zoom in</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1"><italic>pg down</italic>
</td>
<td valign="top" align="left" rowspan="1" colspan="1">Zoom out</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">q</td>
<td valign="top" align="left" rowspan="1" colspan="1">Modulate middle graph value color</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">s</td>
<td valign="top" align="left" rowspan="1" colspan="1">Move the field left</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">t</td>
<td valign="top" align="left" rowspan="1" colspan="1">Change axis text size on Mac systems</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">u</td>
<td valign="top" align="left" rowspan="1" colspan="1">Set maximum z axis limit to maximum field value</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">v</td>
<td valign="top" align="left" rowspan="1" colspan="1">Begin and end video recording</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">y</td>
<td valign="top" align="left" rowspan="1" colspan="1">Change maximum z axis limit</td>
</tr>
<tr><td valign="top" align="left" rowspan="1" colspan="1">z</td>
<td valign="top" align="left" rowspan="1" colspan="1">Modulate maximum graph value color</td>
</tr>
</tbody>
</table>
</table-wrap>
</app>
</app-group>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Lorraine/explor/InforLorV4/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000004 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000004 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Lorraine |area= InforLorV4 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:4611063 |texte= Neural field simulator: two-dimensional spatio-temporal dynamics involving finite transmission speed }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:26539105" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a InforLorV4
![]() | This area was generated with Dilib version V0.6.33. | ![]() |