SmallStrainPlasticity.cpp File Reference

#include <MoFEM.hpp>
#include <boost/numeric/ublas/vector_proxy.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/matrix_proxy.hpp>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/symmetric.hpp>
#include <adolc/adolc.h>
#include <SmallStrainPlasticity.hpp>

Go to the source code of this file.

Functions
PetscErrorCode	SmallStrainPlasticityfRes (SNES snes, Vec chi, Vec r, void *ctx)
PetscErrorCode	SmallStrainPlasticityfJac (SNES snes, Vec chi, Mat A, Mat, void *ctx)

Function Documentation

SmallStrainPlasticityfJac()

PetscErrorCode SmallStrainPlasticityfJac	(	SNES	snes,
		Vec	chi,
		Mat	A,
		Mat	,
		void *	ctx
	)

Definition at line 1207 of file SmallStrainPlasticity.cpp.

                                                                                {
  PetscFunctionBegin;
  SmallStrainPlasticity::ClosestPointProjection *cp;
  cp = (SmallStrainPlasticity::ClosestPointProjection*)ctx;
  
  {
    #if PETSC_VERSION_GE(3,6,0)
    const double *array;
    ierr = VecGetArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    double *array;
    ierr = VecGetArray(chi,&array); CHKERRQ(ierr);
    #endif
    for(int ii = 0;ii<6;ii++) {
      cp->plasticStrain[ii] = array[ii];
    }
    for(unsigned int ii = 0;ii<cp->internalVariables.size();ii++) {
      cp->internalVariables[ii] = array[6+ii];
    }
    cp->deltaGamma = array[6+cp->internalVariables.size()];
    #if PETSC_VERSION_GE(3,6,0)
    ierr = VecRestoreArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    ierr = VecRestoreArray(chi,&array); CHKERRQ(ierr);
    #endif
  }
  ierr = cp->evaluatePotentials(); CHKERRQ(ierr);
  ierr = cp->cAlculateA(); CHKERRQ(ierr);
  PetscFunctionReturn(0);
}

SmallStrainPlasticityfRes()

PetscErrorCode SmallStrainPlasticityfRes	(	SNES	snes,
		Vec	chi,
		Vec	r,
		void *	ctx
	)

Definition at line 1154 of file SmallStrainPlasticity.cpp.

                                                                            {
  PetscFunctionBegin;
  SmallStrainPlasticity::ClosestPointProjection *cp;
  cp = (SmallStrainPlasticity::ClosestPointProjection*)ctx;
  
  //ierr = VecView(chi,PETSC_VIEWER_STDOUT_WORLD); CHKERRQ(ierr);
  {
    #if PETSC_VERSION_GE(3,6,0)
    const double *array;
    ierr = VecGetArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    const double *array;
    ierr = VecGetArray(chi,&array); CHKERRQ(ierr);
    #endif
    for(int ii = 0;ii<6;ii++) {
      cp->plasticStrain[ii] = array[ii];
    }
    for(unsigned int ii = 0;ii<cp->internalVariables.size();ii++) {
      cp->internalVariables[ii] = array[6+ii];
    }
    cp->deltaGamma = array[6+cp->internalVariables.size()];
    #if PETSC_VERSION_GE(3,6,0)
    ierr = VecRestoreArrayRead(chi,&array); CHKERRQ(ierr);
    #else
    ierr = VecRestoreArray(chi,&array); CHKERRQ(ierr);
    #endif
  }
  ierr = cp->evaluatePotentials(); CHKERRQ(ierr);
  ierr = cp->cAlculateR(r); CHKERRQ(ierr);
  // ierr = VecView(r,PETSC_VIEWER_STDOUT_WORLD); CHKERRQ(ierr);
  // cout << "sTress " << cp->sTress << endl;
  // cout << "internalVariables " << cp->internalVariables << endl;
  // cout << "C " << cp->C << endl;
  // cout << "D " << cp->D << endl;
  // cout << "y " << cp->y << endl;
  // cout << "h " << cp->h << endl;
  // cout << "partialYSigma " << cp->partialYSigma << endl;
  // cout << "partialYFlux " << cp->partialYFlux << endl;
  // cout << "partialHSigma " << cp->partialHSigma << endl;
  // cout << "partialHFlux " << cp->partialHFlux << endl;
  // cout << "partial2HSigma " << cp->partial2HSigma << endl;
  // cout << "partial2HFlux " << cp->partial2HFlux << endl;
  // cout << "partial2HSigmaFlux " << cp->partial2HSigmaFlux << endl;
  // cout << "C " << cp->C << endl;
  // cout << "D " << cp->D << endl;
  // cout << "partialWStrainPlasticStrain " << cp->partialWStrainPlasticStrain << endl;
  // cout << "plasticStrain " << cp->plasticStrain << endl;
  // cout << "internalVariables " << cp->internalVariables << endl;
  // cout << "delta plasticStrain " << cp->plasticStrain-cp->plasticStrain0 << endl;
  // cout << "delta internalVariables " << cp->internalVariables-cp->internalVariables << endl;
  PetscFunctionReturn(0);
}