Mercurial > octave

comparison src/DLD-FUNCTIONS/bsxfun.cc @ 14501:60e5cf354d80

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions. * __contourc__.cc, __delaunayn__.cc, __dispatch__.cc, __dsearchn__.cc, __fltk_uigetfile__.cc, __glpk__.cc, __lin_interpn__.cc, __magick_read__.cc, __pchip_deriv__.cc, __qp__.cc, __voronoi__.cc, besselj.cc, betainc.cc, bsxfun.cc, cellfun.cc, chol.cc, conv2.cc, convhulln.cc, dassl.cc, det.cc, dlmread.cc, dmperm.cc, dot.cc, eig.cc, eigs.cc, fft.cc, fft2.cc, filter.cc, find.cc, gammainc.cc, gcd.cc, givens.cc, hess.cc, hex2num.cc, inv.cc, kron.cc, lookup.cc, lsode.cc, lu.cc, luinc.cc, matrix_type.cc, max.cc, mgorth.cc, nproc.cc, qr.cc, quad.cc, quadcc.cc, qz.cc, rand.cc, rcond.cc, regexp.cc, schur.cc, spparms.cc, sqrtm.cc, str2double.cc, strfind.cc, sub2ind.cc, svd.cc, syl.cc, time.cc, tril.cc, tsearch.cc: Update %!tests in DLD-FUNCTIONS/ directory with Octave coding conventions.
author Rik <octave@nomad.inbox5.com>
date Tue, 27 Mar 2012 22:46:45 -0700
parents 72c96de7a403
children 5bc9b9cb4362
comparison
equal deleted inserted replaced
14498:36cfbd23fe9f 14501:60e5cf354d80
660 %!shared a, b, c, f 660 %!shared a, b, c, f
661 %! a = randn (4, 4); 661 %! a = randn (4, 4);
662 %! b = mean (a, 1); 662 %! b = mean (a, 1);
663 %! c = mean (a, 2); 663 %! c = mean (a, 2);
664 %! f = @minus; 664 %! f = @minus;
665 %!error(bsxfun (f)); 665 %!error (bsxfun (f))
666 %!error(bsxfun (f, a)); 666 %!error (bsxfun (f, a))
667 %!error(bsxfun (a, b)); 667 %!error (bsxfun (a, b))
668 %!error(bsxfun (a, b, c)); 668 %!error (bsxfun (a, b, c))
669 %!error(bsxfun (f, a, b, c)); 669 %!error (bsxfun (f, a, b, c))
670 %!error(bsxfun (f, ones(4, 0), ones(4, 4))) 670 %!error (bsxfun (f, ones (4, 0), ones (4, 4)))
671 %!assert(bsxfun (f, ones(4, 0), ones(4, 1)), zeros(4, 0)); 671 %!assert (bsxfun (f, ones (4, 0), ones (4, 1)), zeros (4, 0))
672 %!assert(bsxfun (f, ones(1, 4), ones(4, 1)), zeros(4, 4)); 672 %!assert (bsxfun (f, ones (1, 4), ones (4, 1)), zeros (4, 4))
673 %!assert(bsxfun (f, a, b), a - repmat(b, 4, 1)); 673 %!assert (bsxfun (f, a, b), a - repmat (b, 4, 1))
674 %!assert(bsxfun (f, a, c), a - repmat(c, 1, 4)); 674 %!assert (bsxfun (f, a, c), a - repmat (c, 1, 4))
675 %!assert(bsxfun ("minus", ones(1, 4), ones(4, 1)), zeros(4, 4)); 675 %!assert (bsxfun ("minus", ones (1, 4), ones (4, 1)), zeros (4, 4))
676 676
677 %!shared a, b, c, f 677 %!shared a, b, c, f
678 %! a = randn (4, 4); 678 %! a = randn (4, 4);
679 %! a(1) *= 1i; 679 %! a(1) *= 1i;
680 %! b = mean (a, 1); 680 %! b = mean (a, 1);
681 %! c = mean (a, 2); 681 %! c = mean (a, 2);
682 %! f = @minus; 682 %! f = @minus;
683 %!error(bsxfun (f)); 683 %!error (bsxfun (f))
684 %!error(bsxfun (f, a)); 684 %!error (bsxfun (f, a))
685 %!error(bsxfun (a, b)); 685 %!error (bsxfun (a, b))
686 %!error(bsxfun (a, b, c)); 686 %!error (bsxfun (a, b, c))
687 %!error(bsxfun (f, a, b, c)); 687 %!error (bsxfun (f, a, b, c))
688 %!error(bsxfun (f, ones(4, 0), ones(4, 4))) 688 %!error (bsxfun (f, ones (4, 0), ones (4, 4)))
689 %!assert(bsxfun (f, ones(4, 0), ones(4, 1)), zeros(4, 0)); 689 %!assert (bsxfun (f, ones (4, 0), ones (4, 1)), zeros (4, 0))
690 %!assert(bsxfun (f, ones(1, 4), ones(4, 1)), zeros(4, 4)); 690 %!assert (bsxfun (f, ones (1, 4), ones (4, 1)), zeros (4, 4))
691 %!assert(bsxfun (f, a, b), a - repmat(b, 4, 1)); 691 %!assert (bsxfun (f, a, b), a - repmat (b, 4, 1))
692 %!assert(bsxfun (f, a, c), a - repmat(c, 1, 4)); 692 %!assert (bsxfun (f, a, c), a - repmat (c, 1, 4))
693 %!assert(bsxfun ("minus", ones(1, 4), ones(4, 1)), zeros(4, 4)); 693 %!assert (bsxfun ("minus", ones (1, 4), ones (4, 1)), zeros (4, 4))
694 694
695 %!shared a, b, c, f 695 %!shared a, b, c, f
696 %! a = randn (4, 4); 696 %! a = randn (4, 4);
697 %! a(end) *= 1i; 697 %! a(end) *= 1i;
698 %! b = mean (a, 1); 698 %! b = mean (a, 1);
699 %! c = mean (a, 2); 699 %! c = mean (a, 2);
700 %! f = @minus; 700 %! f = @minus;
701 %!error(bsxfun (f)); 701 %!error (bsxfun (f))
702 %!error(bsxfun (f, a)); 702 %!error (bsxfun (f, a))
703 %!error(bsxfun (a, b)); 703 %!error (bsxfun (a, b))
704 %!error(bsxfun (a, b, c)); 704 %!error (bsxfun (a, b, c))
705 %!error(bsxfun (f, a, b, c)); 705 %!error (bsxfun (f, a, b, c))
706 %!error(bsxfun (f, ones(4, 0), ones(4, 4))) 706 %!error (bsxfun (f, ones (4, 0), ones (4, 4)))
707 %!assert(bsxfun (f, ones(4, 0), ones(4, 1)), zeros(4, 0)); 707 %!assert (bsxfun (f, ones (4, 0), ones (4, 1)), zeros (4, 0))
708 %!assert(bsxfun (f, ones(1, 4), ones(4, 1)), zeros(4, 4)); 708 %!assert (bsxfun (f, ones (1, 4), ones (4, 1)), zeros (4, 4))
709 %!assert(bsxfun (f, a, b), a - repmat(b, 4, 1)); 709 %!assert (bsxfun (f, a, b), a - repmat (b, 4, 1))
710 %!assert(bsxfun (f, a, c), a - repmat(c, 1, 4)); 710 %!assert (bsxfun (f, a, c), a - repmat (c, 1, 4))
711 %!assert(bsxfun ("minus", ones(1, 4), ones(4, 1)), zeros(4, 4)); 711 %!assert (bsxfun ("minus", ones (1, 4), ones (4, 1)), zeros (4, 4))
712 712
713 %!shared a, b, c, f 713 %!shared a, b, c, f
714 %! a = randn (4, 4); 714 %! a = randn (4, 4);
715 %! b = a (1, :); 715 %! b = a (1, :);
716 %! c = a (:, 1); 716 %! c = a (:, 1);
717 %! f = @(x, y) x == y; 717 %! f = @(x, y) x == y;
718 %!error(bsxfun (f)); 718 %!error (bsxfun (f))
719 %!error(bsxfun (f, a)); 719 %!error (bsxfun (f, a))
720 %!error(bsxfun (a, b)); 720 %!error (bsxfun (a, b))
721 %!error(bsxfun (a, b, c)); 721 %!error (bsxfun (a, b, c))
722 %!error(bsxfun (f, a, b, c)); 722 %!error (bsxfun (f, a, b, c))
723 %!error(bsxfun (f, ones(4, 0), ones(4, 4))) 723 %!error (bsxfun (f, ones (4, 0), ones (4, 4)))
724 %!assert(bsxfun (f, ones(4, 0), ones(4, 1)), zeros(4, 0, "logical")); 724 %!assert (bsxfun (f, ones (4, 0), ones (4, 1)), zeros (4, 0, "logical"))
725 %!assert(bsxfun (f, ones(1, 4), ones(4, 1)), ones(4, 4, "logical")); 725 %!assert (bsxfun (f, ones (1, 4), ones (4, 1)), ones (4, 4, "logical"))
726 %!assert(bsxfun (f, a, b), a == repmat(b, 4, 1)); 726 %!assert (bsxfun (f, a, b), a == repmat (b, 4, 1))
727 %!assert(bsxfun (f, a, c), a == repmat(c, 1, 4)); 727 %!assert (bsxfun (f, a, c), a == repmat (c, 1, 4))
728 728
729 %!shared a, b, c, d, f 729 %!shared a, b, c, d, f
730 %! a = randn (4, 4, 4); 730 %! a = randn (4, 4, 4);
731 %! b = mean (a, 1); 731 %! b = mean (a, 1);
732 %! c = mean (a, 2); 732 %! c = mean (a, 2);
733 %! d = mean (a, 3); 733 %! d = mean (a, 3);
734 %! f = @minus; 734 %! f = @minus;
735 %!error(bsxfun (f, ones([4, 0, 4]), ones([4, 4, 4]))); 735 %!error (bsxfun (f, ones ([4, 0, 4]), ones ([4, 4, 4])))
736 %!assert(bsxfun (f, ones([4, 0, 4]), ones([4, 1, 4])), zeros([4, 0, 4])); 736 %!assert (bsxfun (f, ones ([4, 0, 4]), ones ([4, 1, 4])), zeros ([4, 0, 4]))
737 %!assert(bsxfun (f, ones([4, 4, 0]), ones([4, 1, 1])), zeros([4, 4, 0])); 737 %!assert (bsxfun (f, ones ([4, 4, 0]), ones ([4, 1, 1])), zeros ([4, 4, 0]))
738 %!assert(bsxfun (f, ones([1, 4, 4]), ones([4, 1, 4])), zeros([4, 4, 4])); 738 %!assert (bsxfun (f, ones ([1, 4, 4]), ones ([4, 1, 4])), zeros ([4, 4, 4]))
739 %!assert(bsxfun (f, ones([4, 4, 1]), ones([4, 1, 4])), zeros([4, 4, 4])); 739 %!assert (bsxfun (f, ones ([4, 4, 1]), ones ([4, 1, 4])), zeros ([4, 4, 4]))
740 %!assert(bsxfun (f, ones([4, 1, 4]), ones([1, 4, 4])), zeros([4, 4, 4])); 740 %!assert (bsxfun (f, ones ([4, 1, 4]), ones ([1, 4, 4])), zeros ([4, 4, 4]))
741 %!assert(bsxfun (f, ones([4, 1, 4]), ones([1, 4, 1])), zeros([4, 4, 4])); 741 %!assert (bsxfun (f, ones ([4, 1, 4]), ones ([1, 4, 1])), zeros ([4, 4, 4]))
742 %!assert(bsxfun (f, a, b), a - repmat(b, [4, 1, 1])); 742 %!assert (bsxfun (f, a, b), a - repmat (b, [4, 1, 1]))
743 %!assert(bsxfun (f, a, c), a - repmat(c, [1, 4, 1])); 743 %!assert (bsxfun (f, a, c), a - repmat (c, [1, 4, 1]))
744 %!assert(bsxfun (f, a, d), a - repmat(d, [1, 1, 4])); 744 %!assert (bsxfun (f, a, d), a - repmat (d, [1, 1, 4]))
745 %!assert(bsxfun ("minus", ones([4, 0, 4]), ones([4, 1, 4])), zeros([4, 0, 4])); 745 %!assert (bsxfun ("minus", ones ([4, 0, 4]), ones ([4, 1, 4])), zeros ([4, 0, 4]))
746 746
747 %% The below is a very hard case to treat 747 %% The test below is a very hard case to treat
748 %!assert(bsxfun (f, ones([4, 1, 4, 1]), ones([1, 4, 1, 4])), zeros([4, 4, 4, 4])); 748 %!assert (bsxfun (f, ones ([4, 1, 4, 1]), ones ([1, 4, 1, 4])), zeros ([4, 4, 4, 4]));
749 749
750 %!shared a, b, aa, bb 750 %!shared a, b, aa, bb
751 %! a = randn (3, 1, 3); 751 %! a = randn (3, 1, 3);
752 %! aa = a(:, ones (1, 3), :, ones (1, 3)); 752 %! aa = a(:, ones (1, 3), :, ones (1, 3));
753 %! b = randn (1, 3, 3, 3); 753 %! b = randn (1, 3, 3, 3);
754 %! bb = b(ones (1, 3), :, :, :); 754 %! bb = b(ones (1, 3), :, :, :);
755 %!assert (bsxfun (@plus, a, b), aa + bb); 755 %!assert (bsxfun (@plus, a, b), aa + bb)
756 %!assert (bsxfun (@minus, a, b), aa - bb); 756 %!assert (bsxfun (@minus, a, b), aa - bb)
757 %!assert (bsxfun (@times, a, b), aa .* bb); 757 %!assert (bsxfun (@times, a, b), aa .* bb)
758 %!assert (bsxfun (@rdivide, a, b), aa ./ bb); 758 %!assert (bsxfun (@rdivide, a, b), aa ./ bb)
759 %!assert (bsxfun (@ldivide, a, b), aa .\ bb); 759 %!assert (bsxfun (@ldivide, a, b), aa .\ bb)
760 %!assert (bsxfun (@power, a, b), aa .^ bb); 760 %!assert (bsxfun (@power, a, b), aa .^ bb)
761 %!assert (bsxfun (@power, abs (a), b), abs (aa) .^ bb); 761 %!assert (bsxfun (@power, abs (a), b), abs (aa) .^ bb)
762 %!assert (bsxfun (@eq, round (a), round (b)), round (aa) == round (bb)); 762 %!assert (bsxfun (@eq, round (a), round (b)), round (aa) == round (bb))
763 %!assert (bsxfun (@ne, round (a), round (b)), round (aa) != round (bb)); 763 %!assert (bsxfun (@ne, round (a), round (b)), round (aa) != round (bb))
764 %!assert (bsxfun (@lt, a, b), aa < bb); 764 %!assert (bsxfun (@lt, a, b), aa < bb)
765 %!assert (bsxfun (@le, a, b), aa <= bb); 765 %!assert (bsxfun (@le, a, b), aa <= bb)
766 %!assert (bsxfun (@gt, a, b), aa > bb); 766 %!assert (bsxfun (@gt, a, b), aa > bb)
767 %!assert (bsxfun (@ge, a, b), aa >= bb); 767 %!assert (bsxfun (@ge, a, b), aa >= bb)
768 %!assert (bsxfun (@min, a, b), min (aa, bb)); 768 %!assert (bsxfun (@min, a, b), min (aa, bb))
769 %!assert (bsxfun (@max, a, b), max (aa, bb)); 769 %!assert (bsxfun (@max, a, b), max (aa, bb))
770 %!assert (bsxfun (@and, a > 0, b > 0), (aa > 0) & (bb > 0)); 770 %!assert (bsxfun (@and, a > 0, b > 0), (aa > 0) & (bb > 0))
771 %!assert (bsxfun (@or, a > 0, b > 0), (aa > 0) | (bb > 0)); 771 %!assert (bsxfun (@or, a > 0, b > 0), (aa > 0) | (bb > 0))
772 772
773 %% Test automatic bsxfun 773 %% Test automatic bsxfun
774 % 774 %
775 %!test 775 %!test
776 %! funs = {@plus, @minus, @times, @rdivide, @ldivide, @power, @max, @min, \ 776 %! funs = {@plus, @minus, @times, @rdivide, @ldivide, @power, @max, @min, \
777 %! @rem, @mod, @atan2, @hypot, @eq, @ne, @lt, @le, @gt, @ge, \ 777 %! @rem, @mod, @atan2, @hypot, @eq, @ne, @lt, @le, @gt, @ge, \
778 %! @and, @or, @xor }; 778 %! @and, @or, @xor };
779 %! 779 %!
780 %! float_types = {@single, @double}; 780 %! float_types = {@single, @double};
781 %! int_types = {@int8, @int16, @int32, @int64, \ 781 %! int_types = {@int8, @int16, @int32, @int64, \
782 %! @uint8, @uint16, @uint32, @uint64}; 782 %! @uint8, @uint16, @uint32, @uint64};
783 %! 783 %!
784 %! x = rand (3)*10-5; 784 %! x = rand (3) * 10-5;
785 %! y = rand (3,1)*10-5; 785 %! y = rand (3,1) * 10-5;
786 %! 786 %!
787 %! for i=1:length (funs) 787 %! for i=1:length (funs)
788 %! for j = 1:length(float_types) 788 %! for j = 1:length(float_types)
789 %! for k = 1:length(int_types) 789 %! for k = 1:length(int_types)
790 %! 790 %!