SmoothMultiBoundingBoxBaseIC is the base class for IsolatedBoundingBoxIC and NestedBoundingBoxIC. More...

#include <SmoothMultiBoundingBoxBaseIC.h>

Inheritance diagram for SmoothMultiBoundingBoxBaseIC:

Public Member Functions
	SmoothMultiBoundingBoxBaseIC (const InputParameters &parameters)

Real	value (const Point &p)

Protected Attributes
const unsigned int	_nbox
	number of boxes More...

const Real	_int_width
	value of interfacial width More...

const unsigned int	_dim
	dimensionality of the mesh More...

std::vector< Real >	_inside
	values inside the boxes More...


const std::vector< Point >	_c1
	lists of opposite corners More...

const std::vector< Point >	_c2

Private Attributes
const Real	_outside

Detailed Description

SmoothMultiBoundingBoxBaseIC is the base class for IsolatedBoundingBoxIC and NestedBoundingBoxIC.

The boxes can finite interface width.

Definition at line 24 of file SmoothMultiBoundingBoxBaseIC.h.

Constructor & Destructor Documentation

◆ SmoothMultiBoundingBoxBaseIC()

SmoothMultiBoundingBoxBaseIC::SmoothMultiBoundingBoxBaseIC ( const InputParameters & parameters )

Definition at line 37 of file SmoothMultiBoundingBoxBaseIC.C.

   : InitialCondition(parameters),
     _outside(getParam<Real>("outside")),
     _c1(getParam<std::vector<Point>>("smaller_coordinate_corners")),
     _c2(getParam<std::vector<Point>>("larger_coordinate_corners")),
     _nbox(_c1.size()),
     _int_width(getParam<Real>("int_width")),
     _dim(_fe_problem.mesh().dimension()),
     _inside(getParam<std::vector<Real>>("inside"))
 {
 }

Member Function Documentation

◆ value()

Real SmoothMultiBoundingBoxBaseIC::value ( const Point & p )

Definition at line 50 of file SmoothMultiBoundingBoxBaseIC.C.

 {
   Real value = _outside;
  
   // if "inside" vector only has size 1, all the boxes have the same inside value
   if (_inside.size() == 1)
   {
     _inside.assign(_nbox, _inside[0]);
   }
  
   // make sure inputs are the same length
   if (_c2.size() != _nbox || _inside.size() != _nbox)
     paramError("vector inputs must all be the same size");
  
   if (_int_width < 0.0)
     paramError("'int_width' should be non-negative");
  
   if (_int_width == 0.0)
   {
     for (unsigned int b = 0; b < _nbox; ++b)
     {
       for (unsigned int i = 0; i < _dim; ++i)
         if (_c1[b](i) < _c2[b](i) && p(i) >= _c1[b](i) && p(i) <= _c2[b](i))
         {
           if (i != _dim - 1)
             continue;
           value = _inside[b];
           break;
         }
         else
           break;
     }
   }
   return value;
 }