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

Go to the source code of this file.

Functions
void	mafillsmmain_se (double *co, ITG *nk, ITG *kon, ITG *ipkon, char *lakon, ITG *ne, ITG *nodeboun, ITG *ndirboun, double *xboun, ITG *nboun, ITG *ipompc, ITG *nodempc, double *coefmpc, ITG *nmpc, ITG *nodeforc, ITG *ndirforc, double *xforc, ITG *nforc, ITG *nelemload, char *sideload, double *xload, ITG *nload, double *xbody, ITG *ipobody, ITG *nbody, double *cgr, ITG *nactdof, ITG *neq, ITG *nmethod, ITG *ikmpc, ITG *ilmpc, ITG *ikboun, ITG *ilboun, double *elcon, ITG *nelcon, double *rhcon, ITG *nrhcon, double *alcon, ITG *nalcon, double *alzero, ITG *ielmat, ITG *ielorien, ITG *norien, double *orab, ITG *ntmat_, double *t0, double *t1, ITG *ithermal, double *prestr, ITG *iprestr, double *vold, ITG *iperturb, double *sti, double *stx, ITG *iexpl, double *plicon, ITG *nplicon, double *plkcon, ITG *nplkcon, double *xstiff, ITG *npmat_, double *dtime, char *matname, ITG *mi, ITG *ncmat_, ITG *mass, ITG *stiffness, ITG *buckling, ITG *rhsi, ITG *intscheme, double *physcon, double *shcon, ITG *nshcon, double *cocon, ITG *ncocon, double *ttime, double *time, ITG *istep, ITG *iinc, ITG *coriolis, ITG *ibody, double *xloadold, double *reltime, double *veold, double *springarea, ITG *nstate_, double *xstateini, double *xstate, double *thicke, ITG *integerglob, double *doubleglob, char *tieset, ITG *istartset, ITG *iendset, ITG *ialset, ITG *ntie, ITG *nasym, double *pslavsurf, double *pmastsurf, ITG *mortar, double *clearini, ITG *ielprop, double *prop, ITG *ne0, double *fnext, double *distmin, ITG *ndesi, ITG *nodedesi, double *df, ITG *nzss, ITG *jqs, ITG *irows, ITG *icoordinate, double *dxstiff, double *xdesi, ITG *istartelem, ITG *ialelem, double *v, double *sigma, ITG *cyclicsymmetry, char *labmpc, ITG *ics, double *cs, ITG *mcs, ITG *ieigenfrequency)
void *	mafillsmsemt (ITG *i)

