geomcomp.h

Go to the documentation of this file.

#ifndef MIDAS_GEOMCOMP_H
#define MIDAS_GEOMCOMP_H

#include "geometry.h"
#include "attribute.h"

#include "gelib.h"
#include "arrays.h"


namespace midaspace {

class GeometryComponent : public GeometrySubject
{
 private:
  // Property HELP
  // =============
  // - p. slouzi pro prirazeni vlastnosti nejakemu komponentu (uzlu/elementu).
  // - jeden komponent muze mit zadnou, jednu nebo vice hodnot property, nevim ale, jak to bude s prednosti dedeni, zatim pouzivam jednu
  // - property je typu integer
  // - p. by mela byt kladna nenulova, pak zarucuji funkcnost
  //   - muze byt i zaporna (nedoporucuji, neni duvod) i nulova (nedoporucuji, t3d pri neurcene prop priradi prop 0)
  // - zadani jako mesh: zadat muzu property ve VTK na puvodni siti
  // - zadani jako model: property je ID (lidske cislovani od 1) rodicovskeho komponentu z modelu
  //   ==================
  //   UZEL - main  prop. - vertexu            prop z modelu, na kterym uzel lezi; 1 i vice hodnot, pokud T3d slouci 2 vertexy modelu do jednoho uzlu meshe
  //          other prop. - edge, face, volume prop z modelu, na kterym uzel lezi; 1 i vice hodnot zcela bezne
  //        IMPLEMENTACE mprop je zde, tj. geomcomp.h; oprop je v point.h
  //        a) rucni prirazeni - mprop je jedna, podle ni se priradi informace z ctrl.txt; oprop se ignorujou
  //        b) meshovani v midasu, tj. existuje "master model geometry"
  //           - attributes se budou dedit prioritne ze vsech property, zatim se dedi jen podepreni:
  //             if (jedina mpop)     tak se podepreni zdedi z vertexu
  //             else if (edge oprop) tak se podepreni zdedi z z edge (tj. edge existujici v modelu); pokud budou 2 nenulove edge oprop, tak je to chyba, pak musi byt nenulova mprop
  //             else if (face oprop) tak se podepreni zdedi z z face (tj. face existujici v modelu); pokud budou 2 nenulove face oprop, tak je to chyba, pak musi byt nenulove neco predchazejici
  //
  //   ELEMENT - main  property - asi jen jedna (pokud vic tak je to divny, proc by t3d slucovalo elementy)
  //           - other property - zadna
  //        IMPLEMENTACE
  //        a) rucni prirazeni - mprop je jedna, podle ni se priradi informace z ctrl.txt
  //        b) meshovani v midasu, tj. existuje "master model geometry" a attributes se zdedi
  //
  //   FACEDGE - other property - zadna
  //           - main  property - muze mit i vice property, napr. kdyz mam 2 uzly a mezi nemi lezi 2 pruty, tak ja z nich udelam jednu edge
  //        IMPLEMENTACE
  //        a) rucni prirazeni
  //           - edgi priradim property pouze v t3d (mozna v jktk, ale to zatim nepodporuju)
  //           - mprop je jedna, tu kterou jsme uvedl u "curve", pokud jich bude vic, tak tak error = nepodporuju to
  //             pak podle ni se priradi informace z ctrl.txt primo edgi
  //        b) meshovani v midasu, tj. existuje "master model geometry" a attributes se zdedi,, zatim se dedi jen zatizeni, na edgi hranove na face plosne
  //           - mprop  = 0 - zadnej master, nic se nededi
  //             mprop != 0 - dedi se jen z virtualniho mastra, skutecnej master zatizeni ma
  //                          pokud jsou 2 virtualni tak error
  
  // edge zatizeni dam na z attribudu elementu na attribut edge
  
  // dodelam gravity
  
 protected:
  Lvctr mproperty;       // main property of component (NULL if nonactivated)
  long origid;            // original ID from input file
  bool fordelete;         // flag - under sentence of deletion
  
  Attributes *attributes;   // attributes of geometry component
  
