Is Frobenius norm constant?











up vote
0
down vote

favorite












I came across an article that makes me doubt what I have learnt so far.



For example, let us decompose matrix A into L and U (LU factorization)
$A=left(
begin{array}{cccc}
5 & 4 & 1 & 1 \
4 & 5 & 1 & 1 \
1&1&4&2 \
1&1&2&4
end{array}right)$



The Frobenius norm should be $||A-LU||$. Now, the authors gave different answers when they implement this on MATLAB 2010 and 2012 and also different answers on AMD and Intel processor. Why should we have different answers for a consistent norm?










share|cite|improve this question






















  • I think something is not quite right with your question. The Frobenius norm of A is the sum of squared entries of $A$. Is $| cdot |$ denoting the Frobenius norm in your question? We certainly wouldn't have $|A|=|A-LU|$. As for why the norm may be different on different machines or versions of matlab, it may have to do with the precision the computation was run at. If the difference is very small this is likely the case. Otherwise something is wrong.
    – Eric
    2 days ago












  • I do not still get it. When L and U are multiplied and then subtracted from A, the result is A-LU. When we take the Frobenius norm of this result, should we have different answers on different machine or versions of MATLAb?
    – dupsy
    yesterday










  • Your question is somewhat implicitly assuming the existence of an LU decomposition. That means that theoretically A=LU, so and we have $|A-LU|_F= |0|_F=0$. However, you might not get exactly zero in a computer. This is because the entries of $L$ and $U$ may be irrational numbers. A computer cannot store an irrational number perfectly, meaning the $L$ and $U$ stored in your machine will be very slightly different than the correct ones. This leads to small amounts of round off error when computing $A-LU$. So if all answers are very close, the error is probably just round off error.
    – Eric
    yesterday















up vote
0
down vote

favorite












I came across an article that makes me doubt what I have learnt so far.



For example, let us decompose matrix A into L and U (LU factorization)
$A=left(
begin{array}{cccc}
5 & 4 & 1 & 1 \
4 & 5 & 1 & 1 \
1&1&4&2 \
1&1&2&4
end{array}right)$



The Frobenius norm should be $||A-LU||$. Now, the authors gave different answers when they implement this on MATLAB 2010 and 2012 and also different answers on AMD and Intel processor. Why should we have different answers for a consistent norm?










share|cite|improve this question






















  • I think something is not quite right with your question. The Frobenius norm of A is the sum of squared entries of $A$. Is $| cdot |$ denoting the Frobenius norm in your question? We certainly wouldn't have $|A|=|A-LU|$. As for why the norm may be different on different machines or versions of matlab, it may have to do with the precision the computation was run at. If the difference is very small this is likely the case. Otherwise something is wrong.
    – Eric
    2 days ago












  • I do not still get it. When L and U are multiplied and then subtracted from A, the result is A-LU. When we take the Frobenius norm of this result, should we have different answers on different machine or versions of MATLAb?
    – dupsy
    yesterday










  • Your question is somewhat implicitly assuming the existence of an LU decomposition. That means that theoretically A=LU, so and we have $|A-LU|_F= |0|_F=0$. However, you might not get exactly zero in a computer. This is because the entries of $L$ and $U$ may be irrational numbers. A computer cannot store an irrational number perfectly, meaning the $L$ and $U$ stored in your machine will be very slightly different than the correct ones. This leads to small amounts of round off error when computing $A-LU$. So if all answers are very close, the error is probably just round off error.
    – Eric
    yesterday













up vote
0
down vote

favorite









up vote
0
down vote

favorite











I came across an article that makes me doubt what I have learnt so far.



For example, let us decompose matrix A into L and U (LU factorization)
$A=left(
begin{array}{cccc}
5 & 4 & 1 & 1 \
4 & 5 & 1 & 1 \
1&1&4&2 \
1&1&2&4
end{array}right)$



The Frobenius norm should be $||A-LU||$. Now, the authors gave different answers when they implement this on MATLAB 2010 and 2012 and also different answers on AMD and Intel processor. Why should we have different answers for a consistent norm?










share|cite|improve this question













I came across an article that makes me doubt what I have learnt so far.



For example, let us decompose matrix A into L and U (LU factorization)
$A=left(
begin{array}{cccc}
5 & 4 & 1 & 1 \
4 & 5 & 1 & 1 \
1&1&4&2 \
1&1&2&4
end{array}right)$



The Frobenius norm should be $||A-LU||$. Now, the authors gave different answers when they implement this on MATLAB 2010 and 2012 and also different answers on AMD and Intel processor. Why should we have different answers for a consistent norm?







matrices norm matrix-decomposition






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked 2 days ago









dupsy

475




