Skip to main content

list manipulation - Adding points limits


With 4 points I set a path where your total length was subdivided into 300 points


Clear["Global`*"]  
points={{0,0},{5,8.07774721},{10,4.24499363},{20,9.28880172}}  
Needs["GraphUtilities`"]  

n=300;  
pts=LineScaledCoordinate[points,N@(#/n)]&/@Range[0,n];


However with this code were not included in the list points. Then with the function below I've added.


pts=Insert[pts,points[[2]],107];  
pts=Insert[pts,points[[3]],177];

I added these points in the correct position because I searched by hand and applied the function.
This is the first question:
How do you find these positions in an automatic way?


With that done, with these new items i created divisions where these points are at the end of each list.
How could I do that these points are reference for these divisions?



list1=pts[[#,1]]&/@Range[107];  
list2=pts[[#,1]]&/@Range[108,177];  
list3=pts[[#,1]]&/@Range[178,303];

Answer



This is a answer based on the solution of Michael E2.


I would like to offer my attempts without damaging the response that does not belong to me.


My attempts


Here is the list pts without the values from the list points


Clear["Global`*"]
points = {{0, 0}, {5, 8.07774721}, {10, 4.24499363}, {20, 9.28880172}};

Needs["GraphUtilities`"]

n = 300;
pts = LineScaledCoordinate[points, N@(#/n)] & /@ Range[0, n];

With the function below it is possible to structure the points from the list points in pairs with offset value 1.


Partition[points, 2, 1]


$\left( \begin{array}{cc} \{0,0\} & \{5,8.07775\} \\ \{5,8.07775\} & \{10,4.24499\} \\ \{10,4.24499\} & \{20,9.2888\} \\ \end{array} \right)$




With the function below it is possible to get the lengths between the pairs mentioned above


partialPath = EuclideanDistance @@@ Partition[points, 2, 1] // N

enter image description here



$\{9.5,6.3,11.2\}$



With the function below it is possible to make the accumulated form sum


Accumulate[partialPath]



$\{9.5,15.8,27.\}$



With the function below it is possible to delete the last element of the accumulated numbers list


Most@Accumulate[partialPath]


$\{9.5,15.8\}$




With the function below it is possible to use the last element of the accumulated numbers list


Last@Accumulate[partialPath]


$27$



With the function below it is possible to get the percentage of each path with the full path


porcentPath = {#/Last@Accumulate[partialPath]} & /@Most@Accumulate[partialPath]

enter image description here




$\left( \begin{array}{c} 0.351852 \\ 0.585185 \\ \end{array} \right)$



With the function below, it is possible to separate the points in the points list that are not included in the pts list


valuesPositions = Rest@Most@points

enter image description here



$\left( \begin{array}{cc} 5 & 8.07775 \\ 10 & 4.24499 \\ \end{array} \right)$




With the function below it is possible to identify which positions in should be inserted the values {5,8.07774721} and {10,4.24499363}


pos = Flatten@Transpose[Ceiling[n*porcentPath] + 1]

enter image description here



$\{107,177\}$



With the function below it is possible to define a list to be added values, the values to be added and the positions they should be inserted.


InsertValuePosition[list_, values_, positions_] := 
 Fold[Insert[#, #2[[1]], #2[[2]]] &, list, Transpose[{values, positions}]]

pts = InsertValuePosition[pts, valuesPositions, Transpose[{pos}]]

With the function below it is possible to subdivide the pts list in three other lists where the entered values are at the end of the two first lists


FoldPairList[TakeDrop, pts, Flatten@{0, pos, Length@pts} // Differences]

enter image description here


Comments

Post a Comment

Popular posts from this blog

functions - Get leading series expansion term?

Given a function f[x] , I would like to have a function leadingSeries that returns just the leading term in the series around x=0 . For example: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x)] x and leadingSeries[(1/x + 2 + (1 - 1/x^3)/4)/(4 + x)] -(1/(16 x^3)) Is there such a function in Mathematica? Or maybe one can implement it efficiently? EDIT I finally went with the following implementation, based on Carl Woll 's answer: lds[ex_,x_]:=( (ex/.x->(x+O[x]^2))/.SeriesData[U_,Z_,L_List,Mi_,Ma_,De_]:>SeriesData[U,Z,{L[[1]]},Mi,Mi+1,De]//Quiet//Normal) The advantage is, that this one also properly works with functions whose leading term is a constant: lds[Exp[x],x] 1 Answer Update 1 Updated to eliminate SeriesData and to not return additional terms Perhaps you could use: leadingSeries[expr_, x_] := Normal[expr /. x->(x+O[x]^2) /. a_List :> Take[a, 1]] Then for your examples: leadingSeries[(1/x + 2)/(4 + 1/x^2 + x), x] leadingSeries[Exp[x], x] leadingSeries[(1/x + 2 + (1 - 1/x...
Post a Comment
Read more

How to thread a list

I have data in format data = {{a1, a2}, {b1, b2}, {c1, c2}, {d1, d2}} Tableform: I want to thread it to : tdata = {{{a1, b1}, {a2, b2}}, {{a1, c1}, {a2, c2}}, {{a1, d1}, {a2, d2}}} Tableform: And I would like to do better then pseudofunction[n_] := Transpose[{data2[[1]], data2[[n]]}]; SetAttributes[pseudofunction, Listable]; Range[2, 4] // pseudofunction Here is my benchmark data, where data3 is normal sample of real data. data3 = Drop[ExcelWorkBook[[Column1 ;; Column4]], None, 1]; data2 = {a #, b #, c #, d #} & /@ Range[1, 10^5]; data = RandomReal[{0, 1}, {10^6, 4}]; Here is my benchmark code kptnw[list_] := Transpose[{Table[First@#, {Length@# - 1}], Rest@#}, {3, 1, 2}] &@list kptnw2[list_] := Transpose[{ConstantArray[First@#, Length@# - 1], Rest@#}, {3, 1, 2}] &@list OleksandrR[list_] := Flatten[Outer[List, List@First[list], Rest[list], 1], {{2}, {1, 4}}] paradox2[list_] := Partition[Riffle[list[[1]], #], 2] & /@ Drop[list, 1] RM[list_] := FoldList[Transpose[{First@li...
Post a Comment
Read more

front end - keyboard shortcut to invoke Insert new matrix

I frequently need to type in some matrices, and the menu command Insert > Table/Matrix > New... allows matrices with lines drawn between columns and rows, which is very helpful. I would like to make a keyboard shortcut for it, but cannot find the relevant frontend token command (4209405) for it. Since the FullForm[] and InputForm[] of matrices with lines drawn between rows and columns is the same as those without lines, it's hard to do this via 3rd party system-wide text expanders (e.g. autohotkey or atext on mac). How does one assign a keyboard shortcut for the menu item Insert > Table/Matrix > New... , preferably using only mathematica? Thanks! Answer In the MenuSetup.tr (for linux located in the $InstallationDirectory/SystemFiles/FrontEnd/TextResources/X/ directory), I changed the line MenuItem["&New...", "CreateGridBoxDialog"] to read MenuItem["&New...", "CreateGridBoxDialog", MenuKey["m", Modifiers-...
Post a Comment
Read more