[dune-pdelab] Questions about linearelasticity.hh

[Oliver Sander]

Dear pdelab team,
I have another question about linearelasticity.hh.  This one is a bit more
serious.

The file contains the method lambda_boundary which, IIRC computes some sort
of Neumann boundary values.  The method contains the lines

          if (global[0] != 10.0)
            return;
          std::cout << "BC" << std::endl;

and

          y[dim-1] = 0.001;

which are obviously bogus -- they implement somebody's special problem.
Please remove those lines and any other personal cruft that may be lurking
in that file.

<rant>How am I supposed to trust a discretization framework that contains
such pitfalls?</rant>

thanks,
Oliver

-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 551 bytes
Desc: OpenPGP digital signature
URL: <https://lists.dune-project.org/pipermail/dune-pdelab/attachments/20140104/1ad82f81/attachment.sig>