alg.h

Go to the documentation of this file.

/*
%M This module contains the data structure definitions for the solution
   algorithm.
%a Joao Luiz Elias Campos.
%d September 3rd, 1998.
%r $Id: alg.h,v 1.1 2004/06/22 05:29:58 joaoluiz Exp $
%w (C) COPYRIGHT 1995-1996, Eduardo Nobre Lages.
   (C) COPYRIGHT 1997-1998, Joao Luiz Elias Campos.
   All Rights Reserved
   Duplication of this program or any part thereof without the express
   written consent of the author is prohibited.
 *
 *   Modified:    08-Aug-05    Alexandre A. Del Savio
 *     Adicionado "_count_it".
 *
 *   Modified:    13-Sep-05    Alexandre A. Del Savio
 *     Criada a funcao AlgReset.
 *
 *   Modified:    08-Dec-05    Juan Pablo Ibañez
 *     Acrescentou-se o parâmetro numsteptol na struct config.
 *     Este parâmetro define o tamanho do vetor de erros usado
 *     para controle de convergencia.
 *
 *   Modified   16-Dez-05    Juan Pablo Ibañez
 *     Acrescentou-se o parâmetro delta na função _count_it.
 *
 */

#ifndef _ALG_H
#define _ALG_H

#include "rui.h"

/*
%T Coupled state definition
*/
typedef enum _coupstate
{
 UNCOUPLED = 0,
 COUPLED   = 1
} eCoupState;

/*
%T Configuration parameters
*/
typedef struct _config
{
 int    max_iter;
 int    nonlinear;
 int    print_step;
 int    draw_step;
 int    num_step;
 int    num_time_step;
 double tolerance;
 int    numsteptol; 
 double dtfrac;
 double totalTime;
 double timeStep;
 int    num_load_step;
 double loadfactor;
 int    algtype;
 int    numgaussresults;
 char   **gaussresults;
 int    solver;
} sConfig;

void    (*_count_it)(int, double, int); 
void    (*_fstop)(int *); 

/* Configuration parameters
 */
extern sConfig  Config;

extern void (*_count_it)(int, double, int);
extern void (*_fstop)(int *);

/* Flag to see if is a coupled problem
 */
extern eCoupState State;

/* Macro definition
 */
#define IsGeomNonLinear  (Config.nonlinear == 1)
#define IsCoupled        (State == COUPLED)

/* Public functions
 */
int  DRSolver( UI_State *, double *, double *, double *, double *, int );
int  IMPLINEARSolver( UI_State *, double *, double *, double *, double *, int );
int  IMPBFGSSolver( UI_State *, double *, double *, double *, double *, int );
int  IMPNRMSolver( UI_State *, double *, double *, double *, double *, int );
int  HYBRIDSolver(UI_State *R, double *FVector, double *UVector, double *VVector,
               double *MVector, int);
void AlgReset( void );
double ComputeError  ( double *, double *, double);
void _dUVector( double *, double *, double *);
int  NeedIteration ( double *, int *, int );
#endif

---