Ahmad, > thank you very much, but I do not know, it did not gave me any error in > the compilation, but unfortunately it gave me the following error when I > ran it: > > ============================ Running step-4 > Solving with Q1 elements, adaptive refinement > ============================================= > > -------------------------------------------------------- > An error occurred in line <4206> of file > </home/Ahmad/deal.II/deal.II/include/fe/fe_values.h> in function const > dealii::Point<spacedim>& dealii::FEValuesBase<dim, > spacedim>::quadrature_point(unsigned int) const [with int dim = 2, int > spacedim = 2] The violated condition was: > i<this->quadrature_points.size() > The name and call sequence of the exception was: > ExcIndexRange(i, 0, this->quadrature_points.size()) > Additional Information: > Index 3 is not in [0,3[ > > Stacktrace: > ----------- > #0 ./step-4: dealii::FEValuesBase<2, 2>::quadrature_point(unsigned int) > const #1 ./step-4: HelmholtzProblem<2>::assemble_system() > #2 ./step-4: HelmholtzProblem<2>::run() > #3 ./step-4: main > -------------------------------------------------------- > make: *** [run] Aborted
I assume you probably solved this one already but if not: the error happens because you ask the FEValues or FEFaceValues for a quadrature point with an index that is larger than the number of quadrature points you have. You need to check the upper bound of your loop over all quadrature points. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
