Just to follow up on this thread in case anyone else is interested, I found
that the Einsum.jl package is quite close to what I need for my FEM
application. It is lacking a couple of features, so I added them to the
package and have submitted a pull request. Hopefully the author of
https://github.com/ahwillia/Einsum.jl
or
https://github.com/Jutho/TensorOperations.jl
might also be of interest
Op vrijdag 20 mei 2016 17:17:08 UTC+2 schreef Matt Bauman:
>
> I believe something just like that is already implemented in
> Devectorize.jl's @devec macro. Try: `@devec y[:] = y + a
I believe something just like that is already implemented in
Devectorize.jl's @devec macro. Try: `@devec y[:] = y + a .* x[ind]`.
https://github.com/lindahua/Devectorize.jl
On Friday, May 20, 2016 at 11:04:34 AM UTC-4, vav...@uwaterloo.ca wrote:
>
> I'm not sure that fixed-size array solves
I'm not sure that fixed-size array solves the problem. For one thing, the
size of the small array may not be known at compile time. For another, I
don't want the small vectors necessarily to be immutable.
I decided that what would really solve this problem is indicial notation.
This can be
If you have a lot of small arrays, you might consider FixedSizeArrays, which
may address some of these points.
Best,
--Tim
On Thursday, May 19, 2016 07:20:22 PM vava...@uwaterloo.ca wrote:
> Dear Yichao,
>
> Thanks very much for the prompt response. This question arises regarding a
> code for
On Thursday, May 19, 2016 at 10:20:22 PM UTC-4, vav...@uwaterloo.ca wrote:
>
> Dear Yichao,
>
> Thanks very much for the prompt response. This question arises regarding
> a code for finite element stiffness matrix assembly. This computation
> involves an outer loop over elements (possibly
Dear Yichao,
Thanks very much for the prompt response. This question arises regarding a
code for finite element stiffness matrix assembly. This computation
involves an outer loop over elements (possibly millions of them). Inside
this loop is a sequence of operations on 'small' vectors and
