calcresidual.c File Reference

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

Include dependency graph for calcresidual.c:

Go to the source code of this file.

Functions
void	calcresidual (ITG *nmethod, ITG *neq, double *b, double *fext, double *f, ITG *iexpl, ITG *nactdof, double *aux2, double *vold, double *vini, double *dtime, double *accold, ITG *nk, double *adb, double *aub, ITG *jq, ITG *irow, ITG *nzl, double *alpha, double *fextini, double *fini, ITG *islavnode, ITG *nslavnode, ITG *mortar, ITG *ntie, double *f_cm, double *f_cs, ITG *mi, ITG *nzs, ITG *nasym, ITG *idamping, double *veold, double *adc, double *auc, double *cvini, double *cv)

Function Documentation

calcresidual()

void calcresidual	(	ITG *	nmethod,
		ITG *	neq,
		double *	b,
		double *	fext,
		double *	f,
		ITG *	iexpl,
		ITG *	nactdof,
		double *	aux2,
		double *	vold,
		double *	vini,
		double *	dtime,
		double *	accold,
		ITG *	nk,
		double *	adb,
		double *	aub,
		ITG *	jq,
		ITG *	irow,
		ITG *	nzl,
		double *	alpha,
		double *	fextini,
		double *	fini,
		ITG *	islavnode,
		ITG *	nslavnode,
		ITG *	mortar,
		ITG *	ntie,
		double *	f_cm,
		double *	f_cs,
		ITG *	mi,
		ITG *	nzs,
		ITG *	nasym,
		ITG *	idamping,
		double *	veold,
		double *	adc,
		double *	auc,
		double *	cvini,
		double *	cv
	)

                                                                       {

    ITG j,k,mt=mi[1]+1;
    double scal1;
      
    /* residual for a static analysis */
      
    if(*nmethod!=4){
    for(k=0;k<neq[1];++k){
        b[k]=fext[k]-f[k];
    }
    }
      
    /* residual for implicit dynamics */
      
    else if(*iexpl<=1){
    for(k=0;k<*nk;++k){
        if(nactdof[mt*k]>0){
        aux2[nactdof[mt*k]-1]=(vold[mt*k]-vini[mt*k])/(*dtime);}
        for(j=1;j<mt;++j){
        if(nactdof[mt*k+j]>0){aux2[nactdof[mt*k+j]-1]=accold[mt*k+j];}
        }
    }
    if(*nasym==0){
        FORTRAN(op,(&neq[1],aux2,b,adb,aub,jq,irow)); 
    }else{
        FORTRAN(opas,(&neq[1],aux2,b,adb,aub,jq,irow,nzs)); 
    }
    scal1=1.+*alpha;
    for(k=0;k<neq[0];++k){
        b[k]=scal1*(fext[k]-f[k])-*alpha*(fextini[k]-fini[k])-b[k];
    } 
    for(k=neq[0];k<neq[1];++k){
        b[k]=fext[k]-f[k]-b[k];
    } 

    /* correction for damping */

    if(*idamping==1){
        for(k=0;k<*nk;++k){
        if(nactdof[mt*k]>0){aux2[nactdof[mt*k]-1]=0.;}
        for(j=1;j<mt;++j){
            if(nactdof[mt*k+j]>0){
            aux2[nactdof[mt*k+j]-1]=veold[mt*k+j];}
        }
        }
        if(*nasym==0){
        FORTRAN(op,(&neq[1],aux2,cv,adc,auc,jq,irow));
        }else{
        FORTRAN(opas,(&neq[1],aux2,cv,adc,auc,jq,irow,nzs)); 
        }
        for(k=0;k<neq[0];++k){
        b[k]-=scal1*cv[k]-*alpha*cvini[k];
        }
    }
    }

    /* residual for explicit dynamics */
    
    else{
    for(k=0;k<*nk;++k){
        if(nactdof[mt*k]>0){
        aux2[nactdof[mt*k]-1]=(vold[mt*k]-vini[mt*k])/(*dtime);}
        for(j=1;j<mt;++j){
        if(nactdof[mt*k+j]>0){aux2[nactdof[mt*k+j]-1]=accold[mt*k+j];}
        }
    }
    scal1=1.+*alpha;
    for(k=0;k<neq[0];++k){
        b[k]=scal1*(fext[k]-f[k])-*alpha*(fextini[k]-fini[k])
        -adb[k]*aux2[k];
    } 
    for(k=neq[0];k<neq[1];++k){
        b[k]=fext[k]-f[k]-adb[k]*aux2[k];
    } 
    }

    return;
}