SIFEL/PARMEF/SRC/pnewtonraph.cpp File Reference

#include <math.h>
#include <string.h>
#include <stdlib.h>
#include "pglobal.h"
#include "pnewtonraph.h"
#include "nssolver.h"
#include "backupsol.h"
#include "global.h"
#include "globmat.h"
#include "mechprint.h"
#include "gmatrix.h"
#include "gtopology.h"
#include "mathem.h"
#include "vector.h"
#include "matrix.h"

Go to the source code of this file.

Functions
double	par_gnewton_raphson_one_step (long lcid, nonlinman *nlman, double *fa, double *ra, double *fb, double *dr, double *fi, long istep, long &j, long li, long ini, double ierr)

---

Function Documentation

double par_gnewton_raphson_one_step	(	long	lcid,
		nonlinman *	nlman,
		double *	fa,
		double *	ra,
		double *	fb,
		double *	dr,
		double *	fi,
		long	istep,
		long &	j,
		long	li,
		long	ini,
		double	ierr
	)

Function performs calculation of one parallel load/time step of the Newton-Raphson method for the given load case. Solved equation system does not contain time variable.

Parameters:

	lcid	- load case id
	nlman	- pointer to structure conatining setup of the solver
	ra	- vector of attained displacements
	fa	- attained load vector
	fb	- residual vector or load vector increment - right-hand side vector
	dr	- vector of displacement increment - left-hand side vector
	istep	- time/load step id
	j	- inner loop step id (output parameter, set to -1 if no inner loop was performed)
	li	- initial value of time/load step id
	ierr	- required normed error of residual vector

Returns:

The function returns reached lambda parameter.

TKr, 19/02/2013 according to Tomas Koudelka

Definition at line 41 of file pnewtonraph.cpp.

References addv(), assemble_stiffness_matrix(), check_divergency(), fillm(), fillv(), Gtm, internal_forces(), probdesc::jstep, Mespr, Mp, Myrank, Ndofm, Out, Psolm, Smat, and subv().

Referenced by par_one_step().