MoFEM::TimeScaleVector< SPACE_DIM > Struct Template Reference

Force scale operator for reading four columns (time and vector) More...

#include <src/boundary_conditions/ScalingMethod.hpp>

Inheritance diagram for MoFEM::TimeScaleVector< SPACE_DIM >:

Collaboration diagram for MoFEM::TimeScaleVector< SPACE_DIM >:

Public Member Functions
	TimeScaleVector (std::string name="-time_vector_file", bool error_if_file_not_given=false)
	TimeScaleVector (std::string name, int ms_id, bool error_if_file_not_given=false)
virtual FTensor::Tensor1< double, SPACE_DIM >	getVector (const double time)
virtual FTensor::Tensor1< double, SPACE_DIM >	getVectorFromData (const double time)
Public Member Functions inherited from MoFEM::ScalingMethod
virtual double	getScale (const double time)
	Get scaling at given time. More...
	ScalingMethod ()=default
virtual	~ScalingMethod ()=default

Private Member Functions
MoFEMErrorCode	timeData ()

Private Attributes
std::map< double, FTensor::Tensor1< double, SPACE_DIM > >	tSeries
int	readFile
int	debug
string	nAme
bool	errorIfFileNotGiven
PetscBool	fLg
std::function< FTensor::Tensor1< double, SPACE_DIM >(double)>	scalingMethod

Detailed Description

template<int SPACE_DIM>
struct MoFEM::TimeScaleVector< SPACE_DIM >

Force scale operator for reading four columns (time and vector)

Definition at line 97 of file ScalingMethod.hpp.

Constructor & Destructor Documentation

TimeScaleVector() [1/2]

template<int SPACE_DIM>

MoFEM::TimeScaleVector< SPACE_DIM >::TimeScaleVector	(	std::string	name = "-time_vector_file",
		bool	error_if_file_not_given = false
	)

TimeScaleVector() [2/2]

template<int SPACE_DIM>

MoFEM::TimeScaleVector< SPACE_DIM >::TimeScaleVector	(	std::string	name,
		int	ms_id,
		bool	error_if_file_not_given = false
	)

Definition at line 133 of file ScalingMethod.cpp.

                                                : readFile(0), debug(0),
      errorIfFileNotGiven(error_if_file_not_given) {
  nAme = name + std::to_string(ms_id);
  CHK_THROW_MESSAGE(timeData(), "Error in reading time data");
}

Member Function Documentation

getVector()

template<int SPACE_DIM>

FTensor::Tensor1< double, SPACE_DIM > MoFEM::TimeScaleVector< SPACE_DIM >::getVector ( const double  time )

virtual

Reimplemented in BasicBoundaryConditionsInterface::BasicBCVectorConst, and BasicBoundaryConditionsInterface::BasicBCVectorScale.

Definition at line 208 of file ScalingMethod.cpp.

                                                       {
  return scalingMethod(time);
}

getVectorFromData()

template<int SPACE_DIM>

FTensor::Tensor1< double, SPACE_DIM > MoFEM::TimeScaleVector< SPACE_DIM >::getVectorFromData ( const double  time )

virtual

Definition at line 214 of file ScalingMethod.cpp.

                                                               {
 
  // if (readFile == 0) {
  //   CHK_THROW_MESSAGE(MOFEM_OPERATION_UNSUCCESSFUL, "Data file not read");
  // }
 
  FTensor::Tensor1<double, SPACE_DIM> acc;
  FTensor::Tensor1<double, SPACE_DIM> acc0 = tSeries.begin()->second;
  FTensor::Tensor1<double, SPACE_DIM> Nf;
  FTensor::Index<'i', SPACE_DIM> i;
 
  double t0 = 0, t1, dt;
  for (auto &[ts, vec] : tSeries) {
    if (ts > time) {
      t1 = ts;
      dt = time - t0;
      acc(i) = acc0(i) + ((vec(i) - acc0(i)) / (t1 - t0)) * dt;
      break;
    }
    t0 = ts;
    acc0 = vec;
    acc = acc0;
  }
  Nf(i) = acc(i);
  return Nf;
}

timeData()

