Why did not numpy copy the J rank concept?

11•jrank•7h ago

Recently there was a post [1] about numpy being not easy to work with when using arrays of shape greater than 2.

  One of the problems mentioned in [1] is that  you can not use python loops because they are slow.  In J you can solve for example 100 equations using the rank concept, this is a simple  example:

   a=: 2 2 $ 1 1 1 _1
   b=: 10 2 $ ? 20 # 10
   solutions =: b %. "(1 _) a

That code solve the systems a * v_i = b_i for ten random vectors. I think a similar concept could be developed in numpy. The syntax "(1 _)" indicates to take the rows from the left operator and all (_ is infinite) of a apply solve (that is %. in J). In this case the system is x+y=y0, x-y=y1.

So I would suggest somthing like numpy.linalg.solve(a,b,rank=(1,inf))

[1] https://news.ycombinator.com/item?id=43996431

Comments

tester89•5h ago

The question you're asking seems interesting, but I don't understand J code so I don't know what you're talking about. Expanding the explanation would be helpful!

jrank•1h ago

If you have a a Numpy function that takes for example two arguments, my proposal is to add a optional argument that allows that function to be applied to cells of dimensions i and j by function_name(a,b,range=(i,j)) so that the function is applied to subarrays of dimension i of a and subarrays of dimension j of b to create a new array. The broadcast operation and the axis arguments are not a general solution. I J you have such mechanism as the example shows.

The 1 cell of a matrix are the rows, etc.

Pompidou•3h ago

Maybe the internal broadcasting mecanism in numpy don't allow this nativelly ?

jrank•7m ago

The post I referenced before shows that broadcast is not a general solution.

Amiga Graphics Archive adds a big selection of color cycling images

Show HN: An online 2d shooter game in development

Typing Speed Is the Hidden Signal of a 10x Developer

Stop deploying web application firewalls (2023)

Video Cable Becomes Transmitter with Tempest-LoRa

The Tandy Corporation, Part 2

Cursor's 500 requests => "unlimited" => 225 requests

Regulators warned Air India Express about forging records

Dreadnought Hoax

Ask HN: How many communities HN it devs in C language?

Broken AI Discourse with Steve Klabnik [video]

Show HN: Liryo – Book journal app that tracks stats, goals and spending

Loggers close in on uncontacted people in Peruvian Amazon

Speeding up PostgreSQL dump/restore snapshots

The New York Times News Quiz, July 4, 2025

The Calculator-on-a-Chip (2015)

Long-lost Chinese typewriter prototype from the 1940s that changed computing

Thesis Drift

Show HN: I designed a math boardgame. HTML mock-up is now a playable game

EU says it will continue rolling out AI legislation on schedule

Senate Tax Bill Expands QSBS Benefits for Startups

Git experts should try Jujutsu

Cursor Pricing Changes

ApplePay vs. Alternative Payment Services

Is There Is Any PNG and SVG Icons Website That Work Without JavaScript?

Adding a new instruction to RISC-V back end in LLVM

Cops in [Spain] think everyone using a Google Pixel must be a drug dealer

Goodbye to All That – My Resignation from the FBI

Show HN: I Made a Hot or Not Benchmark for AI Design

The curious rise of giant tablets on wheels