← Back to branch summary

~ubuntu-branches/ubuntu/utopic/ffc/utopic

« back to all changes in this revision

Viewing changes to test/regression/references/FacetIntegrals.h

  • Committer: Bazaar Package Importer
  • Author(s): Johannes Ring
  • Date: 2011-05-18 11:22:37 UTC
  • mfrom: (1.1.8 upstream)
  • Revision ID: james.westby@ubuntu.com-20110518112237-22o6lwrno3uf0jyz
Tags: 0.9.10-1
http://bugs.debian.org/616807
* New upstream release. This release introduces some new features and
  some behind-the-scenes improvements. Special quadrature schemes for
  simplices have been added, which will lead to improved performance
  when using the quadrature representation for low-order forms.
* Move from python-central to dh_python2 (closes: #616807).
  - Remove python-central from Build-Depends.
  - Bump minimum required python-all package version to 2.6.6-3~.
  - Remove XB-Python-Version line.
  - Bump minimum required cdbs version to 0.4.90~.
  - Remove DEB_PYTHON_SYSTEM=pycentral from debian/rules.
  - Replace XS-Python-Version with X-Python-Version.
* debian/rules:
  - Avoid compressing files with .ufl extension.
  - Use DEB_COMPRESS_EXCLUDE_ALL instead of deprecated
    DEB_COMPRESS_EXCLUDE.
  - Include cdbs utils.mk rule for automated copyright checks.
* Remove old custom cdbs rules and licensecheck script for copyright
  check.
* Bump Standards-Version to 3.9.2 (no changes needed).
* debian/copyright:
  - Update for upstream license change to LGPLv3.
  - Switch to DEP-5 format.

Show diffs side-by-side

added added

removed removed

Lines of Context:
test/regression/references/FacetIntegrals.h
1
1
// This code conforms with the UFC specification version 2.0.0
2
 
// and was automatically generated by FFC version 0.9.9.
 
2
// and was automatically generated by FFC version 0.9.10.
3
3
// 
4
4
// This code was generated with the following parameters:
5
5
// 
20
20
//   quadrature_rule:                'auto'
21
21
//   representation:                 'auto'
22
22
//   split:                          False
 
23
//   swig_binary:                    'swig'
 
24
//   swig_path:                      ''
23
25
 
24
26
#ifndef __FACETINTEGRALS_H
25
27
#define __FACETINTEGRALS_H
50
52
  /// Return a string identifying the finite element
51
53
  virtual const char* signature() const
52
54
  {
53
 
    return "FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 0)";
 
55
    return "FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 0, None)";
54
56
  }
55
57
 
56
58
  /// Return the cell shape
428
430
  /// Map coordinate xhat from reference cell to coordinate x in cell
429
431
  virtual void map_from_reference_cell(double* x,
430
432
                                       const double* xhat,
431
 
                                       const ufc::cell& c)
 
433
                                       const ufc::cell& c) const
432
434
  {
433
435
    std::cerr << "*** FFC warning: " << "map_from_reference_cell not yet implemented (introduced in UFC 2.0)." << std::endl;
434
436
  }
436
438
  /// Map from coordinate x in cell to coordinate xhat in reference cell
437
439
  virtual void map_to_reference_cell(double* xhat,
438
440
                                     const double* x,
439
 
                                     const ufc::cell& c)
 
441
                                     const ufc::cell& c) const
440
442
  {
441
443
    std::cerr << "*** FFC warning: " << "map_to_reference_cell not yet implemented (introduced in UFC 2.0)." << std::endl;
442
444
  }
482
484
  /// Return a string identifying the finite element
483
485
  virtual const char* signature() const
484
486
  {
485
 
    return "VectorElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 0, 2)";
 
487
    return "VectorElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 0, 2, None)";
486
488
  }
487
489
 
488
490
  /// Return the cell shape
1090
1092
  /// Map coordinate xhat from reference cell to coordinate x in cell
1091
1093
  virtual void map_from_reference_cell(double* x,
1092
1094
                                       const double* xhat,
1093
 
                                       const ufc::cell& c)
 
1095
                                       const ufc::cell& c) const
1094
1096
  {
1095
1097
    std::cerr << "*** FFC warning: " << "map_from_reference_cell not yet implemented (introduced in UFC 2.0)." << std::endl;
1096
1098
  }