475












  • I think something is not quite right with your question. The Frobenius norm of A is the sum of squared entries of $A$. Is $| cdot |$ denoting the Frobenius norm in your question? We certainly wouldn't have $|A|=|A-LU|$. As for why the norm may be different on different machines or versions of matlab, it may have to do with the precision the computation was run at. If the difference is very small this is likely the case. Otherwise something is wrong.
    – Eric
    2 days ago












  • I do not still get it. When L and U are multiplied and then subtracted from A, the result is A-LU. When we take the Frobenius norm of this result, should we have different answers on different machine or versions of MATLAb?
    – dupsy
    yesterday










  • Your question is somewhat implicitly assuming the existence of an LU decomposition. That means that theoretically A=LU, so and we have $|A-LU|_F= |0|_F=0$. However, you might not get exactly zero in a computer. This is because the entries of $L$ and $U$ may be irrational numbers. A computer cannot store an irrational number perfectly, meaning the $L$ and $U$ stored in your machine will be very slightly different than the correct ones. This leads to small amounts of round off error when computing $A-LU$. So if all answers are very close, the error is probably just round off error.
    – Eric
    yesterday


















  • I think something is not quite right with your question. The Frobenius norm of A is the sum of squared entries of $A$. Is $| cdot |$ denoting the Frobenius norm in your question? We certainly wouldn't have $|A|=|A-LU|$. As for why the norm may be different on different machines or versions of matlab, it may have to do with the precision the computation was run at. If the difference is very small this is likely the case. Otherwise something is wrong.
    – Eric
    2 days ago












  • I do not still get it. When L and U are multiplied and then subtracted from A, the result is A-LU. When we take the Frobenius norm of this result, should we have different answers on different machine or versions of MATLAb?
    – dupsy
    yesterday










  • Your question is somewhat implicitly assuming the existence of an LU decomposition. That means that theoretically A=LU, so and we have $|A-LU|_F= |0|_F=0$. However, you might not get exactly zero in a computer. This is because the entries of $L$ and $U$ may be irrational numbers. A computer cannot store an irrational number perfectly, meaning the $L$ and $U$ stored in your machine will be very slightly different than the correct ones. This leads to small amounts of round off error when computing $A-LU$. So if all answers are very close, the error is probably just round off error.
    – Eric
    yesterday
















I think something is not quite right with your question. The Frobenius norm of A is the sum of squared entries of $A$. Is $| cdot |$ denoting the Frobenius norm in your question? We certainly wouldn't have $|A|=|A-LU|$. As for why the norm may be different on different machines or versions of matlab, it may have to do with the precision the computation was run at. If the difference is very small this is likely the case. Otherwise something is wrong.
– Eric
2 days ago






I think something is not quite right with your question. The Frobenius norm of A is the sum of squared entries of $A$. Is $| cdot |$ denoting the Frobenius norm in your question? We certainly wouldn't have $|A|=|A-LU|$. As for why the norm may be different on different machines or versions of matlab, it may have to do with the precision the computation was run at. If the difference is very small this is likely the case. Otherwise something is wrong.
– Eric
2 days ago














I do not still get it. When L and U are multiplied and then subtracted from A, the result is A-LU. When we take the Frobenius norm of this result, should we have different answers on different machine or versions of MATLAb?
– dupsy
yesterday




I do not still get it. When L and U are multiplied and then subtracted from A, the result is A-LU. When we take the Frobenius norm of this result, should we have different answers on different machine or versions of MATLAb?
– dupsy
yesterday












Your question is somewhat implicitly assuming the existence of an LU decomposition. That means that theoretically A=LU, so and we have $|A-LU|_F= |0|_F=0$. However, you might not get exactly zero in a computer. This is because the entries of $L$ and $U$ may be irrational numbers. A computer cannot store an irrational number perfectly, meaning the $L$ and $U$ stored in your machine will be very slightly different than the correct ones. This leads to small amounts of round off error when computing $A-LU$. So if all answers are very close, the error is probably just round off error.
– Eric
yesterday




Your question is somewhat implicitly assuming the existence of an LU decomposition. That means that theoretically A=LU, so and we have $|A-LU|_F= |0|_F=0$. However, you might not get exactly zero in a computer. This is because the entries of $L$ and $U$ may be irrational numbers. A computer cannot store an irrational number perfectly, meaning the $L$ and $U$ stored in your machine will be very slightly different than the correct ones. This leads to small amounts of round off error when computing $A-LU$. So if all answers are very close, the error is probably just round off error.
– Eric
yesterday















active

oldest

votes











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});

function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














 

draft saved


draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3005570%2fis-frobenius-norm-constant%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown






























active

oldest

votes













active

oldest

votes









active

oldest

votes






active

oldest

votes
















 

draft saved


draft discarded



















































 


draft saved


draft discarded














StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3005570%2fis-frobenius-norm-constant%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

MongoDB - Not Authorized To Execute Command

How to fix TextFormField cause rebuild widget in Flutter

in spring boot 2.1 many test slices are not allowed anymore due to multiple @BootstrapWith