Dear Daniel,
I will try to come up with a minimal example. In the meantime, below is the
piece of code and I will try to explain using this. I tried running in
Debug mode as well but not able to figure out the error. For the code
below in combination with the rest of the code, the last line of output on
the console is *Checkpoint 1*.
Note: I am using dealii as an external library and am writing another small
library for my own purpose. Objects such as *mesh* in the piece of code
below refer to my library.
This is a vector valued problem with 2 dofs per support point. *dof_handler
is a pointer for the dof_handler object of hp::DoFHandler<dim> type for a
rectangular triangulation with 8X4 elements of Q2 type and using Q1
mapping. system_rhs is a vector of size 17 * 9 * 2 = 306.
*template <int dim>void FEM<dim>::add_point_source_to_rhs(){
deallog.depth_console (2); *
* unsigned int no_sources = (mesh->point_source_vector).size();
for(unsigned int s = 0; s < no_sources; ++s){ std::vector<double>
load_point = mesh->point_source_vector[s].first; std::vector<double>
load = mesh->point_source_vector[s].second; Assert
(load_point.size() == dim,
ExcDimensionMismatch(load_point.size(), dim)); Assert (load.size()
== dim, ExcDimensionMismatch(load.size(), dim));*
//defining the load and the load point
* Point<dim> ldp, ld; for(unsigned int i = 0; i < dim;
++i){ ldp(i) = load_point[i]; ld(i) = load[i];
} *
* std::cout<<"Checkpoint 1"<<std::endl;
VectorTools::create_point_source_vector(*dof_handler, ldp, ld,
system_rhs); std::cout<<"Checkpoint 2"<<std::endl; }}*
Are you able to figure out any bug with this information?
Thanks!
Best
Deepak
On Thu, Jul 28, 2016 at 2:45 PM, Daniel Arndt <
[email protected]> wrote:
> Deepak,
>
> could you provide us with a minimal example that shows your problem?
> Are you running in debug mode? What exactly does the program tell you when
> aborting?
>
> Best,
> Daniel
>
> --
> 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 [email protected].
> For more options, visit https://groups.google.com/d/optout.
>
--
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
