<html><head><base href="x-msg://4/"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">ah ok, thanks for this, I thought you were using other basis functions, then everything is clear ;-)<div><br></div><div>thanks</div><div><br></div><div>Klaus</div><div><br></div><div><br><div><div>Am 25.10.2012 um 13:24 schrieb Dedner, Andreas:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div><div>Two remarks <div>If you are using a dg lagrange space then 1 0 0 0 would not be correct - locally your discrete function would would be the first lagrange hat function on the entity - 1at node with local number 0 whichever that is and zero at the other nodes. </div><div>If we are talking Legende dg space then your local dofs are correct I think. </div><div><br></div><div>The local function can act like a vector for the local dof. Is that not good enough for your application? </div><div>Greetings </div><div>Andreas </div><div>.</div><br>Klaus Kaiser <<a href="mailto:klaus.kaiser@rwth-aachen.de">klaus.kaiser@rwth-aachen.de</a>> wrote:<br></div><font size="2"><span style="font-size: 10pt; "><div class="PlainText">Hi Tobias,<br><br>thanks for your fast answer. I have done the whole interpolation with LagrangeDiscontinuousGalerkinSpace and AdaptiveDiscreteFunction and set on every entity the dofs with<br><br>f.localFunction( *it )[j]=1;<br><br>where f is my AdaptiveDiscreteFunction. When I understood it right, then I iterate through every lagrange Point (DiscreteDGSpaceType::LagrangePointSetType) and set the value of the function at this point. So for example I got an constant Function (f(x)=1) and Polynomial Order 5, then on an entity the interpolation should look like this:<br><br>f_(x) = 1*phi_0(x)+0*phi_1(x)+0*phi_2(x)…<br><br>where phi_i is the i-the lagrange polynomial. With the function leakPointer(), I get an pointer on an array which gives me the functionvalues at the lagrange points. But I I think you misunderstood my question, I don't want to access the global DoF Vector, I want to access the coefficients of the basefunctions of the local functions. So in the example above I want to access an vector like this:<br><br>(1 0 0 0 0 0)<br><br>How do I get this vector? Or is this impossible?<br><br>Greetings & thanks<br><br>Klaus<br><br><br>Am 25.10.2012 um 08:20 schrieb Tobias Malkmus:<br><br>> Hi Klaus<br>><span class="Apple-converted-space"> </span><br>><span class="Apple-converted-space"> </span><br>> On 10/24/2012 06:48 PM, Klaus Kaiser wrote:<br>>> Hi,<br>>><span class="Apple-converted-space"> </span><br>>> I have question concerning the dune-fem/dg part. I created an interpolation of a function in the discrete discontinuous galerkin space. This was done by setting the DoF Vector on every entity with the corresponding values of the function. Now I want to access the coefficients of the base functions, how do I do that? I searched a lot in the doxygen documentation, but didn't find a function, which gives me the values.<br>>><span class="Apple-converted-space"> </span><br>> This sounds to me like Lagrange Interpolation, for a DG Space a L2 projection would be the suitable interpolation method.<br>><span class="Apple-converted-space"> </span><br>> In Dune-Fem a discrete function is a combination of DofVector + DofMapper + LocalBasefunctions.<br>> This allows you to generate localFunctions which are the local representation of a GridFunction / DiscreteFunction.<br>><span class="Apple-converted-space"> </span><br>> Some discreteFunction implementations allows you to access the DofVector as double*, eg. AdaptiveDiscreteFunction ( leakPointer() ).<br>><span class="Apple-converted-space"> </span><br>> Let me ask you where do you need the global Dof Vector?<br>><span class="Apple-converted-space"> </span><br>> Best Tobias<br>><span class="Apple-converted-space"> </span><br>>> Greetings and thanks for your help<br>>><span class="Apple-converted-space"> </span><br>>> Klaus<br>>><span class="Apple-converted-space"> </span><br>>> _______________________________________________<br>>> dune-fem mailing list<br>>> <a href="mailto:dune-fem@dune-project.org">dune-fem@dune-project.org</a><br>>><span class="Apple-converted-space"> </span><a href="http://lists.dune-project.org/mailman/listinfo/dune-fem">http://lists.dune-project.org/mailman/listinfo/dune-fem</a><br>>><span class="Apple-converted-space"> </span><br>><span class="Apple-converted-space"> </span><br>><span class="Apple-converted-space"> </span><br>> --<span class="Apple-converted-space"> </span><br>> Tobias Malkmus <tomalk at <a href="http://mathematik.uni-freiburg.de">mathematik.uni-freiburg.de</a>><br>><span class="Apple-converted-space"> </span><br>> Mathematisches Institut Tel: +49 (0) 761 203 5651<br>> Abt. für Angewandte Mathematik Universität Freiburg<br>> Hermann-Herder-Str. 10<br>> 79104 Freiburg<br>><span class="Apple-converted-space"> </span><br><br><br>_______________________________________________<br>dune-fem mailing list<br><a href="mailto:dune-fem@dune-project.org">dune-fem@dune-project.org</a><br><a href="http://lists.dune-project.org/mailman/listinfo/dune-fem">http://lists.dune-project.org/mailman/listinfo/dune-fem</a><br><br></div></span></font></div></span></blockquote></div><br></div></body></html>