Fastest way to go from linear index to grid index












3














I'm sure this has been asked before but I'm interested in going from a position in a vector to the index in the grid version of the vector with given strides, for example, say I have the vector:



vec = {58, 94, 19, 68, 54, 77, 1, 18, 49, 20, 90, 44, 91, 89, 15, 0, 

   60, 18, 19, 44, 87, 5, 8, 42, 51, 55, 87, 71, 83, 68, 53, 58, 27, 

   17, 8, 14, 33, 58, 86, 3, 91, 66, 3, 16, 98, 84, 72, 98, 9, 30, 90,

    99, 15, 0, 82, 76, 86, 58, 77, 58};


And say I have strides {5, 4, 3}, the position 35 in the vector would correspond to the index {3, 4, 2}:



vec[[35]]



4



ArrayReshape[vec, {5, 4, 3}][[3, 4, 2]]



4


How can I get this index fast and in a vectorized fashion because I will have potentially many positions to extract?






list-manipulation performance-tuning







share|improve this question



























    3














    I'm sure this has been asked before but I'm interested in going from a position in a vector to the index in the grid version of the vector with given strides, for example, say I have the vector:



    vec = {58, 94, 19, 68, 54, 77, 1, 18, 49, 20, 90, 44, 91, 89, 15, 0, 
    
   60, 18, 19, 44, 87, 5, 8, 42, 51, 55, 87, 71, 83, 68, 53, 58, 27, 
    
   17, 8, 14, 33, 58, 86, 3, 91, 66, 3, 16, 98, 84, 72, 98, 9, 30, 90,
    
    99, 15, 0, 82, 76, 86, 58, 77, 58};


    And say I have strides {5, 4, 3}, the position 35 in the vector would correspond to the index {3, 4, 2}:



    vec[[35]]
    

    
4
    

    
ArrayReshape[vec, {5, 4, 3}][[3, 4, 2]]
    

    
4


    How can I get this index fast and in a vectorized fashion because I will have potentially many positions to extract?






    list-manipulation performance-tuning







    share|improve this question

























      3












      3








      3







      I'm sure this has been asked before but I'm interested in going from a position in a vector to the index in the grid version of the vector with given strides, for example, say I have the vector:



      vec = {58, 94, 19, 68, 54, 77, 1, 18, 49, 20, 90, 44, 91, 89, 15, 0, 
      
   60, 18, 19, 44, 87, 5, 8, 42, 51, 55, 87, 71, 83, 68, 53, 58, 27, 
      
   17, 8, 14, 33, 58, 86, 3, 91, 66, 3, 16, 98, 84, 72, 98, 9, 30, 90,
      
    99, 15, 0, 82, 76, 86, 58, 77, 58};


      And say I have strides {5, 4, 3}, the position 35 in the vector would correspond to the index {3, 4, 2}:



      vec[[35]]
      

      
4
      

      
ArrayReshape[vec, {5, 4, 3}][[3, 4, 2]]
      

      
4


      How can I get this index fast and in a vectorized fashion because I will have potentially many positions to extract?






      list-manipulation performance-tuning







      share|improve this question













      I'm sure this has been asked before but I'm interested in going from a position in a vector to the index in the grid version of the vector with given strides, for example, say I have the vector:



      vec = {58, 94, 19, 68, 54, 77, 1, 18, 49, 20, 90, 44, 91, 89, 15, 0, 
      
   60, 18, 19, 44, 87, 5, 8, 42, 51, 55, 87, 71, 83, 68, 53, 58, 27, 
      
   17, 8, 14, 33, 58, 86, 3, 91, 66, 3, 16, 98, 84, 72, 98, 9, 30, 90,
      
    99, 15, 0, 82, 76, 86, 58, 77, 58};


      And say I have strides {5, 4, 3}, the position 35 in the vector would correspond to the index {3, 4, 2}:



      vec[[35]]
      

      