Variables
static char *	lakon1
static char *	sideload1
static char *	matname1
static char *	tieset1
static char *	labmpc1
static ITG *	nk1
static ITG *	kon1
static ITG *	ipkon1
static ITG *	ne1
static ITG *	nodeboun1
static ITG *	ndirboun1
static ITG *	nboun1
static ITG *	ipompc1
static ITG *	nodempc1
static ITG *	nmpc1
static ITG *	nodeforc1
static ITG *	ndirforc1
static ITG *	nforc1
static ITG *	nelemload1
static ITG *	nload1
static ITG *	ipobody1
static ITG *	nbody1
static ITG *	nactdof1
static ITG *	neq1
static ITG *	nmethod1 =NULL
static ITG *	ikmpc1
static ITG *	ilmpc1
static ITG *	ikboun1
static ITG *	ilboun1
static ITG *	nelcon1
static ITG *	nrhcon1
static ITG *	nalcon1
static ITG *	ielmat1
static ITG *	ielorien1
static ITG *	norien1
static ITG *	ntmat1_
static ITG *	ithermal1
static ITG *	iprestr1
static ITG *	iperturb1
static ITG *	iexpl1
static ITG *	nplicon1
static ITG *	nplkcon1
static ITG *	npmat1_
static ITG *	mi1
static ITG *	ncmat1_
static ITG *	mass1
static ITG *	stiffness1
static ITG *	buckling1
static ITG *	rhsi1
static ITG *	intscheme1
static ITG *	nshcon1
static ITG *	ncocon1
static ITG *	istep1
static ITG *	iinc1
static ITG *	coriolis1
static ITG *	ibody1
static ITG *	nstate1_
static ITG *	integerglob1
static ITG *	istartset1
static ITG *	iendset1
static ITG *	ialset1
static ITG *	ntie1
static ITG *	nasym1
static ITG *	mortar1
static ITG *	ielprop1
static ITG *	ne01
static ITG	num_cpus
static ITG *	ndesi1
static ITG *	nodedesi1
static ITG *	nzss1
static ITG *	jqs1
static ITG *	irows1
static ITG *	icoordinate1
static ITG *	istartelem1
static ITG *	ialelem1
static ITG *	cyclicsymmetry1
static ITG *	ics1
static ITG *	mcs1
static ITG *	ieigenfrequency1
static double *	co1
static double *	xboun1
static double *	coefmpc1
static double *	xforc1
static double *	xload1
static double *	xbody1
static double *	cgr1
static double *	elcon1
static double *	rhcon1
static double *	alcon1
static double *	alzero1
static double *	xdesi1
static double *	v1
static double *	sigma1
static double *	orab1
static double *	t01
static double *	t11
static double *	prestr1
static double *	vold1
static double *	sti1
static double *	stx1
static double *	cs1
static double *	plicon1
static double *	plkcon1
static double *	xstiff1
static double *	dtime1
static double *	physcon1
static double *	shcon1
static double *	cocon1
static double *	ttime1
static double *	time1
static double *	xloadold1
static double *	reltime1
static double *	veold1
static double *	springarea1
static double *	xstateini1
static double *	xstate1
static double *	thicke1
static double *	doubleglob1
static double *	pslavsurf1
static double *	pmastsurf1
static double *	clearini1
static double *	prop1
static double *	distmin1
static double *	df1
static double *	dfl1
static double *	dxstiff1

Function Documentation

mafillsmmain_se()

void mafillsmmain_se	(	double *	co,
		ITG *	nk,
		ITG *	kon,
		ITG *	ipkon,
		char *	lakon,
		ITG *	ne,
		ITG *	nodeboun,
		ITG *	ndirboun,
		double *	xboun,
		ITG *	nboun,
		ITG *	ipompc,
		ITG *	nodempc,
		double *	coefmpc,
		ITG *	nmpc,
		ITG *	nodeforc,
		ITG *	ndirforc,
		double *	xforc,
		ITG *	nforc,
		ITG *	nelemload,
		char *	sideload,
		double *	xload,
		ITG *	nload,
		double *	xbody,
		ITG *	ipobody,
		ITG *	nbody,
		double *	cgr,
		ITG *	nactdof,
		ITG *	neq,
		ITG *	nmethod,
		ITG *	ikmpc,
		ITG *	ilmpc,
		ITG *	ikboun,
		ITG *	ilboun,
		double *	elcon,
		ITG *	nelcon,
		double *	rhcon,
		ITG *	nrhcon,
		double *	alcon,
		ITG *	nalcon,
		double *	alzero,
		ITG *	ielmat,
		ITG *	ielorien,
		ITG *	norien,
		double *	orab,
		ITG *	ntmat_,
		double *	t0,
		double *	t1,
		ITG *	ithermal,
		double *	prestr,
		ITG *	iprestr,
		double *	vold,
		ITG *	iperturb,
		double *	sti,
		double *	stx,
		ITG *	iexpl,
		double *	plicon,
		ITG *	nplicon,
		double *	plkcon,
		ITG *	nplkcon,
		double *	xstiff,
		ITG *	npmat_,
		double *	dtime,
		char *	matname,
		ITG *	mi,
		ITG *	ncmat_,
		ITG *	mass,
		ITG *	stiffness,
		ITG *	buckling,
		ITG *	rhsi,
		ITG *	intscheme,
		double *	physcon,
		double *	shcon,
		ITG *	nshcon,
		double *	cocon,
		ITG *	ncocon,
		double *	ttime,
		double *	time,
		ITG *	istep,
		ITG *	iinc,
		ITG *	coriolis,
		ITG *	ibody,
		double *	xloadold,
		double *	reltime,
		double *	veold,
		double *	springarea,
		ITG *	nstate_,
		double *	xstateini,
		double *	xstate,
		double *	thicke,
		ITG *	integerglob,
		double *	doubleglob,
		char *	tieset,
		ITG *	istartset,
		ITG *	iendset,
		ITG *	ialset,
		ITG *	ntie,
		ITG *	nasym,
		double *	pslavsurf,
		double *	pmastsurf,
		ITG *	mortar,
		double *	clearini,
		ITG *	ielprop,
		double *	prop,
		ITG *	ne0,
		double *	fnext,
		double *	distmin,
		ITG *	ndesi,
		ITG *	nodedesi,
		double *	df,
		ITG *	nzss,
		ITG *	jqs,
		ITG *	irows,
		ITG *	icoordinate,
		double *	dxstiff,
		double *	xdesi,
		ITG *	istartelem,
		ITG *	ialelem,
		double *	v,
		double *	sigma,
		ITG *	cyclicsymmetry,
		char *	labmpc,
		ITG *	ics,
		double *	cs,
		ITG *	mcs,
		ITG *	ieigenfrequency
	)

                                         {
           
    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;
    }

