geom_fill()

geom_fill

(

struct geometry *

geom

)

Fills the geometry structure with points coordinates, topology and vectors, as produced by the respective builders according to the input data.

Parameters

[in,out]

geom

                                     {
  DOUBLE pi;
  int inodb,i,j;
  int ielmb,ivelb;

  struct vector vtmp1,vtmp2;
  struct vector a1c,a2c,a1xa2c;
  DOUBLE sx,sy,sz;

  logger(LOG_DEBUG,"filling geometry arrays\n");
  pi = 4.0 * atanl(1.0);

  // Translating and rotating nodes
  sx = geom->scale.x != 0.0 ? geom->scale.x : 1.0;
  sy = geom->scale.y != 0.0 ? geom->scale.y : 1.0;
  sz = geom->scale.z != 0.0 ? geom->scale.z : 1.0;
  logger(LOG_DEBUG,"X Scale factor:  %f\n",sx);
  logger(LOG_DEBUG,"Y Scale factor:  %f\n",sy);
  logger(LOG_DEBUG,"Z Scale factor:  %f\n",sz);

  logger(LOG_DEBUG,"filling geometry arrays\n");

  for(inodb=0;inodb<geom->nnodb;inodb++){
    geometry->nodes[inodb + _g_inodb].x *= sx;
    geometry->nodes[inodb + _g_inodb].y *= sy;
    geometry->nodes[inodb + _g_inodb].z *= sz;

    vec_rotate(&geometry->nodes[inodb + _g_inodb], geom->rotation.x*pi/180.0, geom->rotation.y*pi/180.0, geom->rotation.z*pi/180.0);

    geometry->nodes[inodb + _g_inodb].x += geom->translation.x;
    geometry->nodes[inodb + _g_inodb].y += geom->translation.y;
    geometry->nodes[inodb + _g_inodb].z += geom->translation.z;
  }


  // topology 
  for(i=0;i<geom->nelmb;i++){
    geometry->jphib[i + _g_ielmb] = i;
    geometry->jpsib[i + _g_ielmb] = i;
    geometry->mvrtb[i + _g_ielmb] = 4;
    geometry->kcreb[i + _g_ielmb] = 1;
    geometry->kelmb[i + _g_ielmb] = 0;
  }

  /* building elements */
  for(i=0;i<geom->nelmb;i++){
    for(j=0;j<4;j++){
      ielmb = i;
      ivelb = j21(j,ielmb,4);
      inodb = geometry->jnodb[ivelb + _g_ielmb];
      vec_copy(&geometry->elements[i + _g_ielmb].x[j], geometry->nodes[inodb + _g_inodb]);
      vec_mul(&vtmp1,geometry->nodes[inodb + _g_inodb],.25);
      geometry->elements[i + _g_ielmb].idx[j]=inodb + _g_inodb;
    }
  }


  /* Building Control points  */
  for(i=0;i<geom->ncntr;i++){
    vec_zero(&geometry->cntr[i + _g_icntr]);
    for(j=0;j<4;j++){
      ielmb = i;
      ivelb = j21(j,ielmb,4);
      inodb = geometry->jnodb[ivelb + _g_ielmb];

      vec_mul(&vtmp1,geometry->nodes[inodb + _g_inodb],.25);
      vec_sum(&geometry->cntr[i + _g_icntr], geometry->cntr[i + _g_icntr], vtmp1);
    }
  }


  // normals
  for(i=_g_ielmb;i<(geom->nelmb + _g_ielmb);i++){

    vec_diff(&vtmp1,geometry->elements[i].x[0],geometry->elements[i].x[2]);
    vec_diff(&vtmp2,geometry->elements[i].x[1],geometry->elements[i].x[3]);

    a1c.x = .25*(+geometry->elements[i].x[0].x + geometry->elements[i].x[1].x - geometry->elements[i].x[2].x - geometry->elements[i].x[3].x);
    a1c.y = .25*(+geometry->elements[i].x[0].y + geometry->elements[i].x[1].y - geometry->elements[i].x[2].y - geometry->elements[i].x[3].y);
    a1c.z = .25*(+geometry->elements[i].x[0].z + geometry->elements[i].x[1].z - geometry->elements[i].x[2].z - geometry->elements[i].x[3].z);

    a2c.x = .25*(-geometry->elements[i].x[0].x + geometry->elements[i].x[1].x + geometry->elements[i].x[2].x - geometry->elements[i].x[3].x);
    a2c.y = .25*(-geometry->elements[i].x[0].y + geometry->elements[i].x[1].y + geometry->elements[i].x[2].y - geometry->elements[i].x[3].y);
    a2c.z = .25*(-geometry->elements[i].x[0].z + geometry->elements[i].x[1].z + geometry->elements[i].x[2].z - geometry->elements[i].x[3].z);

    vec_cross(&a1xa2c,a1c,a2c);

    vec_copy (&geometry->elements[i].n,a1xa2c);
    geometry->elements[i].nmod = vec_mod(geometry->elements[i].n); /* storing modules */
    vec_normalize(&geometry->elements[i].n);

  }

   geometry->ncntrc += geom->ncntr ;
   for(i=_g_icntr;i<(geom->ncntr+_g_icntr);i++){
       vec_copy(&geometry->cntrc[i],geometry->cntr[i]);
   }
  logger(LOG_DEBUG,"geometry arrays filled\n");
}