Subtract matrices of different sizes












1












$begingroup$


I have $2$ matrices of different sizes: $8 times 3$ matrix $A$ and $2 times 3$ matrix $B$. I need to subtract from each row of $B$, a particular row of $A$ and sum these values. This I need to do for all the rows for $A$ and store the result in another matrix $C$. Is there a way to do this avoiding both for loops?



I can avoid one using bsxfun but seems like I need to use one for loop for iterating over the elements of $A$.










share|cite|improve this question











$endgroup$

















    1












    $begingroup$


    I have $2$ matrices of different sizes: $8 times 3$ matrix $A$ and $2 times 3$ matrix $B$. I need to subtract from each row of $B$, a particular row of $A$ and sum these values. This I need to do for all the rows for $A$ and store the result in another matrix $C$. Is there a way to do this avoiding both for loops?



    I can avoid one using bsxfun but seems like I need to use one for loop for iterating over the elements of $A$.










    share|cite|improve this question











    $endgroup$















      1












      1








      1





      $begingroup$


      I have $2$ matrices of different sizes: $8 times 3$ matrix $A$ and $2 times 3$ matrix $B$. I need to subtract from each row of $B$, a particular row of $A$ and sum these values. This I need to do for all the rows for $A$ and store the result in another matrix $C$. Is there a way to do this avoiding both for loops?



      I can avoid one using bsxfun but seems like I need to use one for loop for iterating over the elements of $A$.










      share|cite|improve this question











      $endgroup$




      I have $2$ matrices of different sizes: $8 times 3$ matrix $A$ and $2 times 3$ matrix $B$. I need to subtract from each row of $B$, a particular row of $A$ and sum these values. This I need to do for all the rows for $A$ and store the result in another matrix $C$. Is there a way to do this avoiding both for loops?



      I can avoid one using bsxfun but seems like I need to use one for loop for iterating over the elements of $A$.







      matrices matlab






      share|cite|improve this question















      share|cite|improve this question













      share|cite|improve this question




      share|cite|improve this question








      edited Apr 1 '18 at 10:05









      Rodrigo de Azevedo

      13k41958




      13k41958










      asked Mar 11 '14 at 19:50









      Shrinu KushagraShrinu Kushagra

      61




      61






















          2 Answers
          2






          active

          oldest

          votes


















          0












          $begingroup$

          If I understand correctly, you want to (1) subtract every row of A from every row of B, and then (2) sum all values within each resulting row.





          1. For step (1) you can use a single bsxfun, but you need to permute dimensions:



            subtractedRows = bsxfun(@minus, B, permute(A,[3 2 1]));


            For any m, n, subtractedRows(m,:,n) gives B(m,:)-A(n,:).




          2. For step (2) you only need to sum along second dimension:



            result = squeeze(sum(subtractedRows,2));







          share|cite|improve this answer









          $endgroup$





















            0












            $begingroup$

            It sounds like you can sum across the columns of A and B before subtracting and simplify the computation.



            If broadcasting is available (newest versions of Matlab and Octave):



            C = sum(B, 2) - sum(A, 2)'


            If not:



            C = repmat(sum(B, 2), [1 rows(A)]) - repmat(sum(A, 2)', [rows(B) 1])





            share|cite|improve this answer









            $endgroup$













              Your Answer





              StackExchange.ifUsing("editor", function () {
              return StackExchange.using("mathjaxEditing", function () {
              StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
              StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
              });
              });
              }, "mathjax-editing");

              StackExchange.ready(function() {
              var channelOptions = {
              tags: "".split(" "),
              id: "69"
              };
              initTagRenderer("".split(" "), "".split(" "), channelOptions);

              StackExchange.using("externalEditor", function() {
              // Have to fire editor after snippets, if snippets enabled
              if (StackExchange.settings.snippets.snippetsEnabled) {
              StackExchange.using("snippets", function() {
              createEditor();
              });
              }
              else {
              createEditor();
              }
              });

              function createEditor() {
              StackExchange.prepareEditor({
              heartbeatType: 'answer',
              autoActivateHeartbeat: false,
              convertImagesToLinks: true,
              noModals: true,
              showLowRepImageUploadWarning: true,
              reputationToPostImages: 10,
              bindNavPrevention: true,
              postfix: "",
              imageUploader: {
              brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
              contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
              allowUrls: true
              },
              noCode: true, onDemand: true,
              discardSelector: ".discard-answer"
              ,immediatelyShowMarkdownHelp:true
              });


              }
              });














              draft saved

              draft discarded


















              StackExchange.ready(
              function () {
              StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f708452%2fsubtract-matrices-of-different-sizes%23new-answer', 'question_page');
              }
              );

              Post as a guest















              Required, but never shown

























              2 Answers
              2






              active

              oldest

              votes








              2 Answers
              2






              active

              oldest

              votes









              active

              oldest

              votes






              active

              oldest

              votes









              0












              $begingroup$

              If I understand correctly, you want to (1) subtract every row of A from every row of B, and then (2) sum all values within each resulting row.





              1. For step (1) you can use a single bsxfun, but you need to permute dimensions:



                subtractedRows = bsxfun(@minus, B, permute(A,[3 2 1]));


                For any m, n, subtractedRows(m,:,n) gives B(m,:)-A(n,:).




              2. For step (2) you only need to sum along second dimension:



                result = squeeze(sum(subtractedRows,2));







              share|cite|improve this answer









              $endgroup$


















                0












                $begingroup$

                If I understand correctly, you want to (1) subtract every row of A from every row of B, and then (2) sum all values within each resulting row.





                1. For step (1) you can use a single bsxfun, but you need to permute dimensions:



                  subtractedRows = bsxfun(@minus, B, permute(A,[3 2 1]));


                  For any m, n, subtractedRows(m,:,n) gives B(m,:)-A(n,:).




                2. For step (2) you only need to sum along second dimension:



                  result = squeeze(sum(subtractedRows,2));







                share|cite|improve this answer









                $endgroup$
















                  0












                  0








                  0





                  $begingroup$

                  If I understand correctly, you want to (1) subtract every row of A from every row of B, and then (2) sum all values within each resulting row.





                  1. For step (1) you can use a single bsxfun, but you need to permute dimensions:



                    subtractedRows = bsxfun(@minus, B, permute(A,[3 2 1]));


                    For any m, n, subtractedRows(m,:,n) gives B(m,:)-A(n,:).




                  2. For step (2) you only need to sum along second dimension:



                    result = squeeze(sum(subtractedRows,2));







                  share|cite|improve this answer









                  $endgroup$



                  If I understand correctly, you want to (1) subtract every row of A from every row of B, and then (2) sum all values within each resulting row.





                  1. For step (1) you can use a single bsxfun, but you need to permute dimensions:



                    subtractedRows = bsxfun(@minus, B, permute(A,[3 2 1]));


                    For any m, n, subtractedRows(m,:,n) gives B(m,:)-A(n,:).




                  2. For step (2) you only need to sum along second dimension:



                    result = squeeze(sum(subtractedRows,2));








                  share|cite|improve this answer












                  share|cite|improve this answer



                  share|cite|improve this answer










                  answered Mar 14 '14 at 11:13









                  Luis MendoLuis Mendo

                  1,3321821




                  1,3321821























                      0












                      $begingroup$

                      It sounds like you can sum across the columns of A and B before subtracting and simplify the computation.



                      If broadcasting is available (newest versions of Matlab and Octave):



                      C = sum(B, 2) - sum(A, 2)'


                      If not:



                      C = repmat(sum(B, 2), [1 rows(A)]) - repmat(sum(A, 2)', [rows(B) 1])





                      share|cite|improve this answer









                      $endgroup$


















                        0












                        $begingroup$

                        It sounds like you can sum across the columns of A and B before subtracting and simplify the computation.



                        If broadcasting is available (newest versions of Matlab and Octave):



                        C = sum(B, 2) - sum(A, 2)'


                        If not:



                        C = repmat(sum(B, 2), [1 rows(A)]) - repmat(sum(A, 2)', [rows(B) 1])





                        share|cite|improve this answer









                        $endgroup$
















                          0












                          0








                          0





                          $begingroup$

                          It sounds like you can sum across the columns of A and B before subtracting and simplify the computation.



                          If broadcasting is available (newest versions of Matlab and Octave):



                          C = sum(B, 2) - sum(A, 2)'


                          If not:



                          C = repmat(sum(B, 2), [1 rows(A)]) - repmat(sum(A, 2)', [rows(B) 1])





                          share|cite|improve this answer









                          $endgroup$



                          It sounds like you can sum across the columns of A and B before subtracting and simplify the computation.



                          If broadcasting is available (newest versions of Matlab and Octave):



                          C = sum(B, 2) - sum(A, 2)'


                          If not:



                          C = repmat(sum(B, 2), [1 rows(A)]) - repmat(sum(A, 2)', [rows(B) 1])






                          share|cite|improve this answer












                          share|cite|improve this answer



                          share|cite|improve this answer










                          answered Sep 15 '16 at 14:44









                          NirNir

                          65837




                          65837






























                              draft saved

                              draft discarded




















































                              Thanks for contributing an answer to Mathematics Stack Exchange!


                              • Please be sure to answer the question. Provide details and share your research!

                              But avoid



                              • Asking for help, clarification, or responding to other answers.

                              • Making statements based on opinion; back them up with references or personal experience.


                              Use MathJax to format equations. MathJax reference.


                              To learn more, see our tips on writing great answers.




                              draft saved


                              draft discarded














                              StackExchange.ready(
                              function () {
                              StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f708452%2fsubtract-matrices-of-different-sizes%23new-answer', 'question_page');
                              }
                              );

                              Post as a guest















                              Required, but never shown





















































                              Required, but never shown














                              Required, but never shown












                              Required, but never shown







                              Required, but never shown

































                              Required, but never shown














                              Required, but never shown












                              Required, but never shown







                              Required, but never shown







                              Popular posts from this blog

                              Bressuire

                              Cabo Verde

                              Gyllenstierna