//    sys_cpus=1;

    /* 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_STIFFNESS");
    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(*ne<num_cpus) num_cpus=*ne;
    
    pthread_t tid[num_cpus];

    /* allocating fields for mass and stiffness matrix and sensitivity matrix */
    
    if(!*cyclicsymmetry){
    NNEW(dfl1,double,num_cpus*60**ndesi);
    NNEW(df1,double,num_cpus**nzss);
    }else{
    NNEW(dfl1,double,num_cpus*120**ndesi);
    NNEW(df1,double,num_cpus*2**nzss);
    }

    /* allocating memory for nmethod; if the Jacobian determinant
       in any of the elements is nonpositive, nmethod is set to
       zero */

    NNEW(nmethod1,ITG,num_cpus);
    for(j=0;j<num_cpus;j++){
    nmethod1[j]=*nmethod;
    }

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

    co1=co;nk1=nk;kon1=kon;ipkon1=ipkon;lakon1=lakon;ne1=ne;
    nodeboun1=nodeboun;ndirboun1=ndirboun;xboun1=xboun;
    nboun1=nboun;ipompc1=ipompc;nodempc1=nodempc;coefmpc1=coefmpc;
    nmpc1=nmpc;nodeforc1=nodeforc;ndirforc1=ndirforc;xforc1=xforc;
    nforc1=nforc;nelemload1=nelemload;sideload1=sideload;xload1=xload;
    nload1=nload;xbody1=xbody;ipobody1=ipobody;nbody1=nbody;
    cgr1=cgr;nactdof1=nactdof;neq1=neq;
    ikmpc1=ikmpc;ilmpc1=ilmpc;ikboun1=ikboun;
    ilboun1=ilboun;elcon1=elcon;nelcon1=nelcon;rhcon1=rhcon;
    nrhcon1=nrhcon;alcon1=alcon;nalcon1=nalcon;alzero1=alzero;
    ielmat1=ielmat;ielorien1=ielorien;norien1=norien;orab1=orab;
    ntmat1_=ntmat_;t01=t0;t11=t1;ithermal1=ithermal;prestr1=prestr;
    iprestr1=iprestr;vold1=vold;iperturb1=iperturb;sti1=sti;
    stx1=stx;iexpl1=iexpl;plicon1=plicon;nplicon1=nplicon;
    plkcon1=plkcon;nplkcon1=nplkcon;xstiff1=xstiff;npmat1_=npmat_;
    dtime1=dtime;matname1=matname;mi1=mi;ncmat1_=ncmat_;mass1=mass;
    stiffness1=stiffness;buckling1=buckling;rhsi1=rhsi;intscheme1=intscheme;
    physcon1=physcon;shcon1=shcon;nshcon1=nshcon;cocon1=cocon;
    ncocon1=ncocon;ttime1=ttime;time1=time;istep1=istep;iinc1=iinc;
    coriolis1=coriolis;ibody1=ibody;xloadold1=xloadold;reltime1=reltime;
    veold1=veold;springarea1=springarea;nstate1_=nstate_;xstateini1=xstateini;
    xstate1=xstate;thicke1=thicke;integerglob1=integerglob;
    doubleglob1=doubleglob;tieset1=tieset;istartset1=istartset;
    iendset1=iendset;ialset1=ialset;ntie1=ntie;nasym1=nasym;
    pslavsurf1=pslavsurf;pmastsurf1=pmastsurf;mortar1=mortar;
    clearini1=clearini;ielprop1=ielprop;prop1=prop;ne01=ne0;
    distmin1=distmin;ndesi1=ndesi;nodedesi1=nodedesi;v1=v;
    nzss1=nzss;jqs1=jqs;irows1=irows;icoordinate1=icoordinate;
    dxstiff1=dxstiff;xdesi1=xdesi;istartelem1=istartelem;ialelem1=ialelem;
    sigma1=sigma;cyclicsymmetry1=cyclicsymmetry;labmpc1=labmpc;
    ics1=ics;cs1=cs;mcs1=mcs;ieigenfrequency1=ieigenfrequency;

    /* calculating the stiffness/mass sensitivity */
    
    printf(" Using up to %" ITGFORMAT " cpu(s) for the calculation of the sensitivity of the external forces \n and/or the element stiffness matrices.\n\n", num_cpus);
    
    /* 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 *)mafillsmsemt, (void *)&ithread[i]);
    }
    for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
    
    SFREE(ithread);

    for(j=0;j<num_cpus;j++){
    if(nmethod1[j]==0){
        *nmethod=0;
        break;
    }
    }
    SFREE(nmethod1);
    
    /* passing of df */

    if(!*cyclicsymmetry){
//  if((iperturb[1]==1)&&(*ieigenfrequency!=1)){
    if(*ieigenfrequency!=1){

            /* nonlinear geometric: add df to df from results_se */

        for(i=0;i<*nzss;i++){
        df[i]+=df1[i];
        }
    }else{

        for(i=0;i<*nzss;i++){
        df[i]=df1[i];
        }
    }
    for(i=0;i<*nzss;i++){
        for(j=1;j<num_cpus;j++){
        df[i]+=df1[i+j**nzss];
        }
    }
    }else{
    for(i=0;i<2**nzss;i++){
        df[i]=df1[i];
    }
    for(i=0;i<2**nzss;i++){
        for(j=1;j<num_cpus;j++){
        df[i]+=df1[i+j*2**nzss];
        }
    }
    }
    SFREE(df1);SFREE(dfl1);
        
        return;

}