1098
1100
  /// Map from coordinate x in cell to coordinate xhat in reference cell
1099
1101
  virtual void map_to_reference_cell(double* xhat,
1100
1102
                                     const double* x,
1101
 
                                     const ufc::cell& c)
 
1103
                                     const ufc::cell& c) const
1102
1104
  {
1103
1105
    std::cerr << "*** FFC warning: " << "map_to_reference_cell not yet implemented (introduced in UFC 2.0)." << std::endl;
1104
1106
  }
1158
1160
  /// Return a string identifying the finite element
1159
1161
  virtual const char* signature() const
1160
1162
  {
1161
 
    return "FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1)";
 
1163
    return "FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None)";
1162
1164
  }
1163
1165
 
1164
1166
  /// Return the cell shape
2061
2063
  /// Map coordinate xhat from reference cell to coordinate x in cell
2062
2064
  virtual void map_from_reference_cell(double* x,
2063
2065
                                       const double* xhat,
2064
 
                                       const ufc::cell& c)
 
2066
                                       const ufc::cell& c) const
2065
2067
  {
2066
2068
    std::cerr << "*** FFC warning: " << "map_from_reference_cell not yet implemented (introduced in UFC 2.0)." << std::endl;
2067
2069
  }
2069
2071
  /// Map from coordinate x in cell to coordinate xhat in reference cell
2070
2072
  virtual void map_to_reference_cell(double* xhat,
2071
2073
                                     const double* x,
2072
 
                                     const ufc::cell& c)
 
2074
                                     const ufc::cell& c) const
2073
2075
  {
2074
2076
    std::cerr << "*** FFC warning: " << "map_to_reference_cell not yet implemented (introduced in UFC 2.0)." << std::endl;
2075
2077
  }
2119
2121
  /// Return a string identifying the dofmap
2120
2122
  virtual const char* signature() const
2121
2123
  {
2122
 
    return "FFC dofmap for FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 0)";
 
2124
    return "FFC dofmap for FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 0, None)";
2123
2125
  }
2124
2126
 
2125
2127
  /// Return true iff mesh entities of topological dimension d are needed
2351
2353
  /// Return a string identifying the dofmap
2352
2354
  virtual const char* signature() const
2353
2355
  {
2354
 
    return "FFC dofmap for VectorElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 0, 2)";
 
2356
    return "FFC dofmap for VectorElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 0, 2, None)";
2355
2357
  }
2356
2358
 
2357
2359
  /// Return true iff mesh entities of topological dimension d are needed
2604
2606
  /// Return a string identifying the dofmap
2605
2607
  virtual const char* signature() const
2606
2608
  {
2607
 
    return "FFC dofmap for FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1)";
 
2609
    return "FFC dofmap for FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None)";
2608
2610
  }
2609
2611
 
2610
2612
  /// Return true iff mesh entities of topological dimension d are needed
3063
3065
        A[7] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3064
3066
        A[8] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3065
3067
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1;
3066
 
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.33333333*G8_;
3067
 
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.16666667*G8_;
 
3068
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
 
3069
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3068
3070
        A[12] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1;
3069
3071
        A[13] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3070
3072
        A[14] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3071
3073
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1;
3072
 
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.16666667*G8_;
3073
 
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3074
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
 
3075
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3074
3076
        A[18] = 0.00000000;
3075
3077
        A[19] = -0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3076
3078
        A[20] = -0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3078
3080
        A[22] = -0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3079
3081
        A[23] = -0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3080
3082
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1;
3081
 
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
3082
 
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3083
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.33333333*G8_;
 
3084
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.16666667*G8_;
3083
3085
        A[27] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1;
3084
3086
        A[28] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3085
3087
        A[29] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3086
3088
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1;
3087
 
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
3088
 
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3089
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.16666667*G8_;
 
3090
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3089
3091
        A[33] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1;
3090
3092
        A[34] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3091
3093
        A[35] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3102
3104
        A[6] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1;
3103
3105
        A[7] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3104
3106
        A[8] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3105
 
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.33333333*G8_;
 
3107
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3106
3108
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1;
3107
 
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.16666667*G8_;
 
