Personal tools
You are here: Home Projects C++ Cfront releases Release 1.0 Source cfront lib complex abs.c
Document Actions

abs.c

by Paul McJones last modified 2007-02-02 09:32

Click here to get the file

Size 1 kB - File type text/x-csrc

File contents

/* @(#) abs.c 1.2 1/27/86 17:47:28 */
/*ident	"@(#)cfront:lib/complex/abs.c	1.2"*/
#include "complex.h"
#include "const.h"

double
abs(complex z)
{
	register double temp;
	register double r = z.re;
	register double i = z.im;
	
	if (r < 0) r = -r;
	if (i == 0) return r;

	if (i < 0) i = -i;
	if (r == 0) return i;

	if (i > r) { temp = r; r = i; i = temp; }
	temp = i/r;
	temp = r*sqrt(1.0 + temp*temp);  /*overflow!!*/
	return temp;
}

double
norm(complex z)
{
	register double r = z.re;
	register double i = z.im;

#define	SAFE 6.521908912666391000e+18  /* 0.5*sqrt(GREATEST)*/

	if (r<SAFE && i<SAFE) return r*r+i*i;
	return GREATEST;
}
« November 2014 »
Su Mo Tu We Th Fr Sa
1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30
 

Powered by Plone CMS, the Open Source Content Management System

This site conforms to the following standards: