- Aug 25, 2017
-
-
Dominic Kempf authored
-
Dominic Kempf authored
-
Dominic Kempf authored
-
Dominic Kempf authored
The introduction of FunctionView turned out to be a major problem with more complicated forms. The original idea was to preserver the structure of the finite element in a way, that loops over components of a mixed element are realized by actual loops (treating them with free indices and such). However, this causes quite some nightmares and was never implemented as generically as needed (I even doubt that is possible). However, there is another option, which is to unroll any such loops on a symbolic level. While this may sound like a bad idea at first there is some really positive aspects about it: * ListTensor and ComponentTensor nodes collapse completely (and would otherwise have a big nightmare potential) * Symbolic zeroes do not generate code - important in hyperbolic problems where the system matrices are quite sparse or for axiparallel grids, where geometric quantities have many zeroes. * The compiler would unroll these small loops anyway. * TSFC (and I guess also FFC) do it the same way. Implementing this required me to redo the form splitting algorithm. I rethought it and integrated it into the main ufl->loopy visitor.
-
- Aug 18, 2017
-
-
Marcel Koch authored
-
- Jul 31, 2017
-
-
Marcel Koch authored
-
- Jul 25, 2017
-
-
Marcel Koch authored
Computes tensor index from single index.
-
- Jul 24, 2017
-
-
Marcel Koch authored
-
- Jul 20, 2017
-
-
Marcel Koch authored
-
- Jul 14, 2017
-
-
Marcel Koch authored
-
Marcel Koch authored
-
- Jul 13, 2017
-
-
Marcel Koch authored
-
Marcel Koch authored
-
- Jul 12, 2017
-
-
Marcel Koch authored
-
Marcel Koch authored
Removes micro_to_macro_index from accum_insn and introduces a class BlockstructuredAccumulationSpace to specialize AccumulationSpace for blockstructured FEM
-
- Jul 07, 2017
-
-
Marcel Koch authored
-
- Jun 06, 2017
-
-
René Heß authored
-
- May 26, 2017
-
-
René Heß authored
-
- May 25, 2017
-
-
René Heß authored
-
- May 15, 2017
-
-
René Heß authored
-
- May 03, 2017
-
-
Dominic Kempf authored
-
Dominic Kempf authored
-
- Apr 27, 2017
-
-
Dominic Kempf authored
It was the wrong level of asbtraction!
-
- Apr 26, 2017
-
-
Dominic Kempf authored
-
- Apr 25, 2017
-
-
Dominic Kempf authored
-
- Apr 13, 2017
-
-
René Heß authored
-
- Mar 30, 2017
- Mar 29, 2017
- Mar 28, 2017
-
-
René Heß authored
-
Dominic Kempf authored
-
- Mar 27, 2017
-
-
René Heß authored
We avoid the bug by not applying the transformation and casting the result of jacobianInverseTransposed to a FieldMatrix. This way we do not exploit the diagonal structure of the Transformation. Should get fixed in the near future.
-
- Mar 24, 2017
-
-
René Heß authored
-
René Heß authored
Example: Before this commit the expression \sum_i a_i (\nabla v)_i + \sum_j b_j (\nabla v)_j was split into two accumulation terms: 1) a_k with corresponding test function (\nabla v)_k 2) b_l with corresponding test function (\nabla v)_l Now we split into: a_k + b_k with corresponding test function (\nabla v)_k This is possible since we have linearity in the test function. TODO: Jacobians are not yet working
-
René Heß authored
-
- Feb 03, 2017
-
-
Dominic Kempf authored
But only when they are actually used!
-
- Feb 02, 2017
-
-
Dominic Kempf authored
-
- Jan 25, 2017
-
-
Dominic Kempf authored
-
- Jan 24, 2017
-
-
Dominic Kempf authored
-