3109
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3108
3110
        A[12] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1;
3109
3111
        A[13] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3110
3112
        A[14] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3111
 
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3113
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3112
3114
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1;
3113
 
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3115
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3114
3116
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1;
3115
 
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3116
 
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3117
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.33333333*G8_;
 
3118
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3117
3119
        A[21] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3118
3120
        A[22] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1;
3119
3121
        A[23] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3124
3126
        A[28] = 0.00000000;
3125
3127
        A[29] = 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3126
3128
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1;
3127
 
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
3128
 
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3129
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.16666667*G8_;
 
3130
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3129
3131
        A[33] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3130
3132
        A[34] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1;
3131
3133
        A[35] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3142
3144
        A[6] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1;
3143
3145
        A[7] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3144
3146
        A[8] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3145
 
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.33333333*G8_;
3146
 
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.16666667*G8_;
 
3147
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
 
3148
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3147
3149
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1;
3148
3150
        A[12] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1;
3149
3151
        A[13] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3150
3152
        A[14] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3151
 
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.16666667*G8_;
3152
 
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3153
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
 
3154
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3153
3155
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1;
3154
3156
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1;
3155
 
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3156
 
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3157
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.33333333*G8_;
 
3158
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3157
3159
        A[21] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3158
3160
        A[22] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3159
3161
        A[23] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1;
3160
3162
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1;
3161
 
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
3162
 
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3163
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.16666667*G8_;
 
3164
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.33333333*G8_;
3163
3165
        A[27] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3164
3166
        A[28] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3165
3167
        A[29] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1;
3185
3187
        A[1] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1;
3186
3188
        A[2] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3187
3189
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1;
3188
 
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.33333333*G8_;
3189
 
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3190
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
 
3191
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3190
3192
        A[6] = 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3191
3193
        A[7] = 0.00000000;
3192
3194
        A[8] = 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3197
3199
        A[13] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1;
3198
3200
        A[14] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3199
3201
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1;
3200
 
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.16666667*G8_;
3201
 
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3202
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
 
3203
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3202
3204
        A[18] = -0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3203
3205
        A[19] = 0.00000000;
3204
3206
        A[20] = -0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3205
3207
        A[21] = 0.00000000;
3206
3208
        A[22] = -0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3207
3209
        A[23] = -0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3208
 
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3210
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.33333333*G8_;
3209
3211
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1;
3210
 
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3212
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.16666667*G8_;
3211
3213
        A[27] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1;
3212
3214
        A[28] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3213
3215
        A[29] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3214
 
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3216
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3215
3217
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1;
3216
 
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3218
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3217
3219
        A[33] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1;
3218
3220
        A[34] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3219
3221
        A[35] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3224
3226
          A[0] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3225
3227
        A[1] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1;
3226
3228
        A[2] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3227
 
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3229
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3228
3230
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1;
3229
 
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3231
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3230
3232
        A[6] = 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3231
3233
        A[7] = 0.00000000;
3232
3234
        A[8] = 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3236
3238
        A[12] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3237
3239
        A[13] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1;
3238
3240
        A[14] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3239
 
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3241
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3240
3242
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1;
3241
 
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.33333333*G8_;
3242
 
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3243
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
 
3244
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3243
3245
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1;
3244
 
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3246
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3245
3247
        A[21] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3246
3248
        A[22] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1;
3247
3249
        A[23] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3251
3253
        A[27] = 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3252
3254
        A[28] = 0.00000000;
3253
3255
        A[29] = 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3254
 
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3256
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3255
3257
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1;
3256
 
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3258
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3257
3259
        A[33] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3258
3260
        A[34] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1;
3259
3261
        A[35] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3264
3266
          A[0] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3265
3267
        A[1] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1;
3266
3268
        A[2] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3267
 
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.33333333*G8_;
3268
 
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3269
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
 
3270
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3269
3271
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1;
3270
3272
        A[6] = 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3271
3273
        A[7] = 0.00000000;
3276
3278
        A[12] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3277
3279
        A[13] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1;
3278
3280
        A[14] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1 + 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3279
 
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.16666667*G8_;
3280
 
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3281
        A[15] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
 
3282
        A[16] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3281
3283
        A[17] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1;