4
      

      
ArrayReshape[vec, {5, 4, 3}][[3, 4, 2]]
      

      
4


      How can I get this index fast and in a vectorized fashion because I will have potentially many positions to extract?






      list-manipulation performance-tuning




      list-manipulation performance-tuning






      share|improve this question













      share|improve this question











      share|improve this question




      share|improve this question










      asked 1 hour ago









      b3m2a1

      26.7k257154




      26.7k257154






















          2 Answers
          2






          active

          oldest

          votes


















          2














          I think this is what you want:



          IntegerDigits[35 - 1, MixedRadix[{5, 4, 3}], 3] + 1


          In general:



          gridIndex[n_Integer, shape_List] := 
          
 IntegerDigits[n - 1, MixedRadix[shape], Length@shape] + 1





          share|improve this answer



















          • 1




            @C.E. You're right, it just needs dimension length specification
            – swish
            45 mins ago






          • 1




            This is a very nice solution, +1
            – C. E.
            43 mins ago



















          1














          Here's what I came up with:



          getSubindex[index_, stride_] := {
          
  Mod[index, stride, 1],
          
  Ceiling[index/stride]
          
  }
          
getIndex[index_, strides_] := 
          
 Reverse@FoldPairList[getSubindex, index, Reverse@strides]


          This is comparable to swish's solution speed-wise:



          gridIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



          {0.000061, {3, 2, 3, 4}}




          getIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



          {0.000052, {3, 2, 3, 4}}







          share|improve this answer























            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: "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%2f188806%2ffastest-way-to-go-from-linear-index-to-grid-index%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









            2














            I think this is what you want:



            IntegerDigits[35 - 1, MixedRadix[{5, 4, 3}], 3] + 1


            In general:



            gridIndex[n_Integer, shape_List] := 
            
 IntegerDigits[n - 1, MixedRadix[shape], Length@shape] + 1





            share|improve this answer



















            • 1




              @C.E. You're right, it just needs dimension length specification
              – swish
              45 mins ago






            • 1




              This is a very nice solution, +1
              – C. E.
              43 mins ago
















            2














            I think this is what you want:



            IntegerDigits[35 - 1, MixedRadix[{5, 4, 3}], 3] + 1


            In general:



            gridIndex[n_Integer, shape_List] := 
            
 IntegerDigits[n - 1, MixedRadix[shape], Length@shape] + 1





            share|improve this answer



















            • 1




              @C.E. You're right, it just needs dimension length specification
              – swish
              45 mins ago






            • 1




              This is a very nice solution, +1
              – C. E.
              43 mins ago














            2












            2








            2






            I think this is what you want:



            IntegerDigits[35 - 1, MixedRadix[{5, 4, 3}], 3] + 1


            In general:



            gridIndex[n_Integer, shape_List] := 
            
 IntegerDigits[n - 1, MixedRadix[shape], Length@shape] + 1





            share|improve this answer














            I think this is what you want:



            IntegerDigits[35 - 1, MixedRadix[{5, 4, 3}], 3] + 1


            In general:



            gridIndex[n_Integer, shape_List] := 
            
 IntegerDigits[n - 1, MixedRadix[shape], Length@shape] + 1






            share|improve this answer














            share|improve this answer



            share|improve this answer








            edited 42 mins ago

























            answered 53 mins ago









            swish

            3,9611534




            3,9611534








            • 1




              @C.E. You're right, it just needs dimension length specification
              – swish
              45 mins ago






            • 1




              This is a very nice solution, +1
              – C. E.
              43 mins ago














            • 1




              @C.E. You're right, it just needs dimension length specification
              – swish
              45 mins ago






            • 1




              This is a very nice solution, +1
              – C. E.
              43 mins ago








            1




            1




            @C.E. You're right, it just needs dimension length specification
            – swish
            45 mins ago




            @C.E. You're right, it just needs dimension length specification
            – swish
            45 mins ago




            1




            1




            This is a very nice solution, +1
            – C. E.
            43 mins ago




            This is a very nice solution, +1
            – C. E.
            43 mins ago











            1














            Here's what I came up with:



            getSubindex[index_, stride_] := {
            
  Mod[index, stride, 1],
            
  Ceiling[index/stride]
            
  }
            
getIndex[index_, strides_] := 
            
 Reverse@FoldPairList[getSubindex, index, Reverse@strides]


            This is comparable to swish's solution speed-wise:



            gridIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



            {0.000061, {3, 2, 3, 4}}




            getIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



            {0.000052, {3, 2, 3, 4}}







            share|improve this answer




























              1














              Here's what I came up with:



              getSubindex[index_, stride_] := {
              
  Mod[index, stride, 1],
              
  Ceiling[index/stride]
              
  }
              
