Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] View of /branches/vis12/test/heron/heron.diderot
ViewVC logotype

View of /branches/vis12/test/heron/heron.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2361 - (download) (annotate)
Sun Apr 7 16:33:14 2013 UTC (7 years, 10 months ago) by jhr
File size: 693 byte(s)
  changes from trunk (e.g., 32-bit support)
/*! \file heron.diderot
 *
 * \author John Reppy
 *
 * This example illustrates using Heron's method for computing square roots.  It also
 * serves as a test for loading data sequences.
 */

/*
 * COPYRIGHT (c) 2012 The Diderot Project (http://diderot-language.cs.uchicago.edu)
 * All rights reserved.
 */

input int{} args = load("numbers.nrrd");
//int{} args = load("numbers.nrrd");
int nArgs = length(args);
input real eps = 0.00001;
input bool dummy;

strand SqRoot (real arg)
{
    output real root = arg;

    update {
	if (dummy) root = arg+1;
	root = (root + arg/root) / 2.0;
	if (|root^2 - arg| / arg < eps)
	    stabilize;
    }
}

initially { SqRoot(args{i}) | i in 0 .. nArgs-1 };

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0