mafillsmsemt()

void* mafillsmsemt

(

ITG *

i

)

                          {

    ITG indexdf,indexdfl,nea,neb,nedelta;

    if(!*cyclicsymmetry1){
    indexdf=*i**nzss1;
    indexdfl=*i*60**ndesi1;
    }else{
    indexdf=*i*2**nzss1;
    indexdfl=*i*120**ndesi1;
    }

    nedelta=(ITG)floor(*ne1/(double)num_cpus);
    nea=*i*nedelta+1;
    neb=(*i+1)*nedelta;
    if((*i==num_cpus-1)&&(neb<*ne1)) neb=*ne1;

    if(!*cyclicsymmetry1){
    FORTRAN(mafillsmse,(co1,kon1,ipkon1,lakon1,ne1,ipompc1,nodempc1,
        coefmpc1,nmpc1,nelemload1,sideload1,xload1,nload1,xbody1,
        ipobody1,nbody1,cgr1,nactdof1,neq1,&nmethod1[*i],ikmpc1,
        ilmpc1,elcon1,nelcon1,rhcon1,nrhcon1,alcon1,nalcon1,alzero1,
        ielmat1,ielorien1,norien1,orab1,ntmat1_,t01,t11,ithermal1,
        iprestr1,vold1,iperturb1,sti1,stx1,iexpl1,plicon1,nplicon1,
            plkcon1,nplkcon1,xstiff1,npmat1_,dtime1,matname1,mi1,
            ncmat1_,mass1,stiffness1,buckling1,rhsi1,intscheme1,physcon1,
            ttime1,time1,istep1,iinc1,coriolis1,ibody1,xloadold1,
        reltime1,veold1,springarea1,nstate1_,xstateini1,xstate1,
            thicke1,integerglob1,doubleglob1,tieset1,istartset1,iendset1,
        ialset1,ntie1,nasym1,pslavsurf1,pmastsurf1,mortar1,clearini1,
        ielprop1,prop1,ne01,&nea,&neb,distmin1,ndesi1,nodedesi1,
        &df1[indexdf],jqs1,irows1,&dfl1[indexdfl],
        icoordinate1,dxstiff1,xdesi1,istartelem1,ialelem1,v1,sigma1,
            ieigenfrequency1));
    }else{
    FORTRAN(mafillsmcsse,(co1,kon1,ipkon1,lakon1,ne1,ipompc1,nodempc1,
        coefmpc1,nmpc1,nelemload1,sideload1,xload1,nload1,xbody1,
        ipobody1,nbody1,cgr1,nactdof1,neq1,&nmethod1[*i],ikmpc1,
        ilmpc1,elcon1,nelcon1,rhcon1,nrhcon1,alcon1,nalcon1,alzero1,
        ielmat1,ielorien1,norien1,orab1,ntmat1_,t01,t11,ithermal1,
        iprestr1,vold1,iperturb1,sti1,stx1,iexpl1,plicon1,nplicon1,
            plkcon1,nplkcon1,xstiff1,npmat1_,dtime1,matname1,mi1,
            ncmat1_,mass1,stiffness1,buckling1,rhsi1,intscheme1,physcon1,
            ttime1,time1,istep1,iinc1,coriolis1,ibody1,xloadold1,
        reltime1,veold1,springarea1,nstate1_,xstateini1,xstate1,
            thicke1,integerglob1,doubleglob1,tieset1,istartset1,iendset1,
        ialset1,ntie1,nasym1,pslavsurf1,pmastsurf1,mortar1,clearini1,
        ielprop1,prop1,ne01,&nea,&neb,distmin1,ndesi1,nodedesi1,
        &df1[indexdf],jqs1,irows1,&dfl1[indexdfl],
        icoordinate1,dxstiff1,xdesi1,istartelem1,ialelem1,v1,sigma1,
        labmpc1,ics1,cs1,mcs1,nk1,nzss1));
    }
    
        
    return NULL;
}

Variable Documentation

alcon1

double * alcon1

static

alzero1

double * alzero1

static

buckling1

ITG * buckling1

static

cgr1

double * cgr1

static

clearini1

double * clearini1

static

co1

double* co1

static

cocon1

double * cocon1

static

coefmpc1

double * coefmpc1

static

coriolis1

ITG * coriolis1

static

cs1

double * cs1

static

cyclicsymmetry1

ITG * cyclicsymmetry1

static

df1

double * df1

static

dfl1

double * dfl1

static

distmin1

double * distmin1

static

doubleglob1

double * doubleglob1

static

dtime1

double * dtime1

static

dxstiff1

double * dxstiff1

static

elcon1

double * elcon1

static

ialelem1

ITG * ialelem1

static

ialset1

ITG * ialset1

static

ibody1

ITG * ibody1

static

icoordinate1

ITG * icoordinate1

static

ics1

ITG * ics1

static

ieigenfrequency1

ITG * ieigenfrequency1

static

ielmat1

ITG * ielmat1

static

ielorien1

ITG * ielorien1

static

ielprop1

ITG * ielprop1

static

iendset1

ITG * iendset1

static

iexpl1

ITG * iexpl1

static

iinc1

ITG * iinc1

static

ikboun1

ITG * ikboun1