getIndex[index_, strides_] := 
              
 Reverse@FoldPairList[getSubindex, index, Reverse@strides]


              This is comparable to swish's solution speed-wise:



              gridIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



              {0.000061, {3, 2, 3, 4}}




              getIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



              {0.000052, {3, 2, 3, 4}}







              share|improve this answer


























                1












                1








                1






                Here's what I came up with:



                getSubindex[index_, stride_] := {
                
  Mod[index, stride, 1],
                
  Ceiling[index/stride]
                
  }
                
getIndex[index_, strides_] := 
                
 Reverse@FoldPairList[getSubindex, index, Reverse@strides]


                This is comparable to swish's solution speed-wise:



                gridIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



                {0.000061, {3, 2, 3, 4}}




                getIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



                {0.000052, {3, 2, 3, 4}}







                share|improve this answer














                Here's what I came up with:



                getSubindex[index_, stride_] := {
                
  Mod[index, stride, 1],
                
  Ceiling[index/stride]
                
  }
                
getIndex[index_, strides_] := 
                
 Reverse@FoldPairList[getSubindex, index, Reverse@strides]


                This is comparable to swish's solution speed-wise:



                gridIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



                {0.000061, {3, 2, 3, 4}}




                getIndex[1000, {3, 5, 4, 6}] // RepeatedTiming



                {0.000052, {3, 2, 3, 4}}








                share|improve this answer














                share|improve this answer



                share|improve this answer








                edited 24 mins ago

























                answered 45 mins ago









                C. E.

                49.9k397202




                49.9k397202






























                    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.





                    Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


                    Please pay close attention to the following guidance:


                    • 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.


                    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%2f188806%2ffastest-way-to-go-from-linear-index-to-grid-index%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

                    Eastern Orthodox Church

                    Image
                    Not to be confused with Oriental Orthodoxy, Eastern Christianity, or Eastern Catholic Churches. "Eastern Orthodoxy", "Orthodox Church", "Orthodox Catholic Church", and "Orthodox Christian Church" redirect here. For other uses, see Orthodox (disambiguation). This article may be too long to read and navigate comfortably . The readable prose size is 132 kilobytes. Please consider splitting content into sub-articles, condensing it, or adding or removing subheadings. ( January 2018 ) Part of a series on the Eastern Orthodox Church Mosaic of Christ Pantocrator, Hagia Sophia Overview Structure Theology (History of theology) Liturgy Church history Holy Mysteries View of salvation View of Mary View of icons Background Crucifixion / Resurrection / Ascension of Jesus Christianity Christian Church Apostolic succession Four Marks of the Church Orthodoxy Organization Autocephaly Patriarchate ...
                    Read more

                    Zagreb

                    Image
                    This article is about the Croatian capital city. For other uses, see Zagreb (disambiguation). City in City of Zagreb, Croatia Zagreb City Grad Zagreb City of Zagreb Clockwise, from top: St. Mark's Square, Croatian State Archives, Zagreb Tram, Art Pavilion, Cibona & HOTO towers and Croatian National Theatre. Flag Coat of arms Zagreb Location of Zagreb in Croatia Show map of Croatia Zagreb Zagreb (Europe) Show map of Europe Coordinates: 45°49′N 15°59′E  /  45.817°N 15.983°E  / 45.817; 15.983 Coordinates: 45°49′N 15°59′E  /  45.817°N 15.983°E  / 45.817; 15.983 Country   Croatia County City of Zagreb RC diocese 1094 Free royal city 1242 Unified 1850 Subdivisions 17 city districts 218 local committees (70 settlements) Government  • Type Mayor-Council  • Mayor Milan Bandić (BM 365)  • City Assembly 51 members BM 365, ZL, NS-R (14) SDP, GLAS, HNS, HSS, NH-PS (13) ...
                    Read more

                    Understanding the information contained in the Deep Space Network XML data?

                    Image
                    up vote 1 down vote favorite NASA's Deep Space Network activity can be viewed conveniently at https://eyes.nasa.gov/dsn/dsn.html but if you would like data in a more raw XML form, it can be accessed as described in this answer Below is a small Python script that reads the XML into a list of dictionaries, and saves it to disk as a JSON. The dictionaries are organized at the higher level by dish, then by spacecraft per dish. Question: What are the most-likely meanings of all the data field contained in the Deep Space Network XML data? This is the example url from the other answer, used in the script below: https://eyes.nasa.gov/dsn/data/1365107113.xml def dictify(r,root=True): # https://stackoverflow.com/questions/2148119/how-to-convert-an-xml-string-to-a-dictionary-in-python/30923963#30923963 if root: ret...
                    Read more