|
Warning: this is a htmlized version!
The original is across this link. |
#######
#
# E-scripts on programs for math.
#
# Note 1: use the eev command (defined in eev.el) and the
# ee alias (in my .zshrc) to execute parts of this file.
# Executing this file as a whole makes no sense.
#
# Note 2: be VERY careful and make sure you understand what
# you're doing.
#
# Note 3: If you use a shell other than zsh things like |&
# and the for loops may not work.
#
# Note 4: I always run as root.
#
# Note 5: some parts are too old and don't work anymore. Some
# never worked.
#
# Note 6: the definitions for the find-xxxfile commands are on my
# .emacs.
#
# Note 7: if you see a strange command check my .zshrc -- it may
# be defined there as a function or an alias.
#
# Note 8: the sections without dates are always older than the
# sections with dates.
#
# This file is at <http://angg.twu.net/a/e/math.e>
# or at <http://angg.twu.net/e/math.e.html>.
# See also <http://angg.twu.net/emacs.html>,
# <http://angg.twu.net/a/.emacs[.html]>,
# <http://angg.twu.net/a/.zshrc[.html]>,
# <http://angg.twu.net/escripts.html>,
# and <http://angg.twu.net/>.
#
#######
######
#
# Programação linear
#
######
apti libwn6 libwn-dev octave-sp
apti pcx
# (find-vldifile "libwn6.list")
# (find-fline "/usr/doc/libwn6/")
# (find-vldifile "libwn-dev.list")
# (find-fline "/usr/doc/libwn-dev/")
# (eeman "wnsplx")
# (find-vldifile "octave-sp.list")
# (find-fline "/usr/doc/octave-sp/")
cd /usr/doc/octave-sp/
for i in doc.ps semidef_prog.ps; do gzip -cd < $i.gz > /tmp/$i; done
cd /tmp/
gv doc.ps
gv semidef_prog.ps
#####
#
# Octave
#
#####
# (find-node "(octave)Index Expressions")
# (find-node "(octave)Simple Examples")
# (find-node "(octave)Script Files")
octave <<'---'
a = [1,2;3,4]
a*a
---
cat > /tmp/simple.m <<'---'
1;
function xdot = f (x, t)
r = 0.25;
k = 1.4;
a = 1.5;
b = 0.16;
c = 0.9;
d = 0.8;
xdot(1) = r*x(1)*(1 - x(1)/k) - a*x(1)*x(2)/(1 + b*x(1));
xdot(2) = c*a*x(1)*x(2)/(1 + b*x(1)) - d*x(2);
endfunction
x0 = [1; 2];
t = linspace (0, 50, 200)';
x = lsode ("f", x0, t);
# use_plplot
plot (t, x)
keyboard
---
octave /tmp/simple.m
cd /tmp/
octave
# source "simple.m"
#####
#
# geomview (slink)
# 99nov08
#
#####
# «geomview»
apti geomview
# (find-vldifile "geomview.list")
# (find-node "(geomview)Top")
# (find-fline "/usr/doc/geomview/")
# (find-fline "/usr/doc/geomview/oogltour" "Polylist file format")
# (find-fline "/usr/doc/geomview/examples/")
# (find-fline "/usr/doc/geomview/examples/example4.tcl")
lynx /usr/doc/geomview/html/geomview_toc.html
edrxnetscape /usr/doc/geomview/html/geomview_toc.html &
lynx /usr/doc/geomview/html/geomview_8.html#SEC58
# (find-fline "/usr/doc/geomview/geomview.tex" "C, perl, tcl/tk")
# (find-node "(geomview)Example4")
# (find-fline "/usr/lib/geomview/modules/")
# (find-fline "/usr/lib/geomview/modules/tcl/")
#####
#
# qhull
# 99nov08
#
#####
# «qhull»
rm -Rv /usr/src/qhull/
cd /usr/src/
tar -xvzf $S/ftp/ftp.geom.umn.edu/pub/software/qhull.tar.Z
cd /usr/src/qhull/
make |& tee om
# (find-fline "/usr/src/qhull/")
# (find-fline "/usr/src/qhull/README.txt" "if you have Geomview")
# (find-fline "/usr/src/qhull/qhull.txt")
# (find-fline "/usr/src/qhull/q_eg")
# - try 'rbox 100 | qhull G >a' and load 'a' into Geomview
# - run 'q_eg' for Geomview examples of Qhull output (see qh-eg.htm)
rm -Rv /tmp/qhull/
mkdir /tmp/qhull/
cd /usr/src/qhull/
rbox 100 | tee /tmp/qhull/a3d.rbox | qhull G > /tmp/qhull/a3d
rbox D4 12 | tee /tmp/qhull/a4d.rbox | qhull G > /tmp/qhull/a4d
cd /usr/src/qhull/
q_eg |& tee /tmp/qhull/oqeg
mv -v eg.* /tmp/qhull
# (find-fline "/tmp/qhull/")
cd /tmp/qhull/
geomview
geomview /tmp/qhull/a3d
geomview /tmp/qhull/a4d
# Problema do Leonardo Lustosa (politopo estranho):
# (find-fline "/tmp/qhull/")
# (find-fline "/tmp/qhull/a3d.rbox")
cd /usr/src/qhull/
expect -c '
proc ang {frac} {expr 2*3.14159265*$frac}
set n 200
set oscs 3
puts "3 rbox $n"
puts "$n"
for {set i 0} {$i<$n} {incr i} {
puts "[expr cos([ang $i/$n])] [expr sin([ang $i/$n])] [expr sin([ang $i/$n*$oscs])]"
}
' | tee /tmp/qhull/garagem.rbox \
| qhull G > /tmp/qhull/garagem
geomview /tmp/qhull/garagem
#####
#
# cdd+
# 99nov08
#
#####
apti gmp2 gmp2-dev
apti libg++2.8.2 libg++2.8.2-dev
apti g++
rm -Rv /usr/src/cdd+-076a1/
cd /usr/src/
tar -xvzf $S/ftp/ftp.ifor.math.ethz.ch/pub/fukuda/cdd/cdd+-076a1.tar.gz
cd /usr/src/cdd+-076a1/
# (code-c-d "cdd+" "/usr/src/cdd+-076a1/")
# (find-cdd+file "README.cdd+")
# (find-cdd+file "makefile")
# The default is to use -O3, but I just want to compile fast and run the demos
make CC=g++ \
GMPUSED=TRUE GMPLIB=/usr/lib RATLIB=gmp2 GMPINCLUDEDIR=/usr/include/gmp2 \
OPTFLAGS= all \
|& tee om
latex cddman.tex
latex cddman.tex
latex cddman.tex
rm -Rv /tmp/cdd/
mkdir /tmp/cdd/
cd /usr/src/cdd+-076a1/
cp -v ext/ccc4.ext /tmp/cdd
cp -v ine/ucube.ine /tmp/cdd/
cddr+_gmp /tmp/cdd/ucube.ine |& tee /tmp/cdd/ocu
cddr+_gmp /tmp/cdd/ccc4.ext |& tee /tmp/cdd/occ
# (find-fline "/tmp/cdd/")
# (find-cdd+file "")
# (find-cdd+file "cddman.tex" "qhull")
# (find-cdd+file "ext/")
# (find-cdd+file "ine/")
# (find-cdd+file "om")
# Old notes:
# (find-fline "$S/ftp/ftp.ifor.math.ethz.ch/pub/fukuda/cdd/cddlib-085.tar.gz")
Pgrep m/gmp/ > ~/o
# (find-fline "~/o")
cd ~/SLINK/
echo *g++*
echo *gmp*
getlinks < $S/http/www.ifor.math.ethz.ch/ifor/staff/fukuda/cdd_home/cdd.html
######
#
# Roots of a complex polynomial
#
######
# (code-c-d "oct" "/usr/src/octave-2.0.13.95/" "octave")
# (find-octnode "Polynomial Manipulations" "roots (V)")
# (find-octnode "Functions and Scripts")
# (find-vldifile "octave.list")
# (find-fline "/usr/doc/octave/")
pdsc /big/slinks2/dists/slink/main/source/math/octave_2.0.13.95-1.dsc
# (find-octfile "scripts/polynomial/roots.m")
# (find-octfile "doc/faq/")