← Back to branch summary

~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to dolfin/fem/Assembler.cpp

  • Committer: Anders Logg
  • Date: 2008-06-09 21:06:51 UTC
  • mto: (2668.1.50 trunk)
  • mto: This revision was merged to the branch mainline in revision 2670.
  • Revision ID: logg@simula.no-20080609210651-xrycs6vcqnrpb3as
Reset local tensor A^K before call to tabulate_tensor()

Show diffs side-by-side

added added

removed removed

Lines of Context:
dolfin/fem/Assembler.cpp
5
5
// Modified by Ola Skavhaug, 2007
6
6
//
7
7
// First added:  2007-01-17
8
 
// Last changed: 2008-05-28
 
8
// Last changed: 2008-06-09
9
9
 
10
10
#include <dolfin/log/dolfin_log.h>
11
11
#include <dolfin/common/Array.h>
190
190
      dof_map_set[i].tabulate_dofs(ufc.dofs[i], ufc.cell, cell->index());
191
191
 
192
192
    // Tabulate cell tensor
 
193
    ufc.reset();
193
194
    integral->tabulate_tensor(ufc.A, ufc.w, ufc.cell);
194
195
 
195
196
    // Add entries to global tensor
256
257
      dof_map_set[i].tabulate_dofs(ufc.dofs[i], ufc.cell, mesh_cell.index());
257
258
 
258
259
    // Tabulate exterior facet tensor
 
260
    ufc.reset();
259
261
    integral->tabulate_tensor(ufc.A, ufc.w, ufc.cell, local_facet);
260
 
 
 
262
    
261
263
    // Add entries to global tensor
262
264
    A.add(ufc.A, ufc.local_dimensions, ufc.dofs);
263
265
 
332
334
    }
333
335
 
334
336
    // Tabulate exterior interior facet tensor on macro element
 
337
    ufc.reset();
335
338
    integral->tabulate_tensor(ufc.macro_A, ufc.macro_w, ufc.cell0, ufc.cell1, facet0, facet1);
336
339
 
337
340
    // Add entries to global tensor