← Back to branch summary

~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to dolfin/mesh/ALE.cpp

  • Committer: ilmarw@gogmagog.simula.no
  • Date: 2008-05-29 10:46:14 UTC
  • mfrom: (2122.3.2 trunk)
  • mto: (2668.1.45 trunk)
  • mto: This revision was merged to the branch mainline in revision 2670.
  • Revision ID: ilmarw@gogmagog.simula.no-20080529104614-xjn0a188wtj2wlmw
Repo fix, roll back unwanted Parameter changes.

Show diffs side-by-side

added added

removed removed

Lines of Context:
dolfin/mesh/ALE.cpp
2
2
// Licensed under the GNU LGPL Version 2.1.
3
3
//
4
4
// First added:  2008-05-02
5
 
// Last changed: 2008-05-29
 
5
// Last changed: 2008-05-28
6
6
 
7
7
#include <string.h>
8
8
#include <dolfin/common/Array.h>
37
37
  // If hermite, create dgdn
38
38
  if (method == hermite)
39
39
  {
40
 
    cout<<"hermite"<<endl;
41
40
    hermiteFunction(ghat, dim, new_boundary,
42
41
                    mesh, *vertex_map, *cell_map);
43
42
  }
235
234
  normals(dfdn, dim, new_boundary,
236
235
          mesh, vertex_map, cell_map);
237
236
 
238
 
  real c;
239
 
  if (dim==2)
240
 
    c=2;
241
 
  else
242
 
    c=M_PI;
243
 
 
244
 
  //FAKTOREN c før dfdn, HVA VELGER VI DER?
 
237
  //FAKTOREN før dfdn, HVA VELGER VI DER?
245
238
  for (VertexIterator v(new_boundary); !v.end(); ++v) {
246
239
    ghat[v->index()]=new real [dim];
247
240
    integral(ghat[v->index()], dim, new_boundary,
248
241
             mesh, vertex_map, *v);
249
242
    for (uint i=0; i<dim;i++) 
250
 
      ghat[v->index()][i]=c*dfdn[v->index()][i]-ghat[v->index()][i];
 
243
      ghat[v->index()][i]=2*dfdn[v->index()][i]-ghat[v->index()][i];
251
244
  }
252
245
  for (uint i=0; i<new_boundary.numVertices(); i++)
253
246
    delete [] dfdn[i];
280
273
          Point n=mesh_cell.normal(facet_index);
281
274
          
282
275
          for (uint j=0; j<dim; j++)
283
 
            dfdn[ind][j]-=n[j];
 
276
            dfdn[ind][j]+=n[j];
284
277
          break;
285
278
        }
286
279
      }