Pardon my ignorance, is merging of loops/delayed evaluation available in Fortran for non-native arrays, when the binary operator is overloaded for a user defined type? If yes, then we have a clear edge over C++. If no, then we should look out for a mechanism to incorporate.
Interesting to note! Long ago, when C++11 was under preparation and was known as C++0x and I had not learned Fortran, I had a small communication with Stroustrup on whether he has specific plans for supporting numerics in the upcoming standard. He was quite positive about it. But, it later turned out that most in their committee were not considerate. There argument, I read somewhere, was that numerics people have a dedicated language: Fortran. Core C++ should be kept application agnostic. The trend, it appears, has buckled. But, I no longer feel attracted towards C++.