 public:
  GeometryComponent (const Geometry *owner, long gid, long oid, long prop) : GeometrySubject (gid, owner) {
    if (prop) this->set_mprop(prop);
    origid   = oid;
    fordelete= false;
    attributes = NULL;
  }
  GeometryComponent (const GeometryComponent *src) : GeometrySubject (src->give_ID(), src->Geom) {
    mproperty = src->mproperty;
    origid    = src->origid;
    fordelete = src->fordelete;
    attributes = NULL;
  }
  virtual ~GeometryComponent() { delete attributes; }
  
  virtual void initialize       (void)       { GeometrySubject :: initialize ();       attributes->initialize();       }
  virtual void finitialize      (void)       { GeometrySubject :: finitialize ();      attributes->finitialize();      }
  virtual void checkConsistency (void) const { GeometrySubject :: checkConsistency();  attributes->checkConsistency();
#ifdef DEBUG
    if (mproperty.give_number_of_zeros()) _errorr("zero property");
#endif
  } // if (attributes == NULL)  errol; 
  //void Element :: initialize       (void)       
  //void Element :: finitialize      (void)       
  //void Element :: checkConsistency (void) const 
  
  
# ifndef DEBUG
  virtual void set_mprop (long  val) { if (mproperty.give_size() == 0) { mproperty.resize_ignore_vals(1);  mproperty[0] = val;               }  else _errorr("main property set already"); }
  void       reset_mprop (long  val) { if (mproperty.give_size() == 0)   mproperty.resize_ignore_vals(1);  mproperty[0] = val;                                                            }
  bool    checkset_mprop (long  val) { if (mproperty.give_size() == 0) { mproperty.resize_ignore_vals(1);  mproperty[0] = val;  return true; }  else if (mproperty[0] != val) _errorr("different mproperty[0]");  return false; }
# else
  virtual void set_mprop (long  val) { if (val == 0) errol;  if (mproperty.give_size() == 0) { mproperty.resize_ignore_vals(1);  mproperty[0] = val;               }  else _errorr("main property set already"); }
  void       reset_mprop (long  val) { if (val == 0) errol;  if (mproperty.give_size() == 0)   mproperty.resize_ignore_vals(1);  mproperty[0] = val;                                                            }
  bool    checkset_mprop (long  val) { if (val == 0) errol;  if (mproperty.give_size() == 0) { mproperty.resize_ignore_vals(1);  mproperty[0] = val;  return true; }  else if (mproperty[0] != val) _errorr("different mproperty[0]");  return false; }
# endif
  //
  void set_origID       (long  val) {    origid = val; }
  void set_delete_flag  (bool  val) { fordelete = val; }
  //
  virtual void set_model_prop (long val, const Model *model, bool flag=false) { errol; }
  
  long         give_mproperty     () const {
    if (!mproperty.give_size())
      _errorr ("Property not initiated");
    return mproperty[0];
  }
  const Lvctr* give_mproperty_ptr  (void) const { return &mproperty;  }
  long         give_mpropertyORzero(void) const { return (mproperty.give_size() ? mproperty[0] : 0); }
  long         give_mproperty_cnt  (void) const { return  mproperty.give_size(); }
  bool         has_mproperty  (long prop) const { return (mproperty.give_size() ? (mproperty[0] == prop ? true : false) : false); }
  //
  long  give_origID      () const { return origid;    }
  bool  give_delete_flag () const { return fordelete; }
        Attributes *give_attributes (void)       { return attributes; }
  const Attributes *give_attributes (void) const { return attributes; }
  
  
  bool Parallel           (void) const { return ((Mesh*)Geom)->give_Parallel();          }
  long NumDomains         (void) const { return ((Mesh*)Geom)->give_NumDomains();        }
  void check_connectivity (void) const {
    if (!((Mesh*)Geom)->connectivity_is_assembled())
      _errorr ("Globla connectivity is not assembled");
  }
  void mg_connectivity_assembling (void) { ((Geometry*)Geom)->connectivity_assembling(); }
  
  virtual void print_row (FILE *stream, femFileFormat fff, bool endline=true, long did=0) const = 0;
  
  
};

} // namespace midaspace

#endif // MIDAS_GEOMCOMP_H