userl_solvparam.for

来自「两维弹性边界元程序」· FOR 代码 · 共 28 行

!For User to adjust 
! PARAMETERS used in the ITERATIVE SOLVER for best efficiency
!Alpha, AlphaC: between 0.3 and 1.2||Normally, a smaller value 
! leads to more stable iteration, but less convergence speed. 
! The optimum value varies with configuration including geometry 
! and loading condition.
!UErrAllow, PErrAllow: If differences of displacement and stress in consective
! iteration steps are both smaller than these parameters respectively, 
! the iteration process is stopped.
!IteratShow: every this number of iteration steps is shown the errors 
! in the solution so that one may see the process is diverging or converging
*********************************************************
	Subroutine UserSolvParam(
	1	Alpha,AlphaC,UErrAllow,PErrAllow,IteratShow)
     	Implicit Real*8 (A-H,O-Z)
	Implicit Integer (I-N)
	Common/Example/NthExample
*============================================================
!These values are used for all the examples.
	Alpha=0.9d0
	AlphaC=0.7
	UErrAllow=1.d-7
	PErrAllow=1.d-7
	IteratShow=50
*============================================================
	return
	end