/*Translated by FOR_C, v3.4.2 (-), on 07/09/115 at 08:30:09 */
/*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 "sdasgh.h"
#include <stdlib.h>
void /*FUNCTION*/ sdasgh(
float est,
float estold,
float *ratio,
float h,
float *hnew,
float tk)
{
	float a2, b1, b2, dkappa, rho;
 
	/* Copyright (c) 2006, Math a la Carte, Inc.
	 *>> 2003-03-06 sdasgh Hanson installed Soderlind smoothing logic
	 *     This is Soderlind's 'H211B4' smoothing method.
	 *     It is  his preferred choice and it is our only choice.
	 *--S replaces "?": ?dasgh */
	b1 = 0.25e0/tk;
	b2 = 0.25e0/tk;
	a2 = 0.25e0;
	dkappa = 1.e0;
 
	*ratio = powf(2.0e0*est + 1.0e-8,-b1)*powf(2.0e0*estold + 1.0e-8,-b2)*
	 powf(*ratio,-a2);
 
	rho = 1.0e0 + dkappa*atanf( (*ratio - 1.0e0)/dkappa );
 
	*hnew = h*rho;
 
	return;
} /* end of function */