/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:30:04 */
/*FOR_C Options SET: ftn=u io=c no=p op=aimnv s=dbov str=l x=f - prototypes */
#include <math.h>
#include "fcrt.h"
#include "ddasgh.h"
#include <stdlib.h>
void /*FUNCTION*/ ddasgh(
double est,
double estold,
double *ratio,
double h,
double *hnew,
double tk)
{
	double a2, b1, b2, dkappa, rho;
 
	/* Copyright (c) 2006, Math a la Carte, Inc.
	 *>> 2003-03-06 ddasgh Hanson installed Soderlind smoothing logic
	 *     This is Soderlind's 'H211B4' smoothing method.
	 *     It is  his preferred choice and it is our only choice.
	 *--D replaces "?": ?dasgh */
	b1 = 0.25e0/tk;
	b2 = 0.25e0/tk;
	a2 = 0.25e0;
	dkappa = 1.e0;
 
	*ratio = pow(2.0e0*est + 1.0e-8,-b1)*pow(2.0e0*estold + 1.0e-8,-b2)*
	 pow(*ratio,-a2);
 
	rho = 1.0e0 + dkappa*atan( (*ratio - 1.0e0)/dkappa );
 
	*hnew = h*rho;
 
	return;
} /* end of function */