solverInterface.h

Go to the documentation of this file.

/* SPDX-FileCopyrightText: 2023 SAP SE
 *
 * SPDX-License-Identifier: Apache-2.0
 *
 * This file is part of FEDEM - https://openfedem.org
 */
#ifndef _SOLVER_INTERFACE_H
#define _SOLVER_INTERFACE_H
 
#include <stddef.h>
 
#ifdef __cplusplus
extern "C" {
#else
#include <stdbool.h>
#endif
 
  int solverInit(int argc, char** argv, const char* cfsi = NULL,
                 const double* stateData = NULL, const int ndat = 0,
                 const double* gagesData = NULL, const int ngda = 0,
                 const double* extInput  = NULL, int* nxinp = NULL);
 
  int restartFromState(const double* stateData, const int ndat,
                       const int writeToRDB = 2);
 
  bool solveWindow(int nStep, int nInc, int nIn, int nOut, const int* fOut,
                   const double* times, const double* inputs, double* outputs,
                   int nDat, double* stateData, int* ierr);
 
  int getStateSize();
 
  int getTransformationStateSize();
 
  int getPartDeformationStateSize(int bid);
 
  int getPartStressStateSize(int bid);
 
  int getGagesSize();
 
  bool saveState(double* stateData, const int ndat);
 
  bool saveTransformationState(double* stateData, const int ndat);
 
  bool saveGages(double* stateData, const int ndat);
 
  bool getStrainsFromDisp(const double* disp, const int* gageIDs, double* eps,
                          const int ndof, const int ng);
 
  bool getBeamForcesFromDisp(const double* disp, const int* beamIDs,
                             double* forces, const int ndof, const int nBeam);
 
  bool getRelDisp(const double* disp, const int* Ids, double* relDis,
                  const int ndof, const int nId);
 
  bool getRespVars(const double* disp, const int* Ids, double* resVar,
                   const int ndof, const int nId);
 
  bool savePartDeformationState(int bid, double* stateData, const int ndat);
 
  bool savePartStressState(int bid, double* stateData, const int ndat);
 
  bool solveNext(int* ierr);
 
  bool startStep(int* ierr);
 
  bool solveIteration(int* ierr, bool all = false);
 
  bool solveEigenModes(const int nmod, double* eval, double* evec,
                       bool dofOrder, int useLaPack, int* ierr);
 
  bool solveInverse(const double* x,
                    const int* xeqs, const int* feqs,
                    int ndis, int nfrs, int* ierr);
 
  int solverDone(bool removeSingletons = true);
 
  void solverClose();
 
  bool setExtFunc(int funcId, double value);
 
  double getTime(bool nextStep, int* ierr);
 
  bool setTime(double nextTime);
 
  inline double getCurrentTime(int* ierr) { return getTime(false,ierr); }
 
  inline double getNextTime(int* ierr) { return getTime(true,ierr); }
 
  double evalFunc(int uid, const char* tag = NULL,
                  double x = -1.0, int* ierr = NULL);
 
  double getFunc(int uid, int* ie = NULL) { return evalFunc(uid,NULL,-1.0,ie); }
 
  int getFuncId(const char* tag);
 
  int getEquations(int bid, int* meqn);
 
  int getStateVar(int bid, double* var);
 
  int getSystemSize(bool dofs = false);
 
  bool getSystemMatrix(double* Smat, int iMat = 0);
 
  bool getNewtonMatrix(double* Nmat);
 
  bool getStiffnessMatrix(double* Kmat);
 
  bool getMassMatrix(double* Mmat);
 
  bool getDampingMatrix(double* Cmat);
 
  bool getElementStiffnessMatrix(double* Kmat, int bid);
 
  bool getRhsVector(double* Rvec, int iVec = 0);
 
  bool setRhsVector(const double* Rvec);
 
  bool addRhsVector(const double* Rvec);
 
  bool getJointSprCoeff(double* sprCoeff, int bid);
 
#ifdef __cplusplus
}
#endif
 
#endif