← Back to branch summary

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

« back to all changes in this revision

Viewing changes to test/regression/references/r_quadrature/TensorWeightedPoisson.h

  • Committer: Package Import Robot
  • Author(s): Johannes Ring
  • Date: 2011-10-26 17:52:20 UTC
  • mfrom: (1.1.10)
  • Revision ID: package-import@ubuntu.com-20111026175220-ope1dzqv4jn2b8pq
Tags: 1.0-beta2-1
* New upstream release. This release includes some performance
  improvements for evaluating basis functions. It also adds support
  for Bessel functions and error functions.
* debian/control: Bump version numbers for python-ufc, python-fiat,
  python-instant, and python-ufl.

Show diffs side-by-side

added added

removed removed

Lines of Context:
test/regression/references/r_quadrature/TensorWeightedPoisson.h
1
 
// This code conforms with the UFC specification version 2.0.0
2
 
// and was automatically generated by FFC version 0.9.10.
 
1
// This code conforms with the UFC specification version 2.0.2
 
2
// and was automatically generated by FFC version 1.0-beta+.
3
3
// 
4
4
// This code was generated with the following parameters:
5
5
// 
1595
1595
      double basisvalues[3] = {0.0, 0.0, 0.0};
1596
1596
      
1597
1597
      // Declare helper variables.
1598
 
      unsigned int rr = 0;
1599
 
      unsigned int ss = 0;
1600
1598
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1601
1599
      
1602
1600
      // Compute basisvalues.
1603
1601
      basisvalues[0] = 1.0;
1604
1602
      basisvalues[1] = tmp0;
1605
 
      for (unsigned int r = 0; r < 1; r++)
1606
 
      {
1607
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1608
 
        ss = r*(r + 1)/2;
1609
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1610
 
      }// end loop over 'r'
1611
 
      for (unsigned int r = 0; r < 2; r++)
1612
 
      {
1613
 
        for (unsigned int s = 0; s < 2 - r; s++)
1614
 
        {
1615
 
          rr = (r + s)*(r + s + 1)/2 + s;
1616
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1617
 
        }// end loop over 's'
1618
 
      }// end loop over 'r'
 
1603
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1604
      basisvalues[0] *= std::sqrt(0.5);
 
1605
      basisvalues[2] *= std::sqrt(1.0);
 
1606
      basisvalues[1] *= std::sqrt(3.0);
1619
1607
      
1620
1608
      // Table(s) of coefficients.
1621
1609
      static const double coefficients0[3] = \
1635
1623
      double basisvalues[3] = {0.0, 0.0, 0.0};
1636
1624
      
1637
1625
      // Declare helper variables.
1638
 
      unsigned int rr = 0;
1639
 
      unsigned int ss = 0;
1640
1626
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1641
1627
      
1642
1628
      // Compute basisvalues.
1643
1629
      basisvalues[0] = 1.0;
1644
1630
      basisvalues[1] = tmp0;
1645
 
      for (unsigned int r = 0; r < 1; r++)
1646
 
      {
1647
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1648
 
        ss = r*(r + 1)/2;
1649
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1650
 
      }// end loop over 'r'
1651
 
      for (unsigned int r = 0; r < 2; r++)
1652
 
      {
1653
 
        for (unsigned int s = 0; s < 2 - r; s++)
1654
 
        {
1655
 
          rr = (r + s)*(r + s + 1)/2 + s;
1656
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1657
 
        }// end loop over 's'
1658
 
      }// end loop over 'r'
 
1631
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1632
      basisvalues[0] *= std::sqrt(0.5);
 
1633
      basisvalues[2] *= std::sqrt(1.0);
 
1634
      basisvalues[1] *= std::sqrt(3.0);
1659
1635
      
1660
1636
      // Table(s) of coefficients.
1661
1637
      static const double coefficients0[3] = \
1675
1651
      double basisvalues[3] = {0.0, 0.0, 0.0};
1676
1652
      
1677
1653
      // Declare helper variables.
1678
 
      unsigned int rr = 0;
1679
 
      unsigned int ss = 0;
1680
1654
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1681
1655
      
1682
1656
      // Compute basisvalues.
1683
1657
      basisvalues[0] = 1.0;
1684
1658
      basisvalues[1] = tmp0;
1685
 
      for (unsigned int r = 0; r < 1; r++)
1686
 
      {
1687
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1688
 
        ss = r*(r + 1)/2;
1689
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1690
 
      }// end loop over 'r'
1691
 
      for (unsigned int r = 0; r < 2; r++)
1692
 
      {
1693
 
        for (unsigned int s = 0; s < 2 - r; s++)
1694
 
        {
1695
 
          rr = (r + s)*(r + s + 1)/2 + s;
1696
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1697
 
        }// end loop over 's'
1698
 
      }// end loop over 'r'
 
