Comments (4)
I usually think of diagonal operators (mostly potentials in my case) as simply function multiplication:
g(x) = U(x) f(x)
where U(x)
is the diagonal operator. If you want it inside an integral, then you'd need the Dirac delta function, which of course complicates things if you want to index the operator itself, but is that something you need often?
In CoulombIntegrals.jl, I calculate the repulsion between two electrons, 1/abs(r_1 - r_2)
, by solving Poisson's problem over-and-over again, each time the electrons' wave functions have changed. The solution is a new potential in coefficient space, i.e. the distribution of e.g. electron 1 along r_1
forms a potential for electron 2. This is also a function, but only known at discrete points, e.g. a diagonal matrix acting on the coefficients of electron 2. Should this still be a QuasiDiagonal
? I don't think so, but at the same time, I usually want to merge all potentials into one before acting with the total operator to save FLOPS.
from continuumarrays.jl.
Right, that’s exactly how I’ve been thinking of them. There is no actual need for getindex except to make sure we are conceptualising them.
Perhaps getindex
is best thought of in terms of the delta function. that is K[s,t]
is equivalent to <δ_s, K*δ_t>
or in other words:
<δ_s, \int K(x,y) δ_t(y) dy > = < δ_s, K(x,t) > = K(s,t)
Perhaps the issue is the usage of QuasiDiagonal
/ Diagonal
, as it hints that it puts the quasi-vector on the diagonal, where really we have Inf * v
on the diagonal of QuasiDiagonal(v)
. Is there another name for this operation?
from continuumarrays.jl.
Maybe something with PointWise*
? In physics, we'd call this a local operator, whereas non-diagonal kernels are non-local; I prefer this terminology to point-wise.
from continuumarrays.jl.
Yes these are local operators, but so are derivatives. Pointwise
is too application specific. Perhaps follow ApproxFun and call it Multiplication
?
from continuumarrays.jl.
Related Issues (20)
- simplify macro does not permit templating? HOT 1
- How to handle Linear operators / functionals? HOT 4
- export diagonal?
- Continuum indices HOT 10
- Continous Linear Algebra HOT 1
- Integral of (basis) functions HOT 1
- Expansion short cut HOT 1
- Split out transform from factorize
- Derivative -> Diff? HOT 2
- Inner product between bases on different grids, basis transforms HOT 1
- Support syntax for kernels HOT 14
- Plot quasivector HOT 5
- [FEATURE]: Issue Template
- [FEATURE]: Pull Request Template
- Support (T/T) \ f for expansions HOT 1
- Product of QuasiDiagonals fails
- Support Plotting HOT 10
- Rename or don't export `grid` HOT 1
- Use bases for operations on AbstractQuasiVector
- Computing bounds on a function over an interval HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from continuumarrays.jl.