Obscure vector function translation from fortran

1 visualización (últimos 30 días)
David Winthrop
David Winthrop el 11 de Jul. de 2019
Respondida: Walter Roberson el 11 de Jul. de 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?
Thanks

Respuesta aceptada

Walter Roberson
Walter Roberson el 11 de Jul. de 2019
dot( v1(1:inc1:inc1*(n-1)),...
v2(1:inc2:inc2*(n-1)))
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
Off1=1
Off2=1
Result = 0
for k=1:n
Result = Result + v1(Off1) * v2(off2)
Off1 = Off1 + inc1
Off2 = Off2 + inc2
end
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.

Más respuestas (0)

Categorías

Más información sobre Fortran with MATLAB en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by