← Back to branch summary

~ubuntu-branches/ubuntu/wily/rheolef/wily

« back to all changes in this revision

Viewing changes to nfem/ptst/form_mass_bdr_tst.cc

  • Committer: Package Import Robot
  • Author(s): Pierre Saramito, Pierre Saramito, Sylvestre Ledru
  • Date: 2013-04-30 11:29:51 UTC
  • mfrom: (1.2.1) (11.1.1 experimental)
  • Revision ID: package-import@ubuntu.com-20130430112951-g69nuqin77t75dek
Tags: 6.4-1
[ Pierre Saramito ]
* New upstream release 6.4 (major changes):
  - new powerful c++ expressions for FEM specifications
  - configure script improved (portability increased)
  - minor bug fixes
* control:
  - all architectures are considered (armel & armhf are back)
  - no more circular-dependency librheolef-dev <--> rheolef
* rules: add security compiler flags (hardening and -W,-l,relro)
* watch: include an updated version, thanks to B. Martens <bartm@debian.org>
* copyright: various file format fixes

[ Sylvestre Ledru ]
*  Standards-Version updated to 3.9.4

Show diffs side-by-side

added added

removed removed

Lines of Context:
nfem/ptst/form_mass_bdr_tst.cc
298
298
    // ---------------------------------------------------
299
299
    // first check : compute boundary area with i-th weight
300
300
    // ---------------------------------------------------
301
 
    Float mes_gamma_1 = dot(one, m_h*weight_h);
 
301
    Float mes_gamma_1 = dual(one, m_h*weight_h);
302
302
    cout << "mes1(gamma) = " << double(mes_gamma_1) << endl;
303
303
    // ---------------------------------------------------
304
304
    // 2nd check : use trace
327
327
        //
328
328
        // compute boundary area with i-th weight
329
329
        //
330
 
        Float mes_gamma_2 = dot(one_bdr, m_gamma_h*weight_bdr);
 
330
        Float mes_gamma_2 = dual(one_bdr, m_gamma_h*weight_bdr);
331
331
        cout << "mes2(gamma) = " << double(mes_gamma_2) << endl;
332
332
        // ---------------------------------------------------
333
333
        // 3nd check : compare two results
334
334
        // ---------------------------------------------------
335
335
        if (fabs(mes_gamma_1-mes_gamma_2) > sqrt(epsilon)) {
336
 
          derr << "boundary mass problem detected." << endl;
 
336
          derr << "mes1(gamma) != mes2(gamma): boundary mass problem detected." << endl;
337
337
          status = 1;
338
338
        }
339
339
    }