exoselect.c File Reference

#include <stdlib.h>
#include <math.h>
#include <stdio.h>
#include <string.h>
#include "CalculiX.h"
#include "exodusII.h"
#include "exo.h"

Include dependency graph for exoselect.c:

Go to the source code of this file.

Macros
#define	min(a, b)   ((a) <= (b) ? (a) : (b))
#define	max(a, b)   ((a) >= (b) ? (a) : (b))

Functions
void	exoselect (double *field1, double *field2, ITG *iset, ITG *nkcoords, ITG *inum, ITG *istartset, ITG *iendset, ITG *ialset, ITG *ngraph, ITG *ncomp, ITG *ifield, ITG *icomp, ITG *nfield, ITG *iselect, ITG exoid, ITG time_step, int countvar, ITG nout)

Macro Definition Documentation

max

#define max	(		a,
			b
	)		((a) >= (b) ? (a) : (b))

min

#define min	(		a,
			b
	)		((a) <= (b) ? (a) : (b))

Function Documentation

exoselect()

void exoselect	(	double *	field1,
		double *	field2,
		ITG *	iset,
		ITG *	nkcoords,
		ITG *	inum,
		ITG *	istartset,
		ITG *	iendset,
		ITG *	ialset,
		ITG *	ngraph,
		ITG *	ncomp,
		ITG *	ifield,
		ITG *	icomp,
		ITG *	nfield,
		ITG *	iselect,
		ITG	exoid,
		ITG	time_step,
		int	countvar,
		ITG	nout
	)

                                                 {

  /* storing scalars, components of vectors and tensors without additional
     transformations */

  /* number of components in field1: nfield[0]
     number of components in field2: nfield[1]

     number of entities to store: ncomp
     for each entity i, 0<=i<ncomp:
     - ifield[i]: 1=field1,2=field2
     - icomp[i]: component: 0...,(nfield[0]-1 or nfield[1]-1) */

  ITG i,j,k,l,m,o,nksegment;
  /* When initializing parameter values:
     "g" (or "G")
     For global variables.
     "n" (or "N")
     For nodal variables.
     "e" (or "E")
     For element variables.
     "m" (or "M")
     For nodeset variables.
     "s" (or "S")
     For sideset variables.
  */
  // ITG num_nod_vars = *ncomp;
  float *nodal_var_vals;
  nodal_var_vals = (float *) calloc (nout, sizeof(float));


  for(j=0;j<*ncomp;j++){
    if(*iset==0){
      m=0;
      for(i=0;i<*nkcoords;i++){
    /* check whether output is requested for solid nodes or
       network nodes */

    if(*iselect==1){
      if(inum[i]<=0) continue;
    }else if(*iselect==-1){
      if(inum[i]>=0) continue;
    }else{
      if(inum[i]==0) continue;
    }

    /* storing the entities */
    if(ifield[j]==1){
      nodal_var_vals[m]=field1[i*nfield[0]+icomp[j]];
    }else{
      nodal_var_vals[m]=field2[i*nfield[1]+icomp[j]];
    }
    m++;
      }
    }else{
      nksegment=(*nkcoords)/(*ngraph);
      m=0;
      for(k=istartset[*iset-1]-1;k<iendset[*iset-1];k++){
    if(ialset[k]>0){
    
      for(l=0;l<*ngraph;l++){
        i=ialset[k]+l*nksegment-1;
    
        /* check whether output is requested for solid nodes or
           network nodes */
    
        if(*iselect==1){
          if(inum[i]<=0) continue;
        }else if(*iselect==-1){
          if(inum[i]>=0) continue;
        }else{
          if(inum[i]==0) continue;
        }
    
        /* storing the entities */
        if(ifield[j]==1){
          nodal_var_vals[m]=field1[i*nfield[0]+icomp[j]];
        }else{
          nodal_var_vals[m]=field2[i*nfield[1]+icomp[j]];
        }
        m++;
      }
    }else{
      l=ialset[k-2];
      do{
        l-=ialset[k];
        if(l>=ialset[k-1]) break;
        o=0;
        for(m=0;m<*ngraph;m++){
          i=l+m*nksegment-1;
    
          /* check whether output is requested for solid nodes or
         network nodes */
    
          if(*iselect==1){
        if(inum[i]<=0) continue;
          }else if(*iselect==-1){
        if(inum[i]>=0) continue;
          }else{
        if(inum[i]==0) continue;
          }
    
          /* storing the entities */
          if(ifield[j]==1){
        nodal_var_vals[o]=field1[i*nfield[0]+icomp[j]];
          }else{
        nodal_var_vals[o]=field2[i*nfield[1]+icomp[j]];
          }
        }
        o++;
        printf ("%f\n",nodal_var_vals[m]);
      }while(1);
    }
      }
    }

    // Note: the results for exo tensor format are strangely ordered,
    if (*ncomp==12){
      k=j;
    }else{k=j;}

    int errr = ex_put_nodal_var (exoid, time_step, k+1+countvar, nout, nodal_var_vals);
    if (errr) printf ("ERROR exoselect data for dim %i record %i.\n", j, countvar+j);
  }

  free(nodal_var_vals);
  return;

}