mafillvmain.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 mafillvmain.c:

Go to the source code of this file.

Functions
void	mafillvmain (ITG *nef, ITG *ipnei, ITG *neifa, ITG *neiel, double *vfa, double *xxn, double *area, double *auv, double *adv, ITG *jq, ITG *irow, ITG *nzs, double *bv, double *vel, double *cosa, double *umfa, double *xlet, double *xle, double *gradvfa, double *xxi, double *body, double *volume, ITG *ielfa, char *lakonf, ITG *ifabou, ITG *nbody, double *dtimef, double *velo, double *veloo, double *sel, double *xrlfa, double *gamma, double *xxj, ITG *nactdohinv, double *a1, double *a2, double *a3, double *flux, ITG *icyclic, double *c, ITG *ifatie, ITG *iau6, double *xxni, double *xxnj, ITG *iturbulent, double *gradvel)
void *	mafillvmt (ITG *i)

Variables
static char *	lakonf1
static ITG	num_cpus
static ITG *	nef1
static ITG *	ipnei1
static ITG *	neifa1
static ITG *	neiel1
static ITG *	jq1
static ITG *	irow1
static ITG *	nzs1
static ITG *	ielfa1
static ITG *	ifabou1
static ITG *	nbody1
static ITG *	nactdohinv1
static ITG *	icyclic1
static ITG *	ifatie1
static ITG *	iau61
static ITG *	iturbulent1
static double *	auv1
static double *	adv1
static double *	bv1
static double *	vfa1
static double *	xxn1
static double *	area1
static double *	vel1
static double *	cosa1
static double *	umfa1
static double *	xlet1
static double *	xle1
static double *	gradvfa1
static double *	xxi1
static double *	body1
static double *	volume1
static double *	dtimef1
static double *	velo1
static double *	veloo1
static double *	sel1
static double *	xrlfa1
static double *	gamma1
static double *	xxj1
static double *	a11
static double *	a21
static double *	a31
static double *	flux1
static double *	c1
static double *	xxni1
static double *	xxnj1
static double *	gradvel1

Function Documentation

mafillvmain()

void mafillvmain	(	ITG *	nef,
		ITG *	ipnei,
		ITG *	neifa,
		ITG *	neiel,
		double *	vfa,
		double *	xxn,
		double *	area,
		double *	auv,
		double *	adv,
		ITG *	jq,
		ITG *	irow,
		ITG *	nzs,
		double *	bv,
		double *	vel,
		double *	cosa,
		double *	umfa,
		double *	xlet,
		double *	xle,
		double *	gradvfa,
		double *	xxi,
		double *	body,
		double *	volume,
		ITG *	ielfa,
		char *	lakonf,
		ITG *	ifabou,
		ITG *	nbody,
		double *	dtimef,
		double *	velo,
		double *	veloo,
		double *	sel,
		double *	xrlfa,
		double *	gamma,
		double *	xxj,
		ITG *	nactdohinv,
		double *	a1,
		double *	a2,
		double *	a3,
		double *	flux,
		ITG *	icyclic,
		double *	c,
		ITG *	ifatie,
		ITG *	iau6,
		double *	xxni,
		double *	xxnj,
		ITG *	iturbulent,
		double *	gradvel
	)

                                                                   {

    ITG i,j;
      
    /* variables for multithreading procedure */
    
    ITG sys_cpus,*ithread=NULL;
    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_CFD");
    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;
    }
    }

// next line is to be inserted in a similar way for all other paralell parts

    if(*nef<num_cpus) num_cpus=*nef;
    
    pthread_t tid[num_cpus];

    /* calculating the stiffness and/or mass matrix 
       (symmetric part) */

    nef1=nef;ipnei1=ipnei;neifa1=neifa;neiel1=neiel;vfa1=vfa;xxn1=xxn;
    area1=area;jq1=jq;irow1=irow;nzs1=nzs;vel1=vel;cosa1=cosa;umfa1=umfa;
    xlet1=xlet;xle1=xle;gradvfa1=gradvfa;xxi1=xxi;body1=body;volume1=volume;
    ielfa1=ielfa;lakonf1=lakonf;ifabou1=ifabou;nbody1=nbody;
    dtimef1=dtimef;velo1=velo;veloo1=veloo;sel1=sel;xrlfa1=xrlfa;
    gamma1=gamma;xxj1=xxj;nactdohinv1=nactdohinv;a11=a1;a21=a2;a31=a3;
    flux1=flux;icyclic1=icyclic;c1=c;ifatie1=ifatie;iau61=iau6;
    adv1=adv;auv1=auv;bv1=bv;xxni1=xxni;xxnj1=xxnj;iturbulent1=iturbulent;
    gradvel1=gradvel;
    
    /* 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 *)mafillvmt, (void *)&ithread[i]);
    }
    for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
    
    SFREE(ithread);
  
  return;

}

mafillvmt()

void* mafillvmt

(

ITG *

i

)

                       {

    ITG nefa,nefb,nefdelta;
    
// ceil -> floor

    nefdelta=(ITG)floor(*nef1/(double)num_cpus);
    nefa=*i*nefdelta+1;
    nefb=(*i+1)*nefdelta;
// next line! -> all parallel sections
    if((*i==num_cpus-1)&&(nefb<*nef1)) nefb=*nef1;

    FORTRAN(mafillv,(nef1,ipnei1,neifa1,neiel1,vfa1,xxn1,area1,
        auv1,adv1,jq1,irow1,nzs1,bv1,
            vel1,cosa1,umfa1,xlet1,xle1,gradvfa1,xxi1,
        body1,volume1,ielfa1,lakonf1,ifabou1,nbody1,
        dtimef1,velo1,veloo1,sel1,xrlfa1,gamma1,xxj1,nactdohinv1,
            a11,a21,a31,flux1,&nefa,&nefb,icyclic1,c1,ifatie1,iau61,
             xxni1,xxnj1,iturbulent1,gradvel1));

    return NULL;
}

Variable Documentation

a11

double * a11

static

a21

double * a21

static

a31

double * a31

static

adv1

double * adv1

static

area1

double * area1

static

auv1

double* auv1

static

body1

double * body1

static

bv1

double * bv1

static

c1

double * c1

static

cosa1

double * cosa1

static

dtimef1

double * dtimef1

static

flux1

double * flux1

static

gamma1

double * gamma1

static

gradvel1

double * gradvel1

static

gradvfa1

double * gradvfa1

static

iau61

ITG * iau61

static

icyclic1

ITG * icyclic1

static

ielfa1

ITG * ielfa1

static

ifabou1

ITG * ifabou1

static

ifatie1

ITG * ifatie1

static

ipnei1

ITG * ipnei1

static

irow1

ITG * irow1

static

iturbulent1

ITG * iturbulent1

static

jq1

ITG * jq1

static

lakonf1

char* lakonf1

static

nactdohinv1

ITG * nactdohinv1

static

nbody1

ITG * nbody1

static

nef1

ITG * nef1

static

neiel1

ITG * neiel1

static

neifa1

ITG * neifa1

static

num_cpus

ITG num_cpus

static

nzs1

ITG * nzs1

static

sel1

double * sel1

static

umfa1

double * umfa1

static

vel1

double * vel1

static

velo1

double * velo1

static

veloo1

double * veloo1

static

vfa1

double * vfa1

static

volume1

double * volume1

static

xle1

double * xle1

static

xlet1

double * xlet1

static

xrlfa1

double * xrlfa1

static

xxi1

double * xxi1

static

xxj1

double * xxj1

static

xxn1

double * xxn1

static

xxni1

double * xxni1

static

xxnj1

double * xxnj1

static