<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<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 <klaus.kaiser@rwth-aachen.de> 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>
> <br>
> <br>
> On 10/24/2012 06:48 PM, Klaus Kaiser wrote:<br>
>> Hi,<br>
>> <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>
>> <br>
> This sounds to me like Lagrange Interpolation, for a DG Space a L2 projection would be the suitable interpolation method.<br>
> <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>
> <br>
> Some discreteFunction implementations allows you to access the DofVector as double*, eg. AdaptiveDiscreteFunction ( leakPointer() ).<br>
> <br>
> Let me ask you where do you need the global Dof Vector?<br>
> <br>
> Best Tobias<br>
> <br>
>> Greetings and thanks for your help<br>
>> <br>
>> Klaus<br>
>> <br>
>> _______________________________________________<br>
>> dune-fem mailing list<br>
>> dune-fem@dune-project.org<br>
>> <a href="http://lists.dune-project.org/mailman/listinfo/dune-fem">http://lists.dune-project.org/mailman/listinfo/dune-fem</a><br>
>> <br>
> <br>
> <br>
> -- <br>
> Tobias Malkmus <tomalk at mathematik.uni-freiburg.de><br>
> <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>
> <br>
<br>
<br>
_______________________________________________<br>
dune-fem mailing list<br>
dune-fem@dune-project.org<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>
</body>
</html>