Mercurial > octave
annotate scripts/linear-algebra/cross.m @ 21308:c53bfd6d8e08
maint: Use American spelling for "behavior".
* Array.cc, randmtzig.c, file-ops.cc, oct-inttypes.h, cross.m, bicgstab.m,
kurtosis.m, skewness.m, build-sparse-tests.sh, __ilu__.cc, __magick_read__.cc:
Use American spelling for "behavior".
author | Rik <rik@octave.org> |
---|---|
date | Fri, 19 Feb 2016 09:40:59 -0800 |
parents | 516bb87ea72e |
children | bac0d6f07a3e |
rev | line source |
---|---|
19697
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
17744
diff
changeset
|
1 ## Copyright (C) 1995-2015 Kurt Hornik |
3426 | 2 ## |
3922 | 3 ## This file is part of Octave. |
4 ## | |
5 ## Octave is free software; you can redistribute it and/or modify it | |
6 ## under the terms of the GNU General Public License as published by | |
7016 | 7 ## the Free Software Foundation; either version 3 of the License, or (at |
8 ## your option) any later version. | |
3426 | 9 ## |
3922 | 10 ## Octave is distributed in the hope that it will be useful, but |
2540 | 11 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
3426 | 13 ## General Public License for more details. |
14 ## | |
2540 | 15 ## You should have received a copy of the GNU General Public License |
7016 | 16 ## along with Octave; see the file COPYING. If not, see |
17 ## <http://www.gnu.org/licenses/>. | |
2540 | 18 |
3321 | 19 ## -*- texinfo -*- |
20852
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20160
diff
changeset
|
20 ## @deftypefn {} {} cross (@var{x}, @var{y}) |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20160
diff
changeset
|
21 ## @deftypefnx {} {} cross (@var{x}, @var{y}, @var{dim}) |
20160
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
22 ## Compute the vector cross product of two 3-dimensional vectors @var{x} and |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
23 ## @var{y}. |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
24 ## |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
25 ## If @var{x} and @var{y} are matrices, the cross product is applied along the |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
26 ## first dimension with three elements. |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
27 ## |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
28 ## The optional argument @var{dim} forces the cross product to be calculated |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
29 ## along the specified dimension. |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
30 ## |
03b9d17a2d95
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
31 ## Example Code: |
3418 | 32 ## |
3321 | 33 ## @example |
34 ## @group | |
35 ## cross ([1,1,0], [0,1,1]) | |
36 ## @result{} [ 1; -1; 1 ] | |
37 ## @end group | |
38 ## @end example | |
3883 | 39 ## |
12546
39ca02387a32
Improve docstrings for a number of functions.
Rik <octave@nomad.inbox5.com>
parents:
12520
diff
changeset
|
40 ## @seealso{dot, curl, divergence} |
3321 | 41 ## @end deftypefn |
2540 | 42 |
5428 | 43 ## Author: Kurt Hornik <Kurt.Hornik@wu-wien.ac.at> |
2540 | 44 ## Created: 15 October 1994 |
45 ## Adapted-By: jwe | |
46 | |
4890 | 47 function z = cross (x, y, dim) |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
48 |
4890 | 49 if (nargin != 2 && nargin != 3) |
6046 | 50 print_usage (); |
4890 | 51 endif |
52 | |
53 if (ndims (x) < 3 && ndims (y) < 3 && nargin < 3) | |
21308
c53bfd6d8e08
maint: Use American spelling for "behavior".
Rik <rik@octave.org>
parents:
20852
diff
changeset
|
54 ## COMPATIBILITY -- opposite behavior for cross(row,col) |
c53bfd6d8e08
maint: Use American spelling for "behavior".
Rik <rik@octave.org>
parents:
20852
diff
changeset
|
55 ## Swap x and y in the assignments below to get the matlab behavior. |
4890 | 56 ## Better yet, fix the calling code so that it uses conformant vectors. |
57 if (columns (x) == 1 && rows (y) == 1) | |
58 warning ("cross: taking cross product of column by row"); | |
59 y = y.'; | |
60 elseif (rows (x) == 1 && columns (y) == 1) | |
61 warning ("cross: taking cross product of row by column"); | |
62 x = x.'; | |
63 endif | |
2540 | 64 endif |
65 | |
4890 | 66 if (nargin == 2) |
6024 | 67 dim = find (size (x) == 3, 1); |
11587
c792872f8942
all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
68 if (isempty (dim)) |
4890 | 69 error ("cross: must have at least one dimension with 3 elements"); |
70 endif | |
71 else | |
9877 | 72 if (size (x, dim) != 3) |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
73 error ("cross: dimension DIM must have 3 elements"); |
4890 | 74 endif |
3883 | 75 endif |
3085 | 76 |
4890 | 77 nd = ndims (x); |
78 sz = size (x); | |
9877 | 79 idx2 = idx3 = idx1 = {':'}(ones (1, nd)); |
4890 | 80 idx1(dim) = 1; |
81 idx2(dim) = 2; | |
82 idx3(dim) = 3; | |
83 | |
6157 | 84 if (size_equal (x, y)) |
9877 | 85 x1 = x(idx1{:}); |
86 x2 = x(idx2{:}); | |
87 x3 = x(idx3{:}); | |
88 y1 = y(idx1{:}); | |
89 y2 = y(idx2{:}); | |
90 y3 = y(idx3{:}); | |
91 z = cat (dim, (x2.*y3 - x3.*y2), (x3.*y1 - x1.*y3), (x1.*y2 - x2.*y1)); | |
3085 | 92 else |
11472
1740012184f9
Use uppercase for variable names in error() strings to match Info documentation. Only m-files done.
Rik <octave@nomad.inbox5.com>
parents:
10549
diff
changeset
|
93 error ("cross: X and Y must have the same dimensions"); |
2540 | 94 endif |
95 | |
96 endfunction | |
13048
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
97 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
98 |
13048
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
99 %!test |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
100 %! x = [1 0 0]; |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
101 %! y = [0 1 0]; |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
102 %! r = [0 0 1]; |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
103 %! assert (cross (x, y), r, 2e-8); |
13048
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
104 |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
105 %!test |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
106 %! x = [1 2 3]; |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
107 %! y = [4 5 6]; |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
108 %! r = [(2*6-3*5) (3*4-1*6) (1*5-2*4)]; |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
109 %! assert (cross (x, y), r, 2e-8); |
13048
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
110 |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
111 %!test |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
112 %! x = [1 0 0; 0 1 0; 0 0 1]; |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
113 %! y = [0 1 0; 0 0 1; 1 0 0]; |
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
114 %! r = [0 0 1; 1 0 0; 0 1 0]; |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
115 %! assert (cross (x, y, 2), r, 2e-8); |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
116 %! assert (cross (x, y, 1), -r, 2e-8); |
13048
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
117 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
118 %!error cross (0,0) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
119 %!error cross () |
13048
c5c94b63931f
codesprint: linear algebra tests: cross, housh, planerot, qzhess, rref
Roman Belov <romblv@gmail.com>
parents:
12546
diff
changeset
|
120 |