Any possible optimization?

[Daniel Arteaga]

Hi all

The following function is called many times in my code (always from the 
same point), and it is actually the bottleneck of the computation:

function [s,t] = intersectionPoint(a,b,c,d)

        % Solve the system of equations
        % (1 - s) ax + s bx = (1 - t) cx + t dx
        % (1 - s) ay + s by = (1 - t) cy + t dy

        A = [ b(1) - a(1), c(1) - d(1); b(2) - a(2), c(2) - d(2)];
        B = [ c(1) - a(1); c(2) - a(2)];
        
        v = A\B;
        s = v(1);
        t = v(2);

endfunction


Is there any easy way to get a speed increase? Is it worth transforming 
it into an .oct file? Or maybe including the expressions inline to avoid 
the function call overhead?

Any suggestion is appreciated.

Thank you very much,

Daniel Arteaga