thicknessmain.c File Reference

#include <unistd.h>
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <pthread.h>
#include "CalculiX.h"

Include dependency graph for thicknessmain.c:

Go to the source code of this file.

Functions
void	thicknessmain (double *co, double *dgdx, ITG *nobject, ITG *nk, ITG *nodedesi, ITG *ndesi, char *objectset, ITG *ipkon, ITG *kon, char *lakon, char *set, ITG *nset, ITG *istartset, ITG *iendset, ITG *ialset, ITG *iobject, ITG *nodedesiinv, double *dgdxglob)
void *	thicknessmt (ITG *i)

Variables
static char *	objectset1
static ITG *	nobject1
static ITG *	nodedesiboun1
static ITG	ndesiboun1
static ITG *	nx1
static ITG *	ny1
static ITG *	nz1
static ITG	num_cpus
static ITG	ifree1
static ITG *	iobject1
static ITG *	ndesi1
static ITG *	nk1
static double *	dgdx1
static double *	xo1
static double *	yo1
static double *	zo1
static double *	x1
static double *	yy1
static double *	z1
static double *	co1
static double *	dgdxglob1

Function Documentation

thicknessmain()

void thicknessmain	(	double *	co,
		double *	dgdx,
		ITG *	nobject,
		ITG *	nk,
		ITG *	nodedesi,
		ITG *	ndesi,
		char *	objectset,
		ITG *	ipkon,
		ITG *	kon,
		char *	lakon,
		char *	set,
		ITG *	nset,
		ITG *	istartset,
		ITG *	iendset,
		ITG *	ialset,
		ITG *	iobject,
		ITG *	nodedesiinv,
		double *	dgdxglob
	)

                         {
        
    /* calculation of distance between nodes */

    ITG *nx=NULL,*ny=NULL,*nz=NULL,ifree,i,*ithread=NULL,ndesiboun,
        *nodedesiboun=NULL;
    
    double *xo=NULL,*yo=NULL,*zo=NULL,*x=NULL,*y=NULL,*z=NULL;

    /* prepare for near3d */
    
    NNEW(xo,double,*nk);
    NNEW(yo,double,*nk);
    NNEW(zo,double,*nk);
    NNEW(x,double,*nk);
    NNEW(y,double,*nk);
    NNEW(z,double,*nk);
    NNEW(nx,ITG,*nk);
    NNEW(ny,ITG,*nk);
    NNEW(nz,ITG,*nk);
    NNEW(nodedesiboun,ITG,*ndesi);   
    
    FORTRAN(prethickness,(co,xo,yo,zo,x,y,z,nx,ny,nz,&ifree,nodedesiinv,
                &ndesiboun,nodedesiboun,set,nset,objectset,
                iobject,istartset,iendset,ialset));
                 
    RENEW(nodedesiboun,ITG,ndesiboun);
    
    /* variables for multithreading procedure */
    
    ITG sys_cpus;
    char *env,*envloc,*envsys;
    
    num_cpus = 0;
    sys_cpus=0;
    
    /* explicit user declaration prevails */
    
    envsys=getenv("NUMBER_OF_CPUS");
    if(envsys){
       sys_cpus=atoi(envsys);
       if(sys_cpus<0) sys_cpus=0;
    }
    
    /* automatic detection of available number of processors */
    
    if(sys_cpus==0){
        sys_cpus = getSystemCPUs();
        if(sys_cpus<1) sys_cpus=1;
    }
    
    /* local declaration prevails, if strictly positive */
    
    envloc = getenv("CCX_NPROC_FILTER");
    if(envloc){
       num_cpus=atoi(envloc);
       if(num_cpus<0){
          num_cpus=0;
       }else if(num_cpus>sys_cpus){
          num_cpus=sys_cpus;
       }

    }
    
    /* else global declaration, if any, applies */
    
    env = getenv("OMP_NUM_THREADS");
    if(num_cpus==0){
       if (env)
          num_cpus = atoi(env);
       if (num_cpus < 1) {
          num_cpus=1;
       }else if(num_cpus>sys_cpus){
         num_cpus=sys_cpus;
       }
    }
    
    /* check that the number of cpus does not supercede the number
       of design variables in nodedesiboun */
    
    if(ndesiboun<num_cpus) num_cpus=ndesiboun;
    
    pthread_t tid[num_cpus];

    dgdx1=dgdx;nobject1=nobject;nodedesiboun1=nodedesiboun;
    ndesiboun1=ndesiboun;objectset1=objectset;xo1=xo;yo1=yo;zo1=zo;
    x1=x;yy1=y;z1=z;nx1=nx;ny1=ny;nz1=nz;ifree1=ifree;co1=co;  
    iobject1=iobject;ndesi1=ndesi;dgdxglob1=dgdxglob;nk1=nk;

    /* assessment of actual wallthickness */
    /* create threads and wait */
  
    NNEW(ithread,ITG,num_cpus);
    for(i=0; i<num_cpus; i++)  {
       ithread[i]=i;
       pthread_create(&tid[i], NULL, (void *)thicknessmt, (void *)&ithread[i]);
    }
    for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);

    SFREE(xo);SFREE(yo);SFREE(zo);
    SFREE(x);SFREE(y);SFREE(z);SFREE(nx);SFREE(ny);SFREE(nz);
    SFREE(ithread);SFREE(nodedesiboun);   
                                     
    return;
    
} 

thicknessmt()

void* thicknessmt

(

ITG *

i

)

                         {

    ITG indexr,ndesia,ndesib,ndesidelta;

    indexr=*i*ifree1;
    
    ndesidelta=(ITG)ceil(ndesiboun1/(double)num_cpus);
    ndesia=*i*ndesidelta+1;
    ndesib=(*i+1)*ndesidelta;
    if(ndesib>ndesiboun1) ndesib=ndesiboun1;
    
    //printf("indexr=%" ITGFORMAT","ndesia=%" ITGFORMAT",ndesib=%" ITGFORMAT"\n",indexr,ndesia,ndesib);

    FORTRAN(thickness,(dgdx1,nobject1,nodedesiboun1,&ndesiboun1,objectset1,
                        xo1,yo1,zo1,x1,yy1,z1,nx1,ny1,nz1,co1,&ifree1,
                        &ndesia,&ndesib,iobject1,ndesi1,dgdxglob1,nk1));

    return NULL;
}

Variable Documentation

co1

double * co1

static

dgdx1

double* dgdx1

static

dgdxglob1

double * dgdxglob1

static

ifree1

ITG ifree1

static

iobject1

ITG * iobject1

static

ndesi1

ITG * ndesi1

static

ndesiboun1

ITG ndesiboun1

static

nk1

ITG * nk1

static

nobject1

ITG* nobject1

static

nodedesiboun1

ITG * nodedesiboun1

static

num_cpus

ITG num_cpus

static

nx1

ITG * nx1

static

ny1

ITG * ny1

static

nz1

ITG * nz1

static

objectset1

char* objectset1

static

x1

double * x1

static

xo1

double * xo1

static

yo1

double * yo1

static

yy1

double * yy1

static

z1

double * z1

static

zo1

double * zo1

static