CrackPropagation.cpp File Reference

#include <tetgen.h>
#include <moab/SpatialLocator.hpp>
#include <moab/ElemEvaluator.hpp>
#include <MoFEM.hpp>
#include <cholesky.hpp>
#include <BasicFiniteElements.hpp>
#include <Mortar.hpp>
#include <Hooke.hpp>
#include <AnalyticalFun.hpp>
#include <ComplexConstArea.hpp>
#include <ConstantArea.hpp>
#include <phg-quadrule/quad.h>
#include <NeoHookean.hpp>
#include <Smoother.hpp>
#include <VolumeLengthQuality.hpp>
#include <SurfaceSlidingConstrains.hpp>
#include <CrackFrontElement.hpp>
#include <MWLS.hpp>
#include <GriffithForceElement.hpp>
#include <CrackPropagation.hpp>
#include <petsc/private/pcimpl.h>
#include <CPSolvers.hpp>
#include <CPMeshCut.hpp>
#include <thread>

Go to the source code of this file.

Classes
struct	FractureMechanics::GetSmoothingElementsSkin
struct	FractureMechanics::OpPrint

Namespaces
	FractureMechanics

Macros
#define	SINGULARITY

Functions
static MoFEMErrorCode	snes_monitor_fields (SNES snes, PetscInt its, PetscReal fgnorm, PetscViewerAndFormat *vf)
MoFEMErrorCode	FractureMechanics::clean_pcomms (moab::Interface &moab, boost::shared_ptr< WrapMPIComm > moab_comm_wrap)
MoFEMErrorCode	FractureMechanics::broadcast_entities (moab::Interface &moab, moab::Interface &moab_tmp, const int from_proc, Range &entities, const bool adjacencies, const bool tags)
static MoFEMErrorCode	FractureMechanics::elastic_snes_rhs (SNES snes, Vec x, Vec f, void *ctx)
static MoFEMErrorCode	FractureMechanics::elastic_snes_mat (SNES snes, Vec x, Mat A, Mat B, void *ctx)
static MoFEMErrorCode	FractureMechanics::propagation_snes_rhs (SNES snes, Vec x, Vec f, void *ctx)
static MoFEMErrorCode	FractureMechanics::propagation_snes_mat (SNES snes, Vec x, Mat A, Mat B, void *ctx)

Variables
const char *	FractureMechanics::materials_list []

Macro Definition Documentation

SINGULARITY

#define SINGULARITY

Definition at line 19 of file CrackPropagation.cpp.

Function Documentation

snes_monitor_fields()

static MoFEMErrorCode snes_monitor_fields ( SNES  snes, PetscInt  its, PetscReal  fgnorm, PetscViewerAndFormat *  vf  )

static

Definition at line 70 of file CrackPropagation.cpp.

                                                                    {
  PetscViewer viewer = vf->viewer;
  Vec res;
  DM dm;
  PetscSection s;
  const PetscScalar *r;
  PetscReal *lnorms, *norms;
  PetscInt numFields, f, pStart, pEnd, p;
 
  MoFEMFunctionBegin;
  LogManager::setLog("PETSC");
  MOFEM_LOG_TAG("PETSC", "SnesProp");
 
  PetscValidHeaderSpecific(viewer, PETSC_VIEWER_CLASSID, 4);
  CHKERR SNESGetFunction(snes, &res, 0, 0);
  CHKERR SNESGetDM(snes, &dm);
  CHKERR DMGetDefaultSection(dm, &s);
  CHKERR PetscSectionGetNumFields(s, &numFields);
  CHKERR PetscSectionGetChart(s, &pStart, &pEnd);
  CHKERR PetscCalloc2(numFields, &lnorms, numFields, &norms);
  CHKERR VecGetArrayRead(res, &r);
  for (p = pStart; p < pEnd; ++p) {
    for (f = 0; f < numFields; ++f) {
      PetscInt fdof, foff, d;
 
      CHKERR PetscSectionGetFieldDof(s, p, f, &fdof);
      CHKERR PetscSectionGetFieldOffset(s, p, f, &foff);
      for (d = 0; d < fdof; ++d)
        lnorms[f] += PetscRealPart(PetscSqr(r[foff + d]));
    }
  }
  CHKERR VecRestoreArrayRead(res, &r);
  CHKERR MPIU_Allreduce(lnorms, norms, numFields, MPIU_REAL, MPIU_SUM,
                        PetscObjectComm((PetscObject)dm));
  CHKERR PetscViewerPushFormat(viewer, vf->format);
  CHKERR PetscViewerASCIIAddTab(viewer, ((PetscObject)snes)->tablevel);
  CHKERR PetscViewerASCIIPrintf(viewer, "%3D SNES Function norm %14.12e\n", its,
                                (double)fgnorm);
  for (f = 0; f < numFields; ++f) {
    const char *field_name;
    CHKERR PetscSectionGetFieldName(s, f, &field_name);
    CHKERR PetscViewerASCIIPrintf(viewer, "\t%30.30s\t%14.12e\n", field_name,
                                  (double)PetscSqrtReal(norms[f]));
  }
  CHKERR PetscViewerASCIISubtractTab(viewer, ((PetscObject)snes)->tablevel);
  CHKERR PetscViewerPopFormat(viewer);
  CHKERR PetscFree2(lnorms, norms);
 
  LogManager::setLog("PETSC");
  MoFEMFunctionReturn(0);
}