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

Go to the source code of this file.

Functions
void	rhspmain (ITG *nef, char *lakonf, ITG *ipnei, ITG *neifa, ITG *neiel, double *vfa, double *area, double *advfa, double *xlet, double *cosa, double *volume, double *au, double *ad, ITG *jq, ITG *irow, double *ap, ITG *ielfa, ITG *ifabou, double *xle, double *b, double *xxn, ITG *neq, ITG *nzs, double *hfa, double *gradpel, double *bp, double *xxi, ITG *neij, double *xlen, ITG *iatleastonepressurebc, double *xxicn)
void *	rhspmt (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 *	ielfa1
static ITG *	ifabou1
static ITG *	neq1
static ITG *	nzs1
static ITG *	neij1
static double *	vfa1
static double *	area1
static double *	advfa1
static double *	xlet1
static double *	cosa1
static double *	volume1
static double *	au1
static double *	ad1
static double *	ap1
static double *	xle1
static double *	b1
static double *	xxn1
static double *	hfa1
static double *	gradpel1
static double *	bp1
static double *	xxi1
static double *	xlen1
static double *	cosb1
static double *	xxicn1

Function Documentation

rhspmain()

void rhspmain	(	ITG *	nef,
		char *	lakonf,
		ITG *	ipnei,
		ITG *	neifa,
		ITG *	neiel,
		double *	vfa,
		double *	area,
		double *	advfa,
		double *	xlet,
		double *	cosa,
		double *	volume,
		double *	au,
		double *	ad,
		ITG *	jq,
		ITG *	irow,
		double *	ap,
		ITG *	ielfa,
		ITG *	ifabou,
		double *	xle,
		double *	b,
		double *	xxn,
		ITG *	neq,
		ITG *	nzs,
		double *	hfa,
		double *	gradpel,
		double *	bp,
		double *	xxi,
		ITG *	neij,
		double *	xlen,
		ITG *	iatleastonepressurebc,
		double *	xxicn
	)

                                                  {

    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;lakonf1=lakonf;ipnei1=ipnei;neifa1=neifa;neiel1=neiel;
    vfa1=vfa;area1=area;advfa1=advfa;xlet1=xlet,cosa1=cosa;volume1=volume;
    jq1=jq;irow1=irow;ap1=ap;ielfa1=ielfa;ifabou1=ifabou;xle1=xle;
    xxn1=xxn;neq1=neq;nzs1=nzs;hfa1=hfa;gradpel1=gradpel;bp1=bp;xxi1=xxi;
    neij1=neij;xlen1=xlen;ad1=ad;au1=au;b1=b;xxicn1=xxicn;
    
    /* 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 *)rhspmt, (void *)&ithread[i]);
    }
    for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
    
    SFREE(ithread);

/*  at least one pressure bc is needed. If none is applied,
    the last dof is set to 0 

     a pressure bc is only recognized if not all velocity degrees of
     freedom are prescribed on the same face */

    if(*iatleastonepressurebc==0) b[*nef-1]=0.;
  
  return;

}

rhspmt()

void* rhspmt

(

ITG *

i

)

                    {

    ITG indexad,indexb,nefa,nefb,nefdelta;
    long long indexau;

    indexad=*i**neq1;
    indexau=(long long)*i**nzs1;
    indexb=*i**neq1;
    
// 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(rhsp,(nef1,lakonf1,ipnei1,neifa1,neiel1,vfa1,area1,
          advfa1,xlet1,cosa1,volume1,au1,ad1,
          jq1,irow1,ap1,ielfa1,ifabou1,xle1,b1,xxn1,
          neq1,nzs1,
          hfa1,gradpel1,bp1,xxi1,neij1,xlen1,&nefa,&nefb,
          xxicn1));

    return NULL;
}

Variable Documentation

ad1

double * ad1

static

advfa1

double * advfa1

static

ap1

double * ap1

static

area1

double * area1

static

au1

double * au1

static

b1

double * b1

static

bp1

double * bp1

static

cosa1

double * cosa1

static

cosb1

double * cosb1

static

gradpel1

double * gradpel1

static

hfa1

double * hfa1

static

ielfa1

ITG * ielfa1

static

ifabou1

ITG * ifabou1

static

ipnei1

ITG * ipnei1

static

irow1

ITG * irow1

static

jq1

ITG * jq1

static

lakonf1

char* lakonf1

static

nef1

ITG * nef1

static

neiel1

ITG * neiel1

static

neifa1

ITG * neifa1

static

neij1

ITG * neij1

static

neq1

ITG * neq1

static

num_cpus

ITG num_cpus

static

nzs1

ITG * nzs1

static

vfa1

double* vfa1

static

volume1

double * volume1

static

xle1

double * xle1

static

xlen1

double * xlen1

static

xlet1

double * xlet1

static

xxi1

double * xxi1

static

xxicn1

double * xxicn1

static

xxn1

double * xxn1

static