psmile_ddadd.c File Reference

#include <stddef.h>
#include "PSMILe_f2c.h"

Include dependency graph for psmile_ddadd.c:

Go to the source code of this file.

Functions/Subroutines
void	psmile_ddadd (double *dda, double *ddb)

---

Function Documentation

void psmile_ddadd	(	double *	dda,
		double *	ddb
	)

Definition at line 31 of file psmile_ddadd.c.

Referenced by psmile_ddadd_mpi_callback(), psmile_get_field_dble(), psmile_get_field_real(), psmile_global_sum_compute_dble(), psmile_put_field_dble(), and psmile_put_field_real().

{
   volatile double e, t1, t2, t2_1, t2_2, t2_3;
/* Compute dda + ddb using Knuth's trick.*/
   t1 = dda[0] + ddb[0];
   e = t1 - dda[0];
   t2_1 = t1 - e;
   t2_2 = dda[0] - t2_1;
   t2_3 = ddb[0] - e;
   t2   = t2_3 + t2_2;
   t2  += dda[1];
   t2  += ddb[1];

/* The result is t1 + t2, after normalization.*/
   ddb[0] = t1 + t2;
   e = t1 + t2;
   e -= t1;
   ddb[1] = t2 - e;
   
   return;
}

Here is the caller graph for this function: