rezone.c

Go to the documentation of this file.

/*
%M This modules contains the driver methods to solve a problem using the 
   dynamic relax algorithm with rezone steps.
%a Joao Luiz Elias Campos.
%d July 27th, 1999.
%r $Id: rezone.c,v 1.1 2004/06/22 05:29:59 joaoluiz Exp $
%w (C) COPYRIGHT 1995-1996, Eduardo Nobre Lages.
   (C) COPYRIGHT 1997-1999, 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: 28-Apr-2005    Alexandre A. Del Savio
*    Foram substituídas todas as alocações dinâmicas feitas com malloc 
*    por calloc.
*   
*  Modified: 23-Fev-2006    Juan Pablo Ibañez
*    Modificada a funcão RezoneAnalysis, que passa a retornar un valor int.   
*   
*/

/*
** ------------------------------------------------------------------------
** Global variables and symbols:
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

#include "drv.h"
#include "fem.h"
#include "load.h"
#include "elm.h"
#include "node.h"
#include "alg.h"
#include "rio.h"
#include "nfi.h"
#include "load.h"
#include "xgplib.h"


/*
** ------------------------------------------------------------------------
** Local variables and symbols:
*/

/*
** ------------------------------------------------------------------------
** Local functions:
*/
static void RezoneNew        ( sDriver ** );
static void RezoneFree       ( sDriver * );
static int  RezoneAnalysis   (UI_State *);
static int  RezonePrintResult( int, double *, double * );


/*
%F
*/
static void RezoneNew( sDriver **drv )
{
/* Get memory for the driver descriptor
 */
 (*drv) = (sDriver *)calloc(1, sizeof(sDriver));

/* Fill up driver descriptor
 */
 (*drv)->type = REZONE;

} /* End of RezoneNew */


/*
%F
*/
static void RezoneFree( sDriver *drv )
{
/* Reset driver data
 */
 drv->data = 0L;

} /* End of RezoneFree */


/*
%F This function process the probleam analysis using a dynamic relax 
   algorithm.
*/
static int RezoneAnalysis(UI_State *R)
{
 int     i, esc = 0;
 double *FVector = 0L;
 double *UVector = 0L;
 double *VVector = 0L;
 double *MVector = 0L;
 
/* Get memory for the vectors used during the analysis
 */
 FVector = (double *)calloc( NDof*NumNodes, sizeof(double) );
 UVector = (double *)calloc( NDof*NumNodes, sizeof(double) );
 VVector = (double *)calloc( NDof*NumNodes, sizeof(double) );
 MVector = (double *)calloc( NDof*NumNodes, sizeof(double) );

/* Apply prescribed values
 */
 PrescribedValues( );

/* Start graphic module
 */
 XGPBegin( );

/* Solve the problem
 */
 while( esc < (NumRezones + 1) )
 {
 /* Print feedback message
  */
  if( esc > 0 )
   printf( "\tProcess rezone step..................: %d\n", esc );
  else
   printf( "\tProcess first step without rezone....: %d\n", esc );
  fflush( stdout );

 /* Process rezoning
  */
  if( esc > 0 )
   DoRezone( esc, FVector );

 /* Call the dynamic relaxation algorithm for each rezone step
  */
  if( !DRSolver(R, FVector, UVector, VVector, MVector, 0) ) return 0;

 /* Write results on a temporary file
  */
  if( !IoStartSave( ) ) return 0;

  fwrite( UVector, sizeof(double), NDof*NumNodes, ndlr );
  fflush( ndlr );

  for( i = 0; i < NumElements; i++ )
   ElmWriteStress( ElmList[i], elmr, UVector, VVector );
  fflush( elmr );

 /* Increase the rezone step
  */
  esc++;
 }

/* Finish graphic module
 */
 XGPEnd( );

/* Release memory
 */
 free( FVector );
 free( UVector );
 free( VVector );
 free( MVector );

 return 1;

} /* End of RezoneAnalysis */


/*
%F
*/
static int RezonePrintResult( int iteration, double *U, double *V )
{
 return 1;
} /* End of RezonePrintResult */


/*
** ------------------------------------------------------------------------
** Public functions:
*/

/*
%F This function initializes the sub-class REZONE.
*/
void RezoneInit( void );
void RezoneInit( void )
{
/* Initialize REZONE driver sub-class
 */
 DrvClass[REZONE].new      = RezoneNew;
 DrvClass[REZONE].free     = RezoneFree;
 DrvClass[REZONE].analysis = RezoneAnalysis;
 DrvClass[REZONE].printres = RezonePrintResult;

} /* End of RezoneInit */


/* =======================================================  End of File  */

---