Mercurial > octave-nkf

annotate doc/interpreter/sparseimages.m @ 11033:d9c8916bb9dd

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Untabify a few remaining .m scripts.
author Rik <octave@nomad.inbox5.com>
date Tue, 28 Sep 2010 13:24:21 -0700
parents e76d8c767584
children 4cbb08163a1f
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
8920
eb63fbe60fab update copyright notices
John W. Eaton <jwe@octave.org>
parents: 7540
diff changeset
1 ## Copyright (C) 2006, 2007, 2008 David Bateman
7018
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
2 ##
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
3 ## This file is part of Octave.
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
4 ##
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
5 ## Octave is free software; you can redistribute it and/or modify it
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
6 ## under the terms of the GNU General Public License as published by
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
7 ## the Free Software Foundation; either version 3 of the License, or (at
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
8 ## your option) any later version.
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
9 ##
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
10 ## Octave is distributed in the hope that it will be useful, but
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
13 ## General Public License for more details.
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
14 ##
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
15 ## You should have received a copy of the GNU General Public License
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
16 ## along with Octave; see the file COPYING. If not, see
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
17 ## <http://www.gnu.org/licenses/>.
fd42779a8428 [project @ 2007-10-13 00:52:12 by jwe]
jwe
parents: 6257
diff changeset
18
6126
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
19 function sparseimages (nm, typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
20 set_print_size ();
7257
e8a3e3e33490 [project @ 2007-12-05 01:18:41 by jwe]
jwe
parents: 7256
diff changeset
21 if (strcmp (typ, "png"))
e8a3e3e33490 [project @ 2007-12-05 01:18:41 by jwe]
jwe
parents: 7256
diff changeset
22 set (0, "defaulttextfontname", "*");
e8a3e3e33490 [project @ 2007-12-05 01:18:41 by jwe]
jwe
parents: 7256
diff changeset
23 endif
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
24
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
25 if (! isempty (findstr (octave_config_info ("DEFS"), "HAVE_COLAMD"))
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
26 && ! isempty (findstr (octave_config_info ("DEFS"), "HAVE_CHOLMOD"))
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
27 && ! isempty (findstr (octave_config_info ("DEFS"), "HAVE_UMFPACK")))
6126
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
28 if (strcmp(typ,"txt"))
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
29 txtimages (nm, 15, typ);
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
30 else
6126
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
31 if (strcmp (nm, "gplot"))
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
32 gplotimages ("gplot", typ);
6126
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
33 elseif (strcmp (nm, "grid"))
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
34 femimages ("grid", typ);
6126
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
35 else
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
36 otherimages (nm, 200, typ);
6126
c0011bf8f815 [project @ 2006-10-31 20:08:30 by jwe]
jwe
parents: 6040
diff changeset
37 endif
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
38 endif
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
39 else ## There is no sparse matrix implementation available because
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
40 ## of missing libraries, plot sombreros instead
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
41 sombreroimage (nm, typ);
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
42 endif
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
43 endfunction
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
44
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
45 function set_print_size ()
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
46 image_size = [5.0, 3.5]; # in inches, 16:9 format
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
47 border = 0; # For postscript use 50/72
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
48 set (0, "defaultfigurepapertype", "<custom>");
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
49 set (0, "defaultfigurepaperorientation", "landscape");
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
50 set (0, "defaultfigurepapersize", image_size + 2*border);
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
51 set (0, "defaultfigurepaperposition", [border, border, image_size]);
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
52 endfunction
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
53
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
54 function hide_output ()
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6180
diff changeset
55 f = figure (1);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6180
diff changeset
56 set (f, "visible", "off");
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
57 endfunction
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
58
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
59 function gplotimages (nm, typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
60 hide_output ();
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
61 if (strcmp (typ, "eps"))
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
62 d_typ = "-depsc2";
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
63 else
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
64 d_typ = cstrcat ("-d", typ);
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
65 endif
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
66
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
67 A = sparse ([2,6,1,3,2,4,3,5,4,6,1,5],
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
68 [1,1,2,2,3,3,4,4,5,5,6,6], 1, 6, 6);
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
69 xy = [0,4,8,6,4,2;5,0,5,7,5,7]';
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
70 gplot (A, xy)
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
71 print (cstrcat (nm, ".", typ), d_typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
72 hide_output ();
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
73 endfunction
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
74
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
75 function txtimages(nm,n,typ)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
76 a = 10*speye(n) + sparse(1:n,ceil([1:n]/2),1,n,n) + ...
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
77 sparse(ceil([1:n]/2),1:n,1,n,n);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
78 if (strcmp (nm, "gplot") || strcmp (nm, "grid"))
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
79 fid = fopen (sprintf ("%s.txt", nm), "wt");
7256
be3280acbdb8 [project @ 2007-12-05 01:08:17 by jwe]
jwe
parents: 7018
diff changeset
80 fputs (fid, "\n");
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
81 fputs (fid, "+---------------------------------+\n");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
82 fputs (fid, "| Image unavailable in text mode. |\n");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
83 fputs (fid, "+---------------------------------+\n");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
84 fclose (fid);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
85 elseif (strcmp (nm, "spmatrix"))
7540
3422f39573b1 strcat.m: Matlab compatibility, with cstrcat.m replacing conventional strcat.m.
Ben Abbott <bpabbott@mac.com>
parents: 7262
diff changeset
86 printsparse(a,cstrcat("spmatrix.",typ));
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
87 else
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
88 if (!isempty(findstr(octave_config_info ("DEFS"),"HAVE_COLAMD")) &&
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
89 !isempty(findstr(octave_config_info ("DEFS"),"HAVE_CHOLMOD")))
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
90 if (strcmp (nm, "spchol"))
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
91 r1 = chol(a);
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
92 printsparse(r1,cstrcat("spchol.",typ));
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
93 elseif (strcmp (nm, "spcholperm"))
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
94 [r2,p2,q2]=chol(a);
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
95 printsparse(r2,cstrcat("spcholperm.",typ));
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
96 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
97 ## printf("Text NNZ: Matrix %d, Chol %d, PermChol %d\n",nnz(a),nnz(r1),nnz(r2));
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
98 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
99 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
100 endfunction
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
101
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
102 function otherimages(nm,n,typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
103 hide_output ();
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
104 if (strcmp (typ, "eps"))
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
105 d_typ = "-depsc2";
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
106 else
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
107 d_typ = cstrcat ("-d", typ);
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
108 endif
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
109
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
110 a = 10*speye(n) + sparse(1:n,ceil([1:n]/2),1,n,n) + ...
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
111 sparse(ceil([1:n]/2),1:n,1,n,n);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
112 if (strcmp (nm, "spmatrix"))
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
113 spy(a);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
114 axis("ij")
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
115 print(cstrcat("spmatrix.",typ), d_typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
116 hide_output ();
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
117 else
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
118 if (!isempty(findstr(octave_config_info ("DEFS"),"HAVE_COLAMD")) &&
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
119 !isempty(findstr(octave_config_info ("DEFS"),"HAVE_CHOLMOD")))
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
120 if (strcmp (nm, "spchol"))
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
121 r1 = chol(a);
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
122 spy(r1);
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
123 axis("ij")
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
124 print(cstrcat("spchol.",typ), d_typ)
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
125 hide_output ();
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
126 elseif (strcmp (nm, "spcholperm"))
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
127 [r2,p2,q2]=chol(a);
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
128 spy(r2);
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
129 axis("ij")
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
130 print(cstrcat("spcholperm.",typ), d_typ)
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
131 hide_output ();
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
132 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
133 ## printf("Image NNZ: Matrix %d, Chol %d, PermChol %d\n",nnz(a),nnz(r1),nnz(r2));
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
134 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
135 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
136 endfunction
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
137
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
138 function printsparse(a,nm)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
139 fid = fopen (nm,"wt");
7256
be3280acbdb8 [project @ 2007-12-05 01:08:17 by jwe]
jwe
parents: 7018
diff changeset
140 fputs (fid, "\n");
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
141 for i = 1:size(a,1)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
142 if (rem(i,5) == 0)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
143 fprintf (fid," %2d - ", i);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
144 else
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
145 fprintf (fid," | ");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
146 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
147 for j = 1:size(a,2)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
148 if (a(i,j) == 0)
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
149 fprintf(fid," ")
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
150 else
11033
d9c8916bb9dd Untabify a few remaining .m scripts.
Rik <octave@nomad.inbox5.com>
parents: 9914
diff changeset
151 fprintf(fid," *")
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
152 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
153 endfor
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
154 fprintf(fid,"\n")
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
155 endfor
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
156 fprintf(fid," |-");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
157 for j=1:size(a,2)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
158 if (rem(j,5)==0)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
159 fprintf(fid,"-|");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
160 else
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
161 fprintf(fid,"--");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
162 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
163 endfor
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
164 fprintf(fid,"\n")
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
165 fprintf(fid," ");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
166 for j=1:size(a,2)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
167 if (rem(j,5)==0)
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
168 fprintf(fid,"%2d",j);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
169 else
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
170 fprintf(fid," ");
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
171 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
172 endfor
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
173 fclose(fid);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
174 endfunction
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
175
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
176 function femimages (nm,typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
177 hide_output ();
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
178 if (strcmp (typ, "eps"))
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
179 d_typ = "-depsc2";
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
180 else
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
181 d_typ = cstrcat ("-d", typ);
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
182 endif
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
183
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
184 if (!isempty(findstr(octave_config_info ("DEFS"),"HAVE_COLAMD")) &&
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
185 !isempty(findstr(octave_config_info ("DEFS"),"HAVE_CHOLMOD")) &&
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
186 !isempty(findstr(octave_config_info ("DEFS"),"HAVE_UMFPACK")))
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
187 ## build a rectangle
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
188 node_y = [1;1.2;1.5;1.8;2]*ones(1,11);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
189 node_x = ones(5,1)*[1,1.05,1.1,1.2,1.3,1.5,1.7,1.8,1.9,1.95,2];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
190 nodes = [node_x(:), node_y(:)];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
191
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
192 [h,w] = size(node_x);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
193 elems = [];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
194 for idx = 1:w-1
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
195 widx = (idx-1)*h;
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
196 elems = [elems; widx+[(1:h-1);(2:h);h+(1:h-1)]'];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
197 elems = [elems; widx+[(2:h);h+(2:h);h+(1:h-1)]'];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
198 endfor
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
199
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
200 E = size(elems,1); #No. of elements
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
201 N = size(nodes,1); #No. of elements
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
202 D = size(elems,2); #dimentions+1
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
203
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
204 ## Plot FEM Geometry
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
205 elemx = elems(:,[1,2,3,1])';
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
206 xelems = reshape( nodes(elemx, 1), 4, E);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
207 yelems = reshape( nodes(elemx, 2), 4, E);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
208
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
209 ## Set element conductivity
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
210 conductivity = [1*ones(1,16),2*ones(1,48),1*ones(1,16)];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
211
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
212 ## Dirichlet boundary conditions
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
213 D_nodes = [1:5, 51:55];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
214 D_value = [10*ones(1,5), 20*ones(1,5)];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
215
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
216 ## Neumann boundary conditions
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
217 ## Note that N_value must be normalized by the boundary
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
218 ## length and element conductivity
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
219 N_nodes = [];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
220 N_value = [];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
221
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
222 ## Calculate connectivity matrix
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
223 C = sparse((1:D*E), reshape(elems',D*E,1),1, D*E, N);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
224
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
225 ## Calculate stiffness matrix
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
226 Siidx = floor([0:D*E-1]'/D)*D*ones(1,D) + ones(D*E,1)*(1:D) ;
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
227 Sjidx = [1:D*E]'*ones(1,D);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
228 Sdata = zeros(D*E,D);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
229 dfact = prod(2:(D-1));
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
230 for j = 1:E
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
231 a = inv([ ones(D,1), nodes( elems(j,:), : ) ]);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
232 const = conductivity(j)*2/dfact/abs(det(a));
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
233 Sdata(D*(j-1)+(1:D),:)= const * a(2:D,:)'*a(2:D,:);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
234 endfor
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
235
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
236 ## Element-wise system matrix
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
237 SE = sparse(Siidx,Sjidx,Sdata);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
238 ## Global system matrix
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
239 S = C'* SE *C;
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
240
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
241 ## Set Dirichlet boundary
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
242 V = zeros(N,1);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
243 V(D_nodes) = D_value;
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
244 idx = 1:N;
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
245 idx(D_nodes) = [];
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
246
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
247 ## Set Neumann boundary
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
248 Q = zeros(N,1);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
249 Q(N_nodes) = N_value; # FIXME
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
250
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
251 V(idx) = S(idx,idx)\( Q(idx) - S(idx,D_nodes)*V(D_nodes) );
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
252
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
253 velems = reshape( V(elemx), 4, E);
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
254
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
255 sz = size(xelems,2);
6178
830235f4984f [project @ 2006-11-17 00:16:57 by jwe]
jwe
parents: 6126
diff changeset
256
830235f4984f [project @ 2006-11-17 00:16:57 by jwe]
jwe
parents: 6126
diff changeset
257 plot3 (xelems, yelems, velems);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6180
diff changeset
258 view (10, 10);
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
259 print(cstrcat(nm,".",typ), d_typ)
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
260 hide_output ();
6003
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
261 endif
7f8e4f5e2eb7 [project @ 2006-09-26 21:09:31 by jwe]
jwe
parents:
diff changeset
262 endfunction
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
263
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
264 ## There is no sparse matrix implementation available because of missing
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
265 ## libraries, plot sombreros instead. Also plot a nice title that we are
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
266 ## sorry about that.
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
267 function sombreroimage (nm, typ)
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
268 if (strcmp (typ, "txt"))
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
269 fid = fopen (sprintf ("%s.txt", nm), "wt");
7256
be3280acbdb8 [project @ 2007-12-05 01:08:17 by jwe]
jwe
parents: 7018
diff changeset
270 fputs (fid, "\n");
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
271 fputs (fid, "+---------------------------------------+\n");
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
272 fputs (fid, "| Image unavailable because of a |\n");
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
273 fputs (fid, "| missing sparse matrix implementation. |\n");
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
274 fputs (fid, "+---------------------------------------+\n");
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
275 fclose (fid);
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
276 return;
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
277 else ## if (!strcmp (typ, "txt"))
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
278
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
279 hide_output ();
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
280 if (strcmp (typ, "eps"))
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
281 d_typ = "-depsc2";
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
282 else
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
283 d_typ = cstrcat ("-d", typ);
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
284 endif
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
285
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
286 x = y = linspace (-8, 8, 41)';
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
287 [xx, yy] = meshgrid (x, y);
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
288 r = sqrt (xx .^ 2 + yy .^ 2) + eps;
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
289 z = sin (r) ./ r;
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
290 unwind_protect
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6180
diff changeset
291 mesh (x, y, z);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6180
diff changeset
292 title ("Sorry, graphics not available because octave was\\ncompiled without the sparse matrix implementation.");
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
293 unwind_protect_cleanup
9914
e76d8c767584 Switch to color output for PostScript documentation images
Rik <rdrider0-list@yahoo.com>
parents: 9912
diff changeset
294 print (cstrcat (nm, ".", typ), d_typ);
9912
e9fe12c1b0c0 Set default papersize to the desired imagesize
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
295 hide_output ();
6040
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
296 end_unwind_protect
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
297 endif
2bcb2fa24bbf [project @ 2006-10-09 15:13:29 by jwe]
jwe
parents: 6003
diff changeset
298 endfunction
7262
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
299
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
300 ## generate something for the texinfo @image command to process
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
301 function image_as_txt(nm)
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
302 fid = fopen (sprintf ("%s.txt", nm), "wt");
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
303 fputs (fid, "\n");
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
304 fputs (fid, "+---------------------------------+\n");
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
305 fputs (fid, "| Image unavailable in text mode. |\n");
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
306 fputs (fid, "+---------------------------------+\n");
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
307 fclose (fid);
16e61aecddc1 [project @ 2007-12-05 20:50:47 by jwe]
jwe
parents: 7257
diff changeset
308 endfunction