Comments (6)
This is an interesting case (which I don't remember seeing the likes of). By the time factorization is reached, the matrix is reduced to something like this:
{{-8*eps/(x-1)+8*eps/x,0,0,0,0},{0,-4*eps/(x-1)+6*eps/x,0,0,0},{0,-1/36*(240*eps^3-236*eps^2+56*eps+1)/(eps^2*(x-1))+1/36*(112*eps^3-172*eps^2+88*eps-15)/(eps^2*x),1/4*eps/(x-1)-1/4*eps/x,-25/72*eps/(x-1)+289/72*eps/x,0},{0,1/2*(16*eps^3-20*eps^2+8*eps-1)/(eps^2*(x-1))-1/2*(16*eps^3-20*eps^2+8*eps-1)/(eps^2*x),9/2*eps/(x-1)-9/2*eps/x,-25/4*eps/(x-1)+33/4*eps/x,0},{(4*eps-1)/(2*x-1)-(4*eps-1)/x,-8*(4*eps-1)/(x-1)-4*(12*eps^2-4*eps+1)/(eps*(2*x-1))+4*(8*eps^2-6*eps+1)/(eps*x),-72*eps^3/((4*eps^2-4*eps+1)*(x-1))+18*(2*eps^3+eps^2)/((4*eps^2-4*eps+1)*(2*x-1)),100*eps^3/((4*eps^2-4*eps+1)*(x-1))-17*(2*eps^3+eps^2)/((4*eps^2-4*eps+1)*(2*x-1)),-4*eps/(2*x-1)+4*eps/x}}
This indeed looks like a fully normalized Fuchsian matrix, so no error so far. Next, factorize
tries to find an x-independent transformation to turn this into epsilon-form, but finds none. I've looked into the generated equations a bit, and the conclusion seems correct: there is no constant transformation that can turn this matrix into epsilon-form. There's not much Fuchsia can do about that.
Now, the interesting question is: can there be a non-constant transformation that reduces this matrix? This unfortunately I don't know.
from fuchsia.
Thanks so much for looking into such a case!
I am going to do some try on the normalized Fuchisan result,
and maybe try the solution of DE directly from the pre-factorize matrix ...
Anyway, thanks so much again!
from fuchsia.
@F-Feng, could you, please, shortly describe in a context of which problem have you obtained this matrix?
from fuchsia.
@gituliar, thanks, and sure.
The matrix shows up when I try to solve some master integrals which comes from this paper: https://arxiv.org/abs/1708.01129, the master integrals are listed in Eq.(44) of that paper, and when I try to solve those master integrals by differential equation, I get the matrix I sent to you.
from fuchsia.
@F-Feng, thank you for the details. As @magv mentioned, your example is quite interesting since we have not seen Fuchsia ever fail at the last step, though that is possible in principle (for more details see Section 6 in http://arxiv.org/abs/1411.0911v1). For the moment, the only explanation I have for such behavior, is that your 5x5 matrix is part of a non-homogeneous system, since in the paper you mentioned authors deal with 13 masters. Hence, I would suggest you to try a complete 13x13 matrix. I think that in this case chances that Fuchsia will find a proper transformation are much higher.
from fuchsia.
@gituliar Thanks very much for the suggestion, I will try the complete master integrals, and if I succeed, I will post the detailed information. Thanks again!
from fuchsia.
Related Issues (17)
- Project's name HOT 5
- Performance issue with `fuchsify` and `git_410.mtx` HOT 7
- [Algorithm 2] Avoid tremendous polynoms in eps HOT 4
- CLI issues HOT 5
- Infinite loop in the normalizer HOT 7
- Barkatou examples from ISSAC'10 HOT 5
- Permutation to the block-triangular form needed
- Release #1 HOT 3
- Bug with `solve` in `singularities` HOT 13
- 'fuchsia' binary installed by pip doesn't work HOT 3
- Multivarate dependence HOT 4
- NameError: Unknown variable: 'sqrt' HOT 4
- Some unit tests fail HOT 8
- Different output when using Maple and Maxima HOT 2
- Differential equations with free variables HOT 4
- Transform out spurious poles without solving the polynomials
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 fuchsia.