Wilkinson’s polynomial
If you change the coefficients of a polynomial a little bit, do you change the location of its zeros a little bit? In other words, do the roots of a polynomial depend continuously on its coefficients? You would think so, and you’d be right. Sorta. It’s easy to see that
Read moreHow to Build Scalable Web Apps with React JS
Unlock the best practices for building scalable React apps. Explore strategies for performance, maintainability, state management, code splitting, and real-world success stories. Continue reading How to Build Scalable Web Apps with React JS on SitePoint.
Read moreInterpolation instability
You would think that interpolating at more points would give you a more accurate approximation. There’s a famous example by Runge that proves this is not the case. If you interpolate the function 1/(1 + x²) over the interval [−5, 5], as you add more interpolation points the maximum interpolation
Read moreThe Ampere Porting Advisor Tutorial
This tutorial walks you through building and using the Ampere Porting Advisor and how to mitigate any issues. Continue reading The Ampere Porting Advisor Tutorial on SitePoint.
Read moreDrazin pseudoinverse
The most well-known generalization of the inverse of a matrix is the Moore-Penrose pseudoinverse. But there is another notion of generalized inverse, the Drazin pseudoinverse, for square matrices. If a matrix A has an inverse A−1 then it also has a Moore-Penrose pseudoinverse A+ and a Drazin pseudoinverse AD and A−1 =
Read moreEffective graph resistance
I’ve mentioned the Moore-Penrose pseudoinverse of a matrix a few times, most recently last week. This post will give an application of the pseudoinverse: computing effective graph resistance. Given a graph G, imagine replacing each edge with a one Ohm resistor. The effective resistance between two nodes in G is the electrical
Read moreMultiplying a matrix by its transpose
An earlier post claimed that there practical advantages to partitioning a matrix, thinking of the matrix as a matrix of matrices. This post will give an example. Let M be a square matrix and suppose we need to multiply M by its transpose MT. We can compute this product faster than
Read moreA bit-twiddling marvel
Pop quiz: what does the following code do? bool is_leap_year_fast(uint32_t y) { return ((y * 1073750999) & 3221352463) <= 126976; } It determines whether the year y is a leap year in the Gregorian calendar, of course. :) It’s very efficient, though I don’t image that would ever matter. But
Read moreCEO Update: Exploration and experimentation for bold evolution
As we envision what the ideal future version of Stack Overflow looks like, we’re committed to engaging with our community.
Read moreMatrices of Matrices
It’s often useful to partition a matrix, thinking of it as a matrix whose entries are matrices. For example, you could look at the matrix 6 × 6 as a 2 × 2 matrix whose entries are 3 × 3 matrices. M is not a diagonal matrix of real numbers,
Read more