Apply MapThread to all but one variable












4












$begingroup$


I would like to know what is the most efficient to implement the following computation. Given three lists



    a = {a_1,a_2, a_3, …, a_n}
b = {b_1,b_2, b_3, …, b_n}
c = {c_1,c_2, c_3, …, c_n}


and a function $f(x_1,x_2,x_3)$, obtain



     f(a_1,b_1,c_1)   f(a_1,b_1,c_2)   .....   f(a_1,b_1,c_n)  
f(a_2,b_2,c_1) f(a_2,b_2,c_2) ..... f(a_2,b_2,c_n)
..... ..... ..... .....
f(a_n,b_n,c_1) f(a_n,b_n,c_2) ..... f(a_n,b_n,c_n)


I cannot find a solution not using For.










share|improve this question











$endgroup$

















    4












    $begingroup$


    I would like to know what is the most efficient to implement the following computation. Given three lists



        a = {a_1,a_2, a_3, …, a_n}
    b = {b_1,b_2, b_3, …, b_n}
    c = {c_1,c_2, c_3, …, c_n}


    and a function $f(x_1,x_2,x_3)$, obtain



         f(a_1,b_1,c_1)   f(a_1,b_1,c_2)   .....   f(a_1,b_1,c_n)  
    f(a_2,b_2,c_1) f(a_2,b_2,c_2) ..... f(a_2,b_2,c_n)
    ..... ..... ..... .....
    f(a_n,b_n,c_1) f(a_n,b_n,c_2) ..... f(a_n,b_n,c_n)


    I cannot find a solution not using For.










    share|improve this question











    $endgroup$















      4












      4








      4





      $begingroup$


      I would like to know what is the most efficient to implement the following computation. Given three lists



          a = {a_1,a_2, a_3, …, a_n}
      b = {b_1,b_2, b_3, …, b_n}
      c = {c_1,c_2, c_3, …, c_n}


      and a function $f(x_1,x_2,x_3)$, obtain



           f(a_1,b_1,c_1)   f(a_1,b_1,c_2)   .....   f(a_1,b_1,c_n)  
      f(a_2,b_2,c_1) f(a_2,b_2,c_2) ..... f(a_2,b_2,c_n)
      ..... ..... ..... .....
      f(a_n,b_n,c_1) f(a_n,b_n,c_2) ..... f(a_n,b_n,c_n)


      I cannot find a solution not using For.










      share|improve this question











      $endgroup$




      I would like to know what is the most efficient to implement the following computation. Given three lists



          a = {a_1,a_2, a_3, …, a_n}
      b = {b_1,b_2, b_3, …, b_n}
      c = {c_1,c_2, c_3, …, c_n}


      and a function $f(x_1,x_2,x_3)$, obtain



           f(a_1,b_1,c_1)   f(a_1,b_1,c_2)   .....   f(a_1,b_1,c_n)  
      f(a_2,b_2,c_1) f(a_2,b_2,c_2) ..... f(a_2,b_2,c_n)
      ..... ..... ..... .....
      f(a_n,b_n,c_1) f(a_n,b_n,c_2) ..... f(a_n,b_n,c_n)


      I cannot find a solution not using For.







      list-manipulation






      share|improve this question















      share|improve this question













      share|improve this question




      share|improve this question








      edited 7 hours ago









      corey979

      20.9k64382




      20.9k64382










      asked 7 hours ago









      SmerdjakovSmerdjakov

      1405




      1405






















          4 Answers
          4






          active

          oldest

          votes


















          6












          $begingroup$

          Here's one way to do it with Outer:



          n = 3;
          l1 = Array[a, n];
          l2 = Array[b, n];
          l3 = Array[c, n];

          Outer[
          f[#1[[1]], #1[[2]], #2] &,
          Transpose @ {l1, l2},
          l3,
          1
          ]



          Out[25]= {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
          f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
          f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
          f[a[3], b[3], c[3]]}}







          share|improve this answer









          $endgroup$









          • 1




            $begingroup$
            Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
            $endgroup$
            – Roman
            7 hours ago



















          5












          $begingroup$

          a = {a1, a2, a3, a4, a5};
          b = {b1, b2, b3, b4, b5};
          c = {c1, c2, c3, c4, c5};

          Table[f[a[[j]], b[[j]], c[[k]]], {j, 1, 5}, {k, 1, 5}]


          enter image description here






          share|improve this answer









          $endgroup$





















            3












            $begingroup$

            Another possibility is to use the 3-arg version of Thread. With Sjoerd's example:



            n = 3;
            l1 = Array[a,n];
            l2 = Array[b,n];
            l3 = Array[c,n];


            Using Thread:



            Thread /@ Thread[f[l1, l2, l3], List, 2]



            {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
            f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
            f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
            f[a[3], b[3], c[3]]}}







            share|improve this answer









            $endgroup$





















              1












              $begingroup$

              Another way with Curry and Through.



              Through /@ Apply[Curry[f, {2, 3, 1}] /@ c] /@ Transpose@{a, b}



              Mathematica graphics




              Hope this helps.






              share|improve this answer









              $endgroup$














                Your Answer








                StackExchange.ready(function() {
                var channelOptions = {
                tags: "".split(" "),
                id: "387"
                };
                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: false,
                noModals: true,
                showLowRepImageUploadWarning: true,
                reputationToPostImages: null,
                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
                },
                onDemand: true,
                discardSelector: ".discard-answer"
                ,immediatelyShowMarkdownHelp:true
                });


                }
                });














                draft saved

                draft discarded


















                StackExchange.ready(
                function () {
                StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f197144%2fapply-mapthread-to-all-but-one-variable%23new-answer', 'question_page');
                }
                );

                Post as a guest















                Required, but never shown

























                4 Answers
                4






                active

                oldest

                votes








                4 Answers
                4






                active

                oldest

                votes









                active

                oldest

                votes






                active

                oldest

                votes









                6












                $begingroup$

                Here's one way to do it with Outer:



                n = 3;
                l1 = Array[a, n];
                l2 = Array[b, n];
                l3 = Array[c, n];

                Outer[
                f[#1[[1]], #1[[2]], #2] &,
                Transpose @ {l1, l2},
                l3,
                1
                ]



                Out[25]= {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                f[a[3], b[3], c[3]]}}







                share|improve this answer









                $endgroup$









                • 1




                  $begingroup$
                  Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
                  $endgroup$
                  – Roman
                  7 hours ago
















                6












                $begingroup$

                Here's one way to do it with Outer:



                n = 3;
                l1 = Array[a, n];
                l2 = Array[b, n];
                l3 = Array[c, n];

                Outer[
                f[#1[[1]], #1[[2]], #2] &,
                Transpose @ {l1, l2},
                l3,
                1
                ]



                Out[25]= {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                f[a[3], b[3], c[3]]}}







                share|improve this answer









                $endgroup$









                • 1




                  $begingroup$
                  Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
                  $endgroup$
                  – Roman
                  7 hours ago














                6












                6








                6





                $begingroup$

                Here's one way to do it with Outer:



                n = 3;
                l1 = Array[a, n];
                l2 = Array[b, n];
                l3 = Array[c, n];

                Outer[
                f[#1[[1]], #1[[2]], #2] &,
                Transpose @ {l1, l2},
                l3,
                1
                ]



                Out[25]= {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                f[a[3], b[3], c[3]]}}







                share|improve this answer









                $endgroup$



                Here's one way to do it with Outer:



                n = 3;
                l1 = Array[a, n];
                l2 = Array[b, n];
                l3 = Array[c, n];

                Outer[
                f[#1[[1]], #1[[2]], #2] &,
                Transpose @ {l1, l2},
                l3,
                1
                ]



                Out[25]= {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                f[a[3], b[3], c[3]]}}








                share|improve this answer












                share|improve this answer



                share|improve this answer










                answered 7 hours ago









                Sjoerd SmitSjoerd Smit

                4,630817




                4,630817








                • 1




                  $begingroup$
                  Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
                  $endgroup$
                  – Roman
                  7 hours ago














                • 1




                  $begingroup$
                  Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
                  $endgroup$
                  – Roman
                  7 hours ago








                1




                1




                $begingroup$
                Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
                $endgroup$
                – Roman
                7 hours ago




                $begingroup$
                Or Outer[f[Sequence @@ #1, #2] &, Transpose@{l1, l2}, l3, 1] so you don't need to unravel #1 manually.
                $endgroup$
                – Roman
                7 hours ago











                5












                $begingroup$

                a = {a1, a2, a3, a4, a5};
                b = {b1, b2, b3, b4, b5};
                c = {c1, c2, c3, c4, c5};

                Table[f[a[[j]], b[[j]], c[[k]]], {j, 1, 5}, {k, 1, 5}]


                enter image description here






                share|improve this answer









                $endgroup$


















                  5












                  $begingroup$

                  a = {a1, a2, a3, a4, a5};
                  b = {b1, b2, b3, b4, b5};
                  c = {c1, c2, c3, c4, c5};

                  Table[f[a[[j]], b[[j]], c[[k]]], {j, 1, 5}, {k, 1, 5}]


                  enter image description here






                  share|improve this answer









                  $endgroup$
















                    5












                    5








                    5





                    $begingroup$

                    a = {a1, a2, a3, a4, a5};
                    b = {b1, b2, b3, b4, b5};
                    c = {c1, c2, c3, c4, c5};

                    Table[f[a[[j]], b[[j]], c[[k]]], {j, 1, 5}, {k, 1, 5}]


                    enter image description here






                    share|improve this answer









                    $endgroup$



                    a = {a1, a2, a3, a4, a5};
                    b = {b1, b2, b3, b4, b5};
                    c = {c1, c2, c3, c4, c5};

                    Table[f[a[[j]], b[[j]], c[[k]]], {j, 1, 5}, {k, 1, 5}]


                    enter image description here







                    share|improve this answer












                    share|improve this answer



                    share|improve this answer










                    answered 7 hours ago









                    corey979corey979

                    20.9k64382




                    20.9k64382























                        3












                        $begingroup$

                        Another possibility is to use the 3-arg version of Thread. With Sjoerd's example:



                        n = 3;
                        l1 = Array[a,n];
                        l2 = Array[b,n];
                        l3 = Array[c,n];


                        Using Thread:



                        Thread /@ Thread[f[l1, l2, l3], List, 2]



                        {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                        f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                        f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                        f[a[3], b[3], c[3]]}}







                        share|improve this answer









                        $endgroup$


















                          3












                          $begingroup$

                          Another possibility is to use the 3-arg version of Thread. With Sjoerd's example:



                          n = 3;
                          l1 = Array[a,n];
                          l2 = Array[b,n];
                          l3 = Array[c,n];


                          Using Thread:



                          Thread /@ Thread[f[l1, l2, l3], List, 2]



                          {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                          f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                          f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                          f[a[3], b[3], c[3]]}}







                          share|improve this answer









                          $endgroup$
















                            3












                            3








                            3





                            $begingroup$

                            Another possibility is to use the 3-arg version of Thread. With Sjoerd's example:



                            n = 3;
                            l1 = Array[a,n];
                            l2 = Array[b,n];
                            l3 = Array[c,n];


                            Using Thread:



                            Thread /@ Thread[f[l1, l2, l3], List, 2]



                            {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                            f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                            f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                            f[a[3], b[3], c[3]]}}







                            share|improve this answer









                            $endgroup$



                            Another possibility is to use the 3-arg version of Thread. With Sjoerd's example:



                            n = 3;
                            l1 = Array[a,n];
                            l2 = Array[b,n];
                            l3 = Array[c,n];


                            Using Thread:



                            Thread /@ Thread[f[l1, l2, l3], List, 2]



                            {{f[a[1], b[1], c[1]], f[a[1], b[1], c[2]],
                            f[a[1], b[1], c[3]]}, {f[a[2], b[2], c[1]], f[a[2], b[2], c[2]],
                            f[a[2], b[2], c[3]]}, {f[a[3], b[3], c[1]], f[a[3], b[3], c[2]],
                            f[a[3], b[3], c[3]]}}








                            share|improve this answer












                            share|improve this answer



                            share|improve this answer










                            answered 6 hours ago









                            Carl WollCarl Woll

                            75.9k3100198




                            75.9k3100198























                                1












                                $begingroup$

                                Another way with Curry and Through.



                                Through /@ Apply[Curry[f, {2, 3, 1}] /@ c] /@ Transpose@{a, b}



                                Mathematica graphics




                                Hope this helps.






                                share|improve this answer









                                $endgroup$


















                                  1












                                  $begingroup$

                                  Another way with Curry and Through.



                                  Through /@ Apply[Curry[f, {2, 3, 1}] /@ c] /@ Transpose@{a, b}



                                  Mathematica graphics




                                  Hope this helps.






                                  share|improve this answer









                                  $endgroup$
















                                    1












                                    1








                                    1





                                    $begingroup$

                                    Another way with Curry and Through.



                                    Through /@ Apply[Curry[f, {2, 3, 1}] /@ c] /@ Transpose@{a, b}



                                    Mathematica graphics




                                    Hope this helps.






                                    share|improve this answer









                                    $endgroup$



                                    Another way with Curry and Through.



                                    Through /@ Apply[Curry[f, {2, 3, 1}] /@ c] /@ Transpose@{a, b}



                                    Mathematica graphics




                                    Hope this helps.







                                    share|improve this answer












                                    share|improve this answer



                                    share|improve this answer










                                    answered 4 hours ago









                                    EdmundEdmund

                                    26.8k330103




                                    26.8k330103






























                                        draft saved

                                        draft discarded




















































                                        Thanks for contributing an answer to Mathematica 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%2fmathematica.stackexchange.com%2fquestions%2f197144%2fapply-mapthread-to-all-but-one-variable%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

                                        Bundesstraße 106

                                        Le Mesnil-Réaume

                                        Ida-Boy-Ed-Garten