stream.c

基于格子Boltzmann方法开源可视化软件的源代码 具有很高的实用价值。对学习LBM方法及其软件开发非常游泳

//##############################################################################
//
// stream.c
//

void stream( lattice_ptr lattice)
{
  pdf_ptr pdf;

  double  *f;

  int     *nn;
  
  int     n, a;

#if SAY_HI
  printf("stream() -- Hi!\n");
#endif /* SAY_HI */

  // Get pointer to pdf.
  pdf = lattice->pdf;

  // Get pointer to distribution functions.
  f = pdf->f;

  // Get pointer to neighbors.
  nn = lattice->node->nn;

  // Traverse nodes.
  for( n=0; n<lattice->NumNodes; n++)
  {
    // Copy the rest particle f[0] to ftemp[0].
    *(f+9) = *f++;

    for( a=1; a<=8; a++)
    {
      // S T R E A M   
      //
      //  - f[a] streams to ftemp[a] of its neighboring node.
      //
      if( *nn != INACTIVE_NODE) 
      {
        pdf[ *nn].ftemp[a] = *f;

//printf("STREAM: %d from %d to %d.\n", a, n, *nn);

      } /* if( *nn != INACTIVE_NODE)  */
      else
      {
//printf("STREAM: Skipping %d at %d.\n", a, n);
      }

      // Move nn to neighbor of next element of f.
      nn++;

      // Move next element of f.
      f++;

    } /* for( a=1; a<=8; a++) */

    // Advance nn to neighbor of first element of next node.
    nn+=3;

    // Advance to the next node.
    f+=18; 

  } /* for( n=0; n<NumNodes; n++) */

#if SAY_HI
  printf("stream() -- Bye!\n");
#endif /* SAY_HI */

} /* void stream( lattice_ptr lattice) */