mandelbrotset1 revision da2e3ebdc1edfbc5028edf1354e7dd2fa69a7968
2N/A#!/bin/ksh93
2N/A
2N/A#
2N/A# CDDL HEADER START
2N/A#
2N/A# The contents of this file are subject to the terms of the
2N/A# Common Development and Distribution License (the "License").
2N/A# You may not use this file except in compliance with the License.
2N/A#
2N/A# You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
2N/A# or http://www.opensolaris.org/os/licensing.
2N/A# See the License for the specific language governing permissions
2N/A# and limitations under the License.
2N/A#
2N/A# When distributing Covered Code, include this CDDL HEADER in each
2N/A# file and include the License file at usr/src/OPENSOLARIS.LICENSE.
2N/A# If applicable, add the following below this CDDL HEADER, with the
2N/A# fields enclosed by brackets "[]" replaced with your own identifying
2N/A# information: Portions Copyright [yyyy] [name of copyright owner]
2N/A#
2N/A# CDDL HEADER END
2N/A#
2N/A
2N/A#
2N/A# Copyright 2007 Sun Microsystems, Inc. All rights reserved.
2N/A# Use is subject to license terms.
2N/A#
2N/A# ident "%Z%%M% %I% %E% SMI"
2N/A#
2N/A
2N/A#
2N/A# mandelbrotset1 - a simple mandelbrot set generation and
2N/A# parallel execution demo
2N/A#
2N/A
2N/A# Solaris needs /usr/xpg4/bin/ because the tools in /usr/bin are not POSIX-conformant
2N/Aexport PATH=/usr/xpg4/bin:/bin:/usr/bin:/usr/X11/bin:/usr/X11R6/bin:/usr/openwin/bin
2N/A
2N/Afunction printmsg
{
print -u 2 "$@"
}
function fatal_error
{
print -u 2 "${progname}: $@"
exit 1
}
function print_color
{
print -n "${symbollist:${1}:1}"
}
function mandelbrot
{
float x=$1
float y=$2
float xx
float yy
float x1=$3
float y1=$4
integer iteration=$5
integer max_iteration=$6
float mag
for (( mag=0 ; mag < max_mag && iteration < max_iteration ; iteration++ )) ; do
(( xx=x*x ))
(( yy=y*y ))
(( mag=xx+yy ))
(( y=x*y*2+y1 ))
(( x=xx-yy+x1 ))
done
print ${iteration}
return 0
}
function loop_serial
{
for (( y=y_min ; y < y_max ; y+=stepwidth )) ; do
for (( x=x_min ; x < x_max ; x+=stepwidth )) ; do
print_color $(mandelbrot ${x} ${y} ${x} ${y} 1 ${symbollistlen})
done
print
done
}
function loop_parallel
{
integer numjobs=0
# the following calculation suffers from rounding errors
integer lines_per_job=$(( ((m_height+(numcpus-1)) / numcpus) ))
printmsg "# lines_per_job=${lines_per_job}"
printmsg "# numcpus=${numcpus}"
# "renice" worker jobs
set -o bgnice
if [ "${TMPDIR}" = "" ] ; then
TMPDIR="/tmp"
fi
# try to generate a job identifer prefix which is unique across multiple hosts
jobident="job_host_$(uname -n)pid_$$_ppid${PPID}"
printmsg $"## prepare..."
for (( y=y_min ; y < y_max ; y+=(stepwidth*lines_per_job) )) ; do
rm -f "${TMPDIR}/mandelbrot_${jobident}_child_$y.joboutput"
let numjobs++
done
printmsg $"## running ${numjobs} children..."
for (( y=y_min ; y < y_max ; y+=(stepwidth*lines_per_job) )) ; do
(
for (( ; y < y_max && lines_per_job-- > 0 ; y+=stepwidth )) ; do
for (( x=x_min ; x < x_max ; x+=stepwidth )) ; do
print_color $(mandelbrot ${x} ${y} ${x} ${y} 1 ${symbollistlen})
done
print
done >"${TMPDIR}/mandelbrot_${jobident}_child_$y.joboutput"
) &
done
printmsg $"## waiting for ${numjobs} children..."
wait
printmsg $"## output:"
for (( y=y_min ; y < y_max ; y+=(stepwidth*lines_per_job) )) ; do
print "$(cat "${TMPDIR}/mandelbrot_${jobident}_child_$y.joboutput")"
rm "${TMPDIR}/mandelbrot_${jobident}_child_$y.joboutput"
done
}
function usage
{
OPTIND=0
getopts -a "${progname}" "${USAGE}" OPT '-?'
exit 2
}
# main
builtin printf
builtin cat
builtin rm
builtin sleep
builtin uname # loop_parallel needs the ksh93 builtin version to generate unique job file names
float x_max
float x_min
float y_max
float y_min
float m_width
float m_height
float max_mag
float stepwidth
integer numcpus
# make sure ${COLUMN} and ${LINES} are set
eval $(resize -u)
symbollist=' .:0123456789abcdefghijklmnopqrstuvwxyzABCDEFGHIJKLMNOPQRSTUVWXYZ%#'
symbollistlen=$(( ${#symbollist} - 1))
mode="parallel"
progname="${0}"
max_mag=400
stepwidth=0.1
numcpus=16
let m_width=COLUMNS-1 m_height=LINES-2
progname="${0}"
USAGE=$'
[-?
@(#)\$Id: mandelbrotset1 (Roland Mainz) 2007-06-05 \$
]
[+NAME?mandelbrotset1 - generate mandelbrot set fractals with ksh93]
[+DESCRIPTION?\bmandelbrotset1\b mandelbrot set fractal generator
which runs either in serial or parallel mode (using multiple worker jobs).]
[w:width?Width of fractal.]:[width]
[h:height?Height of fractal.]:[height]
[s:symbols?Symbols to build the fractal from.]:[symbolstring]
[m:mag?Magnification level.]:[magnificationlevel]
[p:stepwidth?Width per step.]:[widthperstep]
[S:serial?Run in serial mode.]
[P:parallel?Run in parallel mode.]
[M:mode?Execution mode.]:[mode]
[C:numcpus?Number of processors used for parallel execution.]:[numcpus]
[+SEE ALSO?\bjuliaset1\b(1), \bksh93\b(1)]
'
while getopts -a "${progname}" "${USAGE}" OPT ; do
# printmsg "## OPT=|${OPT}|, OPTARG=|${OPTARG}|"
case ${OPT} in
w) m_width="${OPTARG}" ;;
h) m_height="${OPTARG}" ;;
s) symbollist="${OPTARG}" ;;
m) max_mag="${OPTARG}" ;;
p) stepwidth="${OPTARG}" ;;
S) mode="serial" ;;
P) mode="parallel" ;;
M) mode="${OPTARG}" ;;
C) numcpus="${OPTARG}" ;;
*) usage ;;
esac
done
shift ${OPTIND}-1
printmsg "# width=${m_width}"
printmsg "# height=${m_height}"
printmsg "# max_mag=${max_mag}"
printmsg "# stepwidth=${stepwidth}"
printmsg "# symbollist='${symbollist}'"
printmsg "# mode=${mode}"
symbollistlen=$(( ${#symbollist} - 1))
let x_max=m_width*stepwidth/2. x_min=-x_max
let y_max=m_height*stepwidth/2. y_min=-y_max
case "${mode}" in
parallel) loop_parallel ;;
serial) loop_serial ;;
*) fatal_error $"Unknown mode \"${mode}\"."
esac
# EOF.