Dear Maxi, Aligned with Wolfgang's thoughts (and considering your recent flurry of posts), it seems like you need to slow down and more carefully consider what you are doing in your work. We all have to work through bugs and conceptual issues - this is part of learning both the deal.II, or any library for that matter. In 99% of cases the documentation (manual + tutorials) does a good job of explaining not only what any particular function does but also how its used in context. Its up to you to ensure that you've done your best to understand that context. Where it fails to be help, we're here and happy to clarify. Not all of the information is in the tutorials that may be directly relevant to your work (such is the philosophy of our tutorial structure) - one often has to do a bit of leg-work to find it.

So, I'm not going to answer your question directly, but rather I'm going to pose a couple of questions to you: 1. In this line VectorTools::project(dof_handler, constraint_matrix, QGauss<dim>(3), InitialValues<dim>(carrier_density), solution); where is the projected function being output/stored? 2. What type of object does fe_values[N_density] return? 3. What information are you asking for here? fe_values[N_density].value(i, q_point) 4. Are what you doing in (1) and (3) in any way related? If not, then why not? What is the disconnect? I'd be happy to direct you to an example in the tutorials that, to the best of my understanding of what you're trying to achieve, demonstrates the mechanism to extract the information that you're trying to get at. Best, Jean-Paul On Saturday, August 5, 2017 at 10:18:44 AM UTC+2, Maxi Miller wrote: > > I tried to debug it (thus the reason for "volatile"), and reduced the > whole code down to (called in SetupSystem()): > > > VectorTools::project(dof_handler, constraint_matrix, QGauss<dim>(3), > InitialValues<dim>(carrier_density), solution); > //For checking of content > const FEValuesExtractors::Vector N_density(N_component); > const FEValuesExtractors::Vector E_temperature(TE_component); > const FEValuesExtractors::Vector L_temperature(TL_component); > const QGauss<dim> quadrature_formula(fe_degree+1); > const unsigned int dofs_per_cell = finite_element.dofs_per_cell; > const unsigned int n_q_points = quadrature_formula.size(); > FEValues<dim> fe_values(finite_element, quadrature_formula, > update_values | update_gradients | > update_JxW_values); > std::vector<Tensor<1, dim>> value_N(dofs_per_cell), value_TE( > dofs_per_cell), value_TL(dofs_per_cell); > for (auto cell = dof_handler.begin_active(); cell!=dof_handler.end > (); ++cell) > { > > fe_values.reinit (cell); > > for (unsigned int q_point=0; q_point<n_q_points; ++q_point) > { > for(size_t i = 0; i < dofs_per_cell; ++i) > { > std::cout << "Value_N[" << i << "] = " << fe_values[ > N_density].value(i, q_point) << '\n'; > getchar(); > value_N[i] = fe_values[N_density].value(i, q_point); > value_TE[i] = fe_values[E_temperature].value(i, > q_point); > value_TL[i] = fe_values[L_temperature].value(i, > q_point); > } > } > } > //Content checking end > Furthermore, I checked (via break points) how the values are set in > InitialValues<>: > > template <int dim> > double InitialValues<dim>::value(const Point<dim> &p, const unsigned > int component) const > { > (void) p; > Assert(component < this->n_components, ExcIndexRange(component, 0, > this->n_components)); > if(component < dim) > return equations.initial_carrier_density; > else > if(component < 2*dim && component >= dim) > return equations.ambient_temperature; > else > return equations.ambient_temperature; > } > > template <int dim> > void InitialValues<dim>::vector_value(const Point<dim> &p, Vector > <double> &values) const > { > for(size_t i = 0; i < this->n_components; ++i) > values(i) = InitialValues<dim>::value(p, i); > } > > and according to the debug output, every value was set correctly in that > function. Nevertheless in the debug output code above I still get (for the > first point, and the first values of i and j: > Value_N[0] = 0.324665 0 0 > > Value_N[1] = 0 0.324665 0 > > Value_N[2] = 0 0 0.324665 > > Value_N[3] = 0 0 0 > > Value_N[4] = 0 0 0 > > Value_N[5] = 0 0 0 > > Value_N[6] = 0 0 0 > > Value_N[7] = 0 0 0 > > Value_N[8] = 0 0 0 > > Value_N[9] = -0.0412379 0 0 > > Value_N[10] = 0 -0.0412379 0 > > Value_N[11] = 0 0 -0.0412379 > > Value_N[12] = 0 0 0 > > Value_N[13] = 0 0 0 > > Value_N[14] = 0 0 0 > > Value_N[15] = 0 0 0 > > Value_N[16] = 0 0 0 > > Value_N[17] = 0 0 0 > > Thus: Is there a misunderstanding in how the data is stored, or how the > data is processed? According to my understanding the values should be 1e18 > at every point. Is that correct? Or is there a better way to debug my code? > > > Am Freitag, 4. August 2017 23:48:08 UTC+2 schrieb Wolfgang Bangerth: >> >> On 08/04/2017 09:49 AM, Maxi Miller wrote: >> > >> > Why is there only one value not equal to zero, and why is it moving? >> Shouldn't >> > all the values be set already? >> >> Maxi -- this is not a question to which you can expect an answer on the >> mailing list. We are happy to help with conceptual questions, but you do >> need >> to debug your code yourself. >> >> (I will note that there is almost never a reason to mark variables as >> 'volatile' as you do. It's not wrong -- it just doesn't have any effect >> in >> this context at all.) >> >> Best >> W. >> >> -- >> ------------------------------------------------------------------------ >> Wolfgang Bangerth email: bang...@colostate.edu >> www: http://www.math.colostate.edu/~bangerth/ >> >> -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.