static

ikmpc1

ITG * ikmpc1

static

ilboun1

ITG * ilboun1

static

ilmpc1

ITG * ilmpc1

static

integerglob1

ITG * integerglob1

static

intscheme1

ITG * intscheme1

static

iperturb1

ITG * iperturb1

static

ipkon1

ITG * ipkon1

static

ipobody1

ITG * ipobody1

static

ipompc1

ITG * ipompc1

static

iprestr1

ITG * iprestr1

static

irows1

ITG * irows1

static

istartelem1

ITG * istartelem1

static

istartset1

ITG * istartset1

static

istep1

ITG * istep1

static

ithermal1

ITG * ithermal1

static

jqs1

ITG * jqs1

static

kon1

ITG * kon1

static

labmpc1

char * labmpc1

static

lakon1

char* lakon1

static

mass1

ITG * mass1

static

matname1

char * matname1

static

mcs1

ITG * mcs1

static

mi1

ITG * mi1

static

mortar1

ITG * mortar1

static

nactdof1

ITG * nactdof1

static

nalcon1

ITG * nalcon1

static

nasym1

ITG * nasym1

static

nbody1

ITG * nbody1

static

nboun1

ITG * nboun1

static

ncmat1_

ITG * ncmat1_

static

ncocon1

ITG * ncocon1

static

ndesi1

ITG * ndesi1

static

ndirboun1

ITG * ndirboun1

static

ndirforc1

ITG * ndirforc1

static

ne01

ITG * ne01

static

ne1

ITG * ne1

static

nelcon1

ITG * nelcon1

static

nelemload1

ITG * nelemload1

static

neq1

ITG * neq1

static

nforc1

ITG * nforc1

static

nk1

ITG* nk1

static

nload1

ITG * nload1

static

nmethod1

ITG * nmethod1 =NULL

static

nmpc1

ITG * nmpc1

static

nodeboun1

ITG * nodeboun1

static

nodedesi1

ITG * nodedesi1

static

nodeforc1

ITG * nodeforc1

static

nodempc1

ITG * nodempc1

static

norien1

ITG * norien1

static

nplicon1

ITG * nplicon1

static

nplkcon1

ITG * nplkcon1

static

npmat1_

ITG * npmat1_

static

nrhcon1

ITG * nrhcon1

static

nshcon1

ITG * nshcon1

static

nstate1_

ITG * nstate1_

static

ntie1

ITG * ntie1

static

ntmat1_

ITG * ntmat1_

static

num_cpus

ITG num_cpus

static

nzss1

ITG * nzss1

static

orab1

double * orab1

static

physcon1

double * physcon1

static

plicon1

double * plicon1

static

plkcon1

double * plkcon1

static

pmastsurf1

double * pmastsurf1

static

prestr1

double * prestr1

static

prop1

double * prop1

static

pslavsurf1

double * pslavsurf1

static

reltime1

double * reltime1

static

rhcon1

double * rhcon1

static

rhsi1

ITG * rhsi1

static

shcon1

double * shcon1

static

sideload1

char * sideload1

static

sigma1

double * sigma1

static

springarea1

double * springarea1

static

sti1

double * sti1

static

stiffness1

ITG * stiffness1

static

stx1

double * stx1

static

t01

double * t01

static

t11

double * t11

static

thicke1

double * thicke1

static

tieset1

char * tieset1

static

time1

double * time1

static

ttime1

double * ttime1

static

v1

double * v1

static

veold1

double * veold1

static

vold1

double * vold1

static

xbody1

double * xbody1

static

xboun1

double * xboun1

static

xdesi1

double * xdesi1

static

xforc1

double * xforc1

static

xload1

double * xload1

static

xloadold1

double * xloadold1

static

xstate1

double * xstate1

static

xstateini1

double * xstateini1

static

xstiff1

double * xstiff1

static