Research Article

Computational Modeling of Microabscess Formation

Algorithm 1

float chemotaxis (float** , float **ch, int , int , float ** )
{
 flux_left = 0;
 flux_right = 0;
 if ( > 0)
 {
    if (( )
    {
     flux_left = − (ch – ch[ ])* [ ] * ( )* )/delta ;
   }
   else
   {
     flux_left = − (ch[ ] – ch[ ])* * ( )* ( )/delta ;
   }
 }
 if ( < (size − 1))
 {
   if ((ch[ – ch[ ]) > 0)
   {
     flux_right = (ch[ ] [ ] – ch[ ])* [ ] [ ]* ( )* ( )/delta ;
   }
   else
   {
     flux_right = (ch[ ] [ ] – ch[ ])* [ ]* ( ])* ( )/delta ;
   }
 }
 res = (flux_left + flux_right)/delta ;
 flux_up = 0;
 flux_down = 0;
 if ( > 0)
 {
   if ((ch[ ] [ ] – ch[ ] [ ]) > 0)
   {
     flux_up = − (ch[ ] [ ] – ch[ ][ ])* [ ][ ]* ( ][ ])* ( )/delta ;
   }
   else
   {
     flux_up = − (ch[ ] [ ] – ch[ ][ ])* [ [ ]* ( )* ( )/delta ;
   }
 }
 if ( < (size − 1))
 {
   if ((ch[ ] [ ] – ch[ ] [ ]) > 0)
   {
     flux_down = (ch[ ] [ ] – ch[ ] [ ])* ] [ ]* )* ] [ ])/delta ;
   }
   else
   {
     flux_down = (ch[ ] [ ] – ch[ ] [ ])* [ ] [ ]* ( ] [ ])* /delta ;
   }
 }
 res = (flux_up + flux_down)/delta ;
 return (res + res );
}