Obscure vector function translation from fortran
2 views (last 30 days)
David Winthrop on 11 Jul 2019
I have been tasked with translating some code into MATLAB. I am puzzled at this function. It is from a FEA software called ANSYS, and the only documentation is shown in the image below. I am having trouble figuring out what exactly this function does. Anyone ever heard of this before? It looks like the input argument inc1 just tells the function the number of rows in the vector, and similarly for inc2. Then the input argument n tells the number of columns (are we talking about vectors or arrays??) What is going on here?
Walter Roberson on 11 Jul 2019
This matlab rendition creates temporary index vectors and creates temporary vectors of extracted elements. That works OK, but is not as efficient at the low level as compiled code along the lines of
Result = 0
Result = Result + v1(Off1) * v2(off2)
Off1 = Off1 + inc1
Off2 = Off2 + inc2
Now let us consider
A * B %matrix multiply
That involves a whole series of
dot( A(:, k), B(k, :) )
Which is vidot(&A(1,k),1,&B(k,1),size(B,1),size(A,1))
Where here & is intended to indicate "address of"
You could implement a matrix multiply with a lot of temporary index vectors and temporary extraction of vectors, but that involves a lot of temporary operations and storage management that you can see are not needed if you know the distance between vector elements.