mnode.cpp

Go to the documentation of this file.

//********************************************************************************************************
// code:                          ###  ###  #####   ####  ##  ##
//                       ##  ##   ########  ##     ##  ## ##  ##
//                       ##  ##   ## ## ##  ###    ##     ######
//                       ##  ##   ##    ##  ##     ##  ## ##  ##
//                       #####    ##    ##  #####   ####  ##  ##
//                       ##
//
// name:           point.cpp
// author(s):      Ladislav Svoboda
// language:       C, C++
// license:        This program is free software; you can redistribute it and/or modify
//                 it under the terms of the GNU Lesser General Public License as published by
//                 the Free Software Foundation; either version 2 of the License, or
//                 (at your option) any later version.
//
//                 This program is distributed in the hope that it will be useful,
//                 but WITHOUT ANY WARRANTY; without even the implied warranty of
//                 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//                 GNU Lesser General Public License for more details.
//
//                 You should have received a copy of the GNU Lesser General Public License
//                 along with this program; if not, write to the Free Software
//                 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
//********************************************************************************************************

#include "mnode.h"

#include "problem.h"
#include "mesh.h"


namespace mumech {

mNode :: mNode (long i, const Mesh *m)
{
  id = i;
  M = m;
  
  //coords = new double[M->P->give_VECT_RANGE()];
  
  displc = strain = stress = NULL;
}
mNode :: mNode (long i, const Mesh *m, const mNode *src)
{
  id = i;
  M = m;

  coords.copy(&src->coords);

  displc = strain = stress = NULL;
  if (src->displc != NULL  ||  src->strain != NULL  ||  src->stress != NULL)  _errorr("results data not supported in copy constructor");
}

mNode :: ~mNode ()
{
  int nP  = M->npa;
  int nLC = M->nLC;
  
  DeleteArray3D(displc, nP, nLC);
  DeleteArray3D(strain, nP, nLC);
  DeleteArray3D(stress, nP, nLC);
}

void mNode :: allocate_fields (int pid)
{
  int nLC = M->nLC;
  
  if (M->no_dspl_rf) { if (!displc)  displc = Allocate1Ddpp (M->npa);  if (M->no_dspl_rf[pid]) { if (!displc[pid])  displc[pid] = Allocate1Ddp(nLC);  for (int i=0; i<nLC; i++) { if (M->no_dspl_rf[pid][i] && displc[pid][i] == NULL)  displc[pid][i] = new double[M->P[pid]->ndisplc_one()]; }  }  }
  if (M->no_strn_rf) { if (!strain)  strain = Allocate1Ddpp (M->npa);  if (M->no_strn_rf[pid]) { if (!strain[pid])  strain[pid] = Allocate1Ddp(nLC);  for (int i=0; i<nLC; i++) { if (M->no_strn_rf[pid][i] && strain[pid][i] == NULL)  strain[pid][i] = new double[M->P[pid]->nstrain_one()]; }  }  }
  if (M->no_strs_rf) { if (!stress)  stress = Allocate1Ddpp (M->npa);  if (M->no_strs_rf[pid]) { if (!stress[pid])  stress[pid] = Allocate1Ddp(nLC);  for (int i=0; i<nLC; i++) { if (M->no_strs_rf[pid][i] && stress[pid][i] == NULL)  stress[pid][i] = new double[M->P[pid]->nstrain_one()]; }  }  }
}


} // end of namespace mumech

/*end of file*/