1659
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1660
      basisvalues[0] *= std::sqrt(0.5);
 
1661
      basisvalues[2] *= std::sqrt(1.0);
 
1662
      basisvalues[1] *= std::sqrt(3.0);
1699
1663
      
1700
1664
      // Table(s) of coefficients.
1701
1665
      static const double coefficients0[3] = \
1834
1798
      double basisvalues[3] = {0.0, 0.0, 0.0};
1835
1799
      
1836
1800
      // Declare helper variables.
1837
 
      unsigned int rr = 0;
1838
 
      unsigned int ss = 0;
1839
1801
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1840
1802
      
1841
1803
      // Compute basisvalues.
1842
1804
      basisvalues[0] = 1.0;
1843
1805
      basisvalues[1] = tmp0;
1844
 
      for (unsigned int r = 0; r < 1; r++)
1845
 
      {
1846
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
1847
 
        ss = r*(r + 1)/2;
1848
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
1849
 
      }// end loop over 'r'
1850
 
      for (unsigned int r = 0; r < 2; r++)
1851
 
      {
1852
 
        for (unsigned int s = 0; s < 2 - r; s++)
1853
 
        {
1854
 
          rr = (r + s)*(r + s + 1)/2 + s;
1855
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
1856
 
        }// end loop over 's'
1857
 
      }// end loop over 'r'
 
1806
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1807
      basisvalues[0] *= std::sqrt(0.5);
 
1808
      basisvalues[2] *= std::sqrt(1.0);
 
1809
      basisvalues[1] *= std::sqrt(3.0);
1858
1810
      
1859
1811
      // Table(s) of coefficients.
1860
1812
      static const double coefficients0[3] = \
1992
1944
      double basisvalues[3] = {0.0, 0.0, 0.0};
1993
1945
      
1994
1946
      // Declare helper variables.
1995
 
      unsigned int rr = 0;
1996
 
      unsigned int ss = 0;
1997
1947
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
1998
1948
      
1999
1949
      // Compute basisvalues.
2000
1950
      basisvalues[0] = 1.0;
2001
1951
      basisvalues[1] = tmp0;
2002
 
      for (unsigned int r = 0; r < 1; r++)
2003
 
      {
2004
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
2005
 
        ss = r*(r + 1)/2;
2006
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
2007
 
      }// end loop over 'r'
2008
 
      for (unsigned int r = 0; r < 2; r++)
2009
 
      {
2010
 
        for (unsigned int s = 0; s < 2 - r; s++)
2011
 
        {
2012
 
          rr = (r + s)*(r + s + 1)/2 + s;
2013
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
2014
 
        }// end loop over 's'
2015
 
      }// end loop over 'r'
 
1952
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
1953
      basisvalues[0] *= std::sqrt(0.5);
 
1954
      basisvalues[2] *= std::sqrt(1.0);
 
1955
      basisvalues[1] *= std::sqrt(3.0);
2016
1956
      
2017
1957
      // Table(s) of coefficients.
2018
1958
      static const double coefficients0[3] = \
2150
2090
      double basisvalues[3] = {0.0, 0.0, 0.0};
2151
2091
      
2152
2092
      // Declare helper variables.
2153
 
      unsigned int rr = 0;
2154
 
      unsigned int ss = 0;
2155
2093
      double tmp0 = (1.0 + Y + 2.0*X)/2.0;
2156
2094
      
2157
2095
      // Compute basisvalues.
2158
2096
      basisvalues[0] = 1.0;
2159
2097
      basisvalues[1] = tmp0;
2160
 
      for (unsigned int r = 0; r < 1; r++)
2161
 
      {
2162
 
        rr = (r + 1)*(r + 1 + 1)/2 + 1;
2163
 
        ss = r*(r + 1)/2;
2164
 
        basisvalues[rr] = basisvalues[ss]*(0.5 + r + Y*(1.5 + r));
2165
 
      }// end loop over 'r'
2166
 
      for (unsigned int r = 0; r < 2; r++)
2167
 
      {
2168
 
        for (unsigned int s = 0; s < 2 - r; s++)
2169
 
        {
2170
 
          rr = (r + s)*(r + s + 1)/2 + s;
2171
 
          basisvalues[rr] *= std::sqrt((0.5 + r)*(1.0 + r + s));
2172
 
        }// end loop over 's'
2173
 
      }// end loop over 'r'
 
2098
      basisvalues[2] = basisvalues[0]*(0.5 + 1.5*Y);
 
2099
      basisvalues[0] *= std::sqrt(0.5);
 
2100
      basisvalues[2] *= std::sqrt(1.0);
 
2101
      basisvalues[1] *= std::sqrt(3.0);
2174
2102
      
2175
2103
      // Table(s) of coefficients.
2176
2104
      static const double coefficients0[3] = \