3282
 
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3284
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3283
3285
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1;
3284
 
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3286
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3285
3287
        A[21] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3286
3288
        A[22] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3287
3289
        A[23] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1;
3288
 
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3290
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.16666667*G8_;
3289
3291
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1;
3290
 
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3292
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.33333333*G8_;
3291
3293
        A[27] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3292
3294
        A[28] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3293
3295
        A[29] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1;
3313
3315
        A[1] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3314
3316
        A[2] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1;
3315
3317
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1;
3316
 
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.33333333*G8_;
3317
 
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3318
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
 
3319
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3318
3320
        A[6] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3319
3321
        A[7] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3320
3322
        A[8] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1;
3321
3323
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1;
3322
 
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.16666667*G8_;
3323
 
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.33333333*G8_;
 
3324
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
 
3325
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3324
3326
        A[12] = 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3325
3327
        A[13] = 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3326
3328
        A[14] = 0.00000000;
3333
3335
        A[21] = 0.00000000;
3334
3336
        A[22] = -0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3335
3337
        A[23] = -0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3336
 
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
3337
 
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3338
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.33333333*G8_;
 
3339
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.16666667*G8_;
3338
3340
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1;
3339
3341
        A[27] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1;
3340
3342
        A[28] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3341
3343
        A[29] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3342
 
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
3343
 
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3344
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.16666667*G8_;
 
3345
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3344
3346
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1;
3345
3347
        A[33] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1;
3346
3348
        A[34] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3352
3354
          A[0] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3353
3355
        A[1] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3354
3356
        A[2] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1;
3355
 
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.33333333*G8_;
 
3357
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3356
3358
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1;
3357
 
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3359
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3358
3360
        A[6] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3359
3361
        A[7] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3360
3362
        A[8] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1;
3361
 
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.16666667*G8_;
 
3363
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3362
3364
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1;
3363
 
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.33333333*G8_;
 
3365
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3364
3366
        A[12] = 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3365
3367
        A[13] = 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3366
3368
        A[14] = 0.00000000;
3367
3369
        A[15] = 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3368
3370
        A[16] = 0.00000000;
3369
3371
        A[17] = 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3370
 
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3371
 
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3372
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.33333333*G8_;
 
3373
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3372
3374
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1;
3373
3375
        A[21] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3374
3376
        A[22] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1;
3379
3381
        A[27] = 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3380
3382
        A[28] = 0.00000000;
3381
3383
        A[29] = 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3382
 
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
3383
 
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1;
 
3384
        A[30] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.16666667*G8_;
 
3385
        A[31] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_1_0_0 + 0.25000000*G5_1_1_1 + 0.33333333*G8_;
3384
3386
        A[32] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1;
3385
3387
        A[33] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 + 0.25000000*G4_1_2_0 + 0.25000000*G4_1_3_1;
3386
3388
        A[34] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1;
3392
3394
          A[0] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3393
3395
        A[1] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 - 0.25000000*G7_0_0_0 - 0.25000000*G7_0_1_1 - 0.25000000*G7_1_0_0 - 0.25000000*G7_1_1_1;
3394
3396
        A[2] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1;
3395
 
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.33333333*G8_;
3396
 
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1 + 0.16666667*G8_;
 
3397
        A[3] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
 
3398
        A[4] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 - 0.25000000*G6_0_2_0 - 0.25000000*G6_0_3_1 - 0.25000000*G6_1_2_0 - 0.25000000*G6_1_3_1;
3397
3399
        A[5] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1;
3398
3400
        A[6] = -0.25000000*G3_0_0_0 - 0.25000000*G3_0_1_1 - 0.25000000*G3_1_0_0 - 0.25000000*G3_1_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3399
3401
        A[7] = 0.25000000*G3_0_0_0 + 0.25000000*G3_0_1_1 + 0.25000000*G7_0_0_0 + 0.25000000*G7_0_1_1;
3400
3402
        A[8] = 0.25000000*G3_1_0_0 + 0.25000000*G3_1_1_1;
3401
 
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.16666667*G8_;
3402
 
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1 + 0.33333333*G8_;
 
3403
        A[9] = -0.25000000*G1_0_0_0 - 0.25000000*G1_0_1_1 - 0.25000000*G1_1_0_0 - 0.25000000*G1_1_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
 