template<int SPACE_DIM>

MoFEMErrorCode MoFEM::TimeScaleVector< SPACE_DIM >::timeData

private

Definition at line 141 of file ScalingMethod.cpp.

                                                    {
 
  MoFEMFunctionBeginHot;
 
  char time_file_name[255];
  PetscBool flg = PETSC_FALSE;
  if (!nAme.empty())
    CHKERR PetscOptionsGetString(PETSC_NULL, PETSC_NULL, nAme.c_str(),
                                 time_file_name, 255, &flg);
 
  if (!flg) {
    if (errorIfFileNotGiven)
      SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "*** ERROR %s (time_data FILE NEEDED)", nAme.c_str());
    MoFEMFunctionReturnHot(0);
  }
 
  FILE *time_data = fopen(time_file_name, "r");
  if (time_data == NULL) {
    SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
             "*** ERROR data file < %s > open unsuccessful", time_file_name);
  }
  double no1 = 0.0;
  FTensor::Index<'i', SPACE_DIM> i;
  FTensor::Tensor1<double, SPACE_DIM> no2;
  no2(i) = 0.;
  tSeries[no1] = no2;
  while (!feof(time_data)) {
    int n =
        fscanf(time_data, "%lf %lf %lf %lf", &no1, &no2(0), &no2(1), &no2(2));
    if (n < 0) {
      fgetc(time_data);
      continue;
    }
    if (n != 4) {
      SETERRQ1(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
               "*** ERROR read data file error (check input time data file) "
               "{ n = %d }",
               n);
    }
    tSeries[no1] = no2;
  }
  int r = fclose(time_data);
  if (debug) {
 
    for (auto &[ts, vec] : tSeries) {
      PetscPrintf(PETSC_COMM_WORLD,
                  "** read vector %3.2e time %3.2e %3.2e %3.2e\n",
                  ts, vec(0), vec(1), vec(2));
    }
  }
  if (r != 0) {
    SETERRQ(PETSC_COMM_SELF, 1, "*** ERROR file close unsuccessful");
  }
  readFile = 1;
 
  if (readFile == 1)
    scalingMethod = [this](double time) {
      return this->getVectorFromData(time);
    };
 
  MoFEMFunctionReturnHot(0);
}

Member Data Documentation

debug

template<int SPACE_DIM>

int MoFEM::TimeScaleVector< SPACE_DIM >::debug

private

Definition at line 113 of file ScalingMethod.hpp.

errorIfFileNotGiven

template<int SPACE_DIM>

bool MoFEM::TimeScaleVector< SPACE_DIM >::errorIfFileNotGiven

private

Definition at line 115 of file ScalingMethod.hpp.

fLg

template<int SPACE_DIM>

PetscBool MoFEM::TimeScaleVector< SPACE_DIM >::fLg

private

Definition at line 117 of file ScalingMethod.hpp.

nAme

template<int SPACE_DIM>

string MoFEM::TimeScaleVector< SPACE_DIM >::nAme

private

Definition at line 114 of file ScalingMethod.hpp.

readFile

template<int SPACE_DIM>

int MoFEM::TimeScaleVector< SPACE_DIM >::readFile

private

Definition at line 113 of file ScalingMethod.hpp.

scalingMethod

template<int SPACE_DIM>

std::function<FTensor::Tensor1<double, SPACE_DIM>(double)> MoFEM::TimeScaleVector< SPACE_DIM >::scalingMethod

private

Initial value:

=
      [this](double time) {
        FTensor::Index<'i', SPACE_DIM> i;
        FTensor::Tensor1<double, SPACE_DIM> s;
        s(i) = time;
        return s;
      }

Definition at line 118 of file ScalingMethod.hpp.

tSeries

template<int SPACE_DIM>

std::map<double, FTensor::Tensor1<double, SPACE_DIM> > MoFEM::TimeScaleVector< SPACE_DIM >::tSeries

private

Definition at line 112 of file ScalingMethod.hpp.

---

The documentation for this struct was generated from the following files:

• src/boundary_conditions/ScalingMethod.hpp
• src/boundary_conditions/impl/ScalingMethod.cpp