← Back to branch summary

~fenics-core/ffc/spatder2grad

« back to all changes in this revision

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

  • Committer: Marie E. Rognes
  • Date: 2013-01-08 12:10:22 UTC
  • mfrom: (1760.1.82 ffc-manifolds-trunk)
  • Revision ID: meg@simula.no-20130108121022-svq8jno62hvwp8eu
Merge work on code generation for manifolds. You can now define
elements and forms over cells with different geometric and topological
dimensions.

Show diffs side-by-side

added added

removed removed

Lines of Context:
test/regression/references/r_quadrature/Poisson.h
107
107
    const double J_11 = x[2][1] - x[0][1];
108
108
    
109
109
    // Compute determinant of Jacobian
110
 
    double detJ = J_00*J_11 - J_01*J_10;
 
110
    const double detJ = J_00*J_11 - J_01*J_10;
111
111
    
112
112
    // Compute inverse of Jacobian
113
113
    
244
244
    const double J_11 = x[2][1] - x[0][1];
245
245
    
246
246
    // Compute determinant of Jacobian
247
 
    double detJ = J_00*J_11 - J_01*J_10;
 
247
    const double detJ = J_00*J_11 - J_01*J_10;
248
248
    
249
249
    // Compute inverse of Jacobian
250
250
    const double K_00 =  J_11 / detJ;
1210
1210
    const double J_11 = x[2][1] - x[0][1];
1211
1211
    
1212
1212
    // Compute determinant of Jacobian
1213
 
    double detJ = J_00*J_11 - J_01*J_10;
 
1213
    const double detJ = J_00*J_11 - J_01*J_10;
1214
1214
    
1215
1215
    // Compute inverse of Jacobian
1216
1216
    const double K_00 =  J_11 / detJ;
1221
1221
    // Set scale factor
1222
1222
    const double det = std::abs(detJ);
1223
1223
    
1224
 
    // Cell Volume.
 
1224
    // Cell volume.
1225
1225
    
1226
 
    // Compute circumradius, assuming triangle is embedded in 2D.
 
1226
    // Compute circumradius of triangle in 2D.
1227
1227
    
1228
1228
    
1229
1229
    // Facet Area.
1312
1312
    const double J_11 = x[2][1] - x[0][1];
1313
1313
    
1314
1314
    // Compute determinant of Jacobian
1315
 
    double detJ = J_00*J_11 - J_01*J_10;
 
1315
    const double detJ = J_00*J_11 - J_01*J_10;
1316
1316
    
1317
1317
    // Compute inverse of Jacobian
1318
1318
    
1319
1319
    // Set scale factor
1320
1320
    const double det = std::abs(detJ);
1321
1321
    
1322
 
    // Cell Volume.
 
1322
    // Cell volume.
1323
1323
    
1324
 
    // Compute circumradius, assuming triangle is embedded in 2D.
 
1324
    // Compute circumradius of triangle in 2D.
1325
1325
    
1326
1326
    
1327
1327
    // Facet Area.