gmNfFileFormat.h

Go to the documentation of this file.

/************************************************************************
**
** Copyright (C) 2014 by Carlos Augusto Teixera Mendes
** All rights reserved.
**
** This file is part of the "GeMA" software. It's use should respect
** the terms in the license agreement that can be found together
** with this source code.
** It is provided AS IS, with NO WARRANTY OF ANY KIND,
** INCLUDING THE WARRANTY OF DESIGN, MERCHANTABILITY AND FITNESS FOR
** A PARTICULAR PURPOSE.
**
************************************************************************/

#ifndef _GM_NF_FILE_FORMAT_H_
#define _GM_NF_FILE_FORMAT_H_

#include "gmFileFormat.h"

#include "gmValueInfo.h"
#include <nfr.h>

class GmNfFileFormat : public GmFileFormat
{
public:
  // See comments on the base class
  virtual const char* formatName() const { return "Neutral file"; }

  // See comments on the base class
  virtual bool supportsMultipleIterations() const { return true; }

  // See comments on the base class
  virtual bool supportsSplitVectors() const { return true; }

  // See comments on the base class
  virtual bool acceptsNodeDimension(int ndim) const { return ndim == 2 || ndim == 3; }

  // See comments on the base class
  virtual bool acceptsDataType(const GmValueInfo* info, int nodeDim, int dimFilter = -1) const
  {
    // The neutral file format accepts scalar, vector and tensor values
    Q_UNUSED(info); Q_UNUSED(nodeDim); Q_UNUSED(dimFilter);
    return true;
  }

  // See comments on the base class
  virtual bool acceptsMultipleGaussProfiles() const { return false; }

  virtual bool acceptsGaussRuleForElement(GmCellType type, const GmIntegrationRule* ir) const;

  // See comments on the base class
  virtual bool acceptsDiscontinuitySet(GmDiscontinuitySet::DiscontinuityType type) const { Q_UNUSED(type); return false; }

  // TODO
  virtual int acceptsCellType(GmCellType type) const { Q_UNUSED(type); return 0; }
  virtual const QVector<int>& cellNumbering(GmCellType type) const { Q_UNUSED(type); static QVector<int> dummy;  return dummy; }
  virtual int fileCellType(GmCellType type) const { Q_UNUSED(type); return 0; }
  virtual const char* fileCellTypeName(GmCellType type) const { Q_UNUSED(type); return NULL; }

  // See comments on the base class
  virtual bool supportsMeshChanges() const { return false; }

protected:
  // Supported element types and file section names
  struct SupportedTypesT
  {
    GmCellType  type;
    const char* topoSection;
    int         ndummy1;       // Number of dummy values after the material label and before nodes ids
    bool        integrationId; // Should we add an integration aid after the dummy values?
    int         ndummy2;       // Number of dummy values after nodes ids
  };

  static const SupportedTypesT supportedTypes[];
  static const int supportedTypesSize;

  struct nfrType {
    NfrResultType  nfType;
    GmValueSetKind gmKind;
    GmDimType      gmType;
    int            nfSize;
  };

  struct nfrElement {
    NfrElementType    nfEleType;
    NfrQuadratureType nfQuadType;
    double            nfQuadNatCrd[200];
  };

  static const nfrType    nfrTypeTable[];
  static const nfrElement nfrElementTable[];

  int nf2gmCellNodeIndex(GmCellType type, int nfId);
  int nf2gmCellGaussIndex(GmCellType type, int ngp, int nfId);

  unsigned int offset(unsigned int n) { return n == 1 ? 0 : n-1 + offset(n-1); }
};


#endif // _GM_NF_FILE_FORMAT_H_