3404
        A[10] = 0.25000000*G1_0_0_0 + 0.25000000*G1_0_1_1 + 0.25000000*G6_0_2_0 + 0.25000000*G6_0_3_1;
3403
3405
        A[11] = 0.25000000*G1_1_0_0 + 0.25000000*G1_1_1_1;
3404
3406
        A[12] = 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3405
3407
        A[13] = 0.25000000*G7_1_0_0 + 0.25000000*G7_1_1_1;
3407
3409
        A[15] = 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3408
3410
        A[16] = 0.25000000*G6_1_2_0 + 0.25000000*G6_1_3_1;
3409
3411
        A[17] = 0.00000000;
3410
 
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
3411
 
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1;
 
3412
        A[18] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.33333333*G8_;
 
3413
        A[19] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 - 0.25000000*G5_0_0_0 - 0.25000000*G5_0_1_1 - 0.25000000*G5_1_0_0 - 0.25000000*G5_1_1_1 + 0.16666667*G8_;
3412
3414
        A[20] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1;
3413
3415
        A[21] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3414
3416
        A[22] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 - 0.25000000*G4_0_2_0 - 0.25000000*G4_0_3_1 - 0.25000000*G4_1_2_0 - 0.25000000*G4_1_3_1;
3415
3417
        A[23] = 0.25000000*G0_1_2_0 + 0.25000000*G0_1_3_1;
3416
 
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
3417
 
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1;
 
3418
        A[24] = -0.25000000*G2_0_2_0 - 0.25000000*G2_0_3_1 - 0.25000000*G2_1_2_0 - 0.25000000*G2_1_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.16666667*G8_;
 
3419
        A[25] = 0.25000000*G2_0_2_0 + 0.25000000*G2_0_3_1 + 0.25000000*G5_0_0_0 + 0.25000000*G5_0_1_1 + 0.33333333*G8_;
3418
3420
        A[26] = 0.25000000*G2_1_2_0 + 0.25000000*G2_1_3_1;
3419
3421
        A[27] = -0.25000000*G0_0_2_0 - 0.25000000*G0_0_3_1 - 0.25000000*G0_1_2_0 - 0.25000000*G0_1_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3420
3422
        A[28] = 0.25000000*G0_0_2_0 + 0.25000000*G0_0_3_1 + 0.25000000*G4_0_2_0 + 0.25000000*G4_0_3_1;
3483
3485
  /// Return a string identifying the form
3484
3486
  virtual const char* signature() const
3485
3487
  {
3486
 
    return "Form([Integral(Product(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 0), Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 1)), Measure('exterior_facet', 0, None)), Integral(Sum(IndexSum(Product(Indexed(ComponentTensor(Product(FloatValue(0.5, (), (), {}), Indexed(Sum(NegativeRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 1), MultiIndex((Index(0),), {Index(0): 2})), MultiIndex((Index(0),), {Index(0): 2}))), PositiveRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 1), MultiIndex((Index(1),), {Index(1): 2})), MultiIndex((Index(1),), {Index(1): 2})))), MultiIndex((Index(2),), {Index(2): 2}))), MultiIndex((Index(2),), {Index(2): 2})), MultiIndex((Index(3),), {Index(3): 2})), Indexed(Sum(ComponentTensor(Product(Indexed(NegativeRestricted(VectorConstant(Cell('triangle', 1, Space(2)), 2, 0)), MultiIndex((Index(4),), {Index(4): 2})), NegativeRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 0))), MultiIndex((Index(4),), {Index(4): 2})), ComponentTensor(Product(Indexed(PositiveRestricted(VectorConstant(Cell('triangle', 1, Space(2)), 2, 0)), MultiIndex((Index(5),), {Index(5): 2})), PositiveRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 0))), MultiIndex((Index(5),), {Index(5): 2}))), MultiIndex((Index(3),), {Index(3): 2}))), MultiIndex((Index(3),), {Index(3): 2})), Sum(IndexSum(Product(Indexed(ComponentTensor(Product(FloatValue(0.5, (), (), {}), Indexed(Sum(NegativeRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 0), MultiIndex((Index(6),), {Index(6): 2})), MultiIndex((Index(6),), {Index(6): 2}))), PositiveRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 0), MultiIndex((Index(7),), {Index(7): 2})), MultiIndex((Index(7),), {Index(7): 2})))), MultiIndex((Index(8),), {Index(8): 2}))), MultiIndex((Index(8),), {Index(8): 2})), MultiIndex((Index(9),), {Index(9): 2})), Indexed(Sum(ComponentTensor(Product(Indexed(NegativeRestricted(VectorConstant(Cell('triangle', 1, Space(2)), 2, 0)), MultiIndex((Index(10),), {Index(10): 2})), NegativeRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 1))), MultiIndex((Index(10),), {Index(10): 2})), ComponentTensor(Product(Indexed(PositiveRestricted(VectorConstant(Cell('triangle', 1, Space(2)), 2, 0)), MultiIndex((Index(11),), {Index(11): 2})), PositiveRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 1))), MultiIndex((Index(11),), {Index(11): 2}))), MultiIndex((Index(9),), {Index(9): 2}))), MultiIndex((Index(9),), {Index(9): 2})), Product(NegativeRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 1)), PositiveRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', 1, Space(2)), 1), 0))))), Measure('interior_facet', 0, None))])";
 
3488
    return "Form([Integral(Product(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 0), Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 1)), Measure('exterior_facet', 0, None)), Integral(Sum(IndexSum(Product(Indexed(ComponentTensor(Product(FloatValue(0.50000000, (), (), {}), Indexed(Sum(NegativeRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 1), MultiIndex((Index(0),), {Index(0): 2})), MultiIndex((Index(0),), {Index(0): 2}))), PositiveRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 1), MultiIndex((Index(1),), {Index(1): 2})), MultiIndex((Index(1),), {Index(1): 2})))), MultiIndex((Index(2),), {Index(2): 2}))), MultiIndex((Index(2),), {Index(2): 2})), MultiIndex((Index(3),), {Index(3): 2})), Indexed(Sum(ComponentTensor(Product(Indexed(NegativeRestricted(VectorConstant(Cell('triangle', Space(2)), 2, 0)), MultiIndex((Index(4),), {Index(4): 2})), NegativeRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 0))), MultiIndex((Index(4),), {Index(4): 2})), ComponentTensor(Product(Indexed(PositiveRestricted(VectorConstant(Cell('triangle', Space(2)), 2, 0)), MultiIndex((Index(5),), {Index(5): 2})), PositiveRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 0))), MultiIndex((Index(5),), {Index(5): 2}))), MultiIndex((Index(3),), {Index(3): 2}))), MultiIndex((Index(3),), {Index(3): 2})), Sum(IndexSum(Product(Indexed(ComponentTensor(Product(FloatValue(0.50000000, (), (), {}), Indexed(Sum(NegativeRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 0), MultiIndex((Index(6),), {Index(6): 2})), MultiIndex((Index(6),), {Index(6): 2}))), PositiveRestricted(ComponentTensor(SpatialDerivative(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 0), MultiIndex((Index(7),), {Index(7): 2})), MultiIndex((Index(7),), {Index(7): 2})))), MultiIndex((Index(8),), {Index(8): 2}))), MultiIndex((Index(8),), {Index(8): 2})), MultiIndex((Index(9),), {Index(9): 2})), Indexed(Sum(ComponentTensor(Product(Indexed(NegativeRestricted(VectorConstant(Cell('triangle', Space(2)), 2, 0)), MultiIndex((Index(10),), {Index(10): 2})), NegativeRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 1))), MultiIndex((Index(10),), {Index(10): 2})), ComponentTensor(Product(Indexed(PositiveRestricted(VectorConstant(Cell('triangle', Space(2)), 2, 0)), MultiIndex((Index(11),), {Index(11): 2})), PositiveRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 1))), MultiIndex((Index(11),), {Index(11): 2}))), MultiIndex((Index(9),), {Index(9): 2}))), MultiIndex((Index(9),), {Index(9): 2})), Product(NegativeRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 0)), PositiveRestricted(Argument(FiniteElement('Discontinuous Lagrange', Cell('triangle', Space(2)), 1, None), 1))))), Measure('interior_facet', 0, None))])";
3487
3489
  }
3488
3490
 
3489
3491
  /// Return the rank of the global tensor (r)