How do I properly subtract two column vectors of class 'Double'

Hi!

Im having trouble subtracting two column vectors of class = 'Double'

I have the equation:

v = A*x - l

A*x is very similar to L (both being 595x1), so i expect very small differences in the output variable v.

However, when running this in my function the output results are not as expected.

I tried isolating A*x and L and get them out of my function alone. I then did the subtraction manually in the command window afterwards and it computed the result i was expecting.

What am I missing?

[Torsten]

We can't tell since we don't know your code. Please include the relevant part so that we can check. Maybe you used "l" instead of "L" as in your question ?

[Max Alger-Meyer]

@Mohammed Sahan As Tosten mentioned, this code is trivial, and correct. So if the output isn't as expected, that likely means that either you are mistaken and the values of A*x and L are not similar, or at some point between when that operation takes place and when you look at the output, 'v' is getting changed. Be sure to put a breakpoint on the line that the subtraction takes place, and check the values of A*x and L in your command window. Then, put another breakpoint on the next line and check the values of 'v'. If the values of 'v' are as expected, then 'v' is likely getting modified later on in your function or script. You'll have to post the function and inputs if you need anything more than that for us to help you.

[Mohammed Sahan]

Ill copy from v = A*x - l to the end.

All variables below are mostly constants or calculated based on constants. Only v, sgg and s0 are of interest.

v = A*x-1; % Residuals

% 9
g_first = g0 - gamma.*z0; % Reference height correction

% 10
dg_transfer = (h_factory + h_su- h_transfer)*gamma; % Transfer height correction
g_transfer = g_first + dg_transfer; 

% 11
g = g_transfer + dg_set + dg_otl + dg_atl + dg_pt; % Corrections for SET, OTL, ATL, PT

% 12
sgg = sqrt((v'*v)/(n-e)) * sqrt(Q(3,3))

%13
s0 = sqrt((v'*v)/(n-e));

As seen above, i cannot find any place where i modify my 'v' after computing it.

As an example of what i expect ill use the two first rows from A*x and l:

Ax = [0.007912389029202;
      0.008228885538678]

l = [0.007912389934100;
     0.008228885531464]

The values are very similar but there should be some difference at around the 9th decimal place.

When computing v in the function i get the following results for the two first rows:

v = [-0.992087610970798;
     -0.991771114461323]

This is clearly not correct, when im expecting something like:

% row1 = 0.007912389029202 - 0.007912389934100 = -9.04898e-10
% row2 = 0.008228885538678 - 0.008228885531464 = 7.2140002e-12

I end up getting the output i expect when calling out A*x and l from the function beforehand and then computing 'v', but not when i calculate it inside the function.

Answer 1

[Torsten]

You calculate

v = A*x - 1

in words: A*x - one

instead of

v = A*x - l

in words: A*x - el

[Mohammed Sahan]

Well, thats awkward :) Sorry for wasting everyones time and thank you for catching that!