cell_hex27.c

一个用来实现偏微分方程中网格的计算库

	    conn[5] = this->node(21)+1;
	    conn[6] = this->node(26)+1;
	    conn[7] = this->node(24)+1;

	    return;
      
	  case 1:

	    conn[0] = this->node(8)+1;
	    conn[1] = this->node(1)+1;
	    conn[2] = this->node(9)+1;
	    conn[3] = this->node(20)+1;
	    conn[4] = this->node(21)+1;
	    conn[5] = this->node(13)+1;
	    conn[6] = this->node(22)+1;
	    conn[7] = this->node(26)+1;

	    return;
      
	  case 2:

	    conn[0] = this->node(11)+1;
	    conn[1] = this->node(20)+1;
	    conn[2] = this->node(10)+1;
	    conn[3] = this->node(3)+1; 
	    conn[4] = this->node(24)+1;
	    conn[5] = this->node(26)+1;
	    conn[6] = this->node(23)+1;
	    conn[7] = this->node(15)+1;

	    return;
      
	  case 3:

	    conn[0] = this->node(20)+1;
	    conn[1] = this->node(9)+1;
	    conn[2] = this->node(2)+1;
	    conn[3] = this->node(10)+1;
	    conn[4] = this->node(26)+1;
	    conn[5] = this->node(22)+1;
	    conn[6] = this->node(14)+1;
	    conn[7] = this->node(23)+1;

	    return;
      
	  case 4:

	    conn[0] = this->node(12)+1;
	    conn[1] = this->node(21)+1;
	    conn[2] = this->node(26)+1;
	    conn[3] = this->node(24)+1;
	    conn[4] = this->node(4)+1;
	    conn[5] = this->node(16)+1;
	    conn[6] = this->node(25)+1;
	    conn[7] = this->node(19)+1;

	    return;
      
	  case 5:

	    conn[0] = this->node(21)+1;
	    conn[1] = this->node(13)+1;
	    conn[2] = this->node(22)+1;
	    conn[3] = this->node(26)+1;
	    conn[4] = this->node(16)+1;
	    conn[5] = this->node(5)+1;
	    conn[6] = this->node(17)+1;
	    conn[7] = this->node(25)+1;

	    return;
      
	  case 6:

	    conn[0] = this->node(24)+1;
	    conn[1] = this->node(26)+1;
	    conn[2] = this->node(23)+1;
	    conn[3] = this->node(15)+1;
	    conn[4] = this->node(19)+1;
	    conn[5] = this->node(25)+1;
	    conn[6] = this->node(18)+1;
	    conn[7] = this->node(7)+1; 

	    return;
      
	  case 7:
      
	    conn[0] = this->node(26)+1;
	    conn[1] = this->node(22)+1;
	    conn[2] = this->node(14)+1;
	    conn[3] = this->node(23)+1;
	    conn[4] = this->node(25)+1;
	    conn[5] = this->node(17)+1;
	    conn[6] = this->node(6)+1;
	    conn[7] = this->node(18)+1;

	    return;

	  default:
	    libmesh_error();
	  }
      }
      
    case VTK:
      {
	switch (sc)
	  {
	  case 0:

	    conn[0] = this->node(0);
	    conn[1] = this->node(8);
	    conn[2] = this->node(20);
	    conn[3] = this->node(11);
	    conn[4] = this->node(12);
	    conn[5] = this->node(21);
	    conn[6] = this->node(26);
	    conn[7] = this->node(24);

	    return;
      
	  case 1:

	    conn[0] = this->node(8);
	    conn[1] = this->node(1);
	    conn[2] = this->node(9);
	    conn[3] = this->node(20);
	    conn[4] = this->node(21);
	    conn[5] = this->node(13);
	    conn[6] = this->node(22);
	    conn[7] = this->node(26);

	    return;
      
	  case 2:

	    conn[0] = this->node(11);
	    conn[1] = this->node(20);
	    conn[2] = this->node(10);
	    conn[3] = this->node(3); 
	    conn[4] = this->node(24);
	    conn[5] = this->node(26);
	    conn[6] = this->node(23);
	    conn[7] = this->node(15);

	    return;
      
	  case 3:

	    conn[0] = this->node(20);
	    conn[1] = this->node(9);
	    conn[2] = this->node(2);
	    conn[3] = this->node(10);
	    conn[4] = this->node(26);
	    conn[5] = this->node(22);
	    conn[6] = this->node(14);
	    conn[7] = this->node(23);

	    return;
      
	  case 4:

	    conn[0] = this->node(12);
	    conn[1] = this->node(21);
	    conn[2] = this->node(26);
	    conn[3] = this->node(24);
	    conn[4] = this->node(4);
	    conn[5] = this->node(16);
	    conn[6] = this->node(25);
	    conn[7] = this->node(19);

	    return;
      
	  case 5:

	    conn[0] = this->node(21);
	    conn[1] = this->node(13);