Dear Daniel,
Thank you for you replies.
The same error happens with the automatically configured version, but on
64%:
[ 64%] Building CXX object source/CMakeFiles/deal_II.dir/dummy.cc.o
Linking CXX shared library ../lib/libdeal_II.so
collect2: error: ld terminated with signal 7 [Bus error],
Dear Wolfgang,
I do compute the solution at many points along the profile, and while Ex and Ey
look reasonable, the Ez does not. Is it possible that there is a problem in
point_value function for z component for Nedelec elements?
Anna
--
The deal.II project is located at
Dear Wolfgang,
Thank you for your replies.
I am still struggling with outputs of the solution to the file.
I have two type of outputs: 1) output vtu files for each processor on the
whole mesh and 2) creating a simple txt file with values of the solution at
receiver locations
1) For the first
Dear Timo and Daniel,
I have tried to submit the job, but it stopped again at the same place:
[ 57%] Linking CXX shared library ../lib/libdeal_II.so
Attached is the error file.
Any ideas?
Anna
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
Dear Daniel,
If I type "free -h" I get the following info:
totalusedfree shared buff/cache
available
Mem: 251G 44G 23G 11G182G
164G
Swap: 11G535M 11G
I can try to submit to the
Dear all,
I am trying to install deal.ii on a cluster.
in my University, that I can run bigger jobs with parallel computations.
I am installing library locally in my home directory.
Before configuring deal.ii, I have installed petsc, p4est-2.0 and
netcdf-cxx-4.2.
Other necessary libraries
Dear Timo,
Sorry, I did not mean to write to your personal e-mail. Just pressed the
wrong reply button by mistake.
Thank you very much for your reply. I will check the way how I set
mpi_communicator now. Hopefully, will find the problem.
I have changed step-40 couple of days ago myself to
Dear Wolfgang,
to compute values of the electric field at the receivers I follow the
strategy of ASPECT code as you suggested
To do this I sum the current_point_values across processors and divide by
the number of processors that contain point p as following
// Reduce all collected values into
Dear Timo,
The main reason I do this is that I do not understand how to reuse this
decomposition in deal.ii.
I am relatively new to deal.ii and C++, and I have never used MUMPS before.
The way I set it up with SparseDirectUMFPACK was to use InnerPreconditioner
structure:
template
struct
Dear Timo,
Because I do not understand how to reuse this decomposition in deal.ii.
I am relatively new to deal.ii and C++, and I have never used mumps before.
I would appreciate any advice on this. Maybe there is some example in deal.ii
that use mumps to construct a preconditioner.
Thank you
Dear Timo,
Thank you for your reply.
I am still having trouble with implementing my code with MUMPS.
I will briefly describe the problem:
I am solving 2 systems of Maxwell's equations in the following way:
the systems are Ax1=b1 and Ax2=b2;
A is a sparse block symmetric matrix (C -M; -M
Dear Timo,
Thank you for your reply.
I am still having trouble with implementing my code with MUMPS.
I will briefly describe the problem:
I am solving 2 systems of Maxwell's equations in the following way:
the systems are Ax1=b1 and Ax2=b2;
A is a sparse block symmetric matrix (C -M; -M
Dear Wolfgang,
thank you very much for the link. I have followed their approach, but
before I can check the result, I have to substitute the solver I used for a
parallel solver.
Anna
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
Dear All,
I have modified step-22 to solve system of Maxwell's equations, and then
used step-40(and 55) to parallelize the code.
In the serial version I used
template
struct InnerPreconditioner;
template<>
struct InnerPreconditioner<0>
{
typedef SparseDirectUMFPACK type;
};
template<>
struct
Dear Wolfgang,
Yes I would like to create a LA::MPI::BlockVector solution_at_receiver. And
then copy it to a
BlockVector localized_solution_at_receiver(solution_at_receiver)
And output localized_solution_at_receiver only on the 0 processor.
I am having trouble initializing and filling
Dear All,
I am not fluent with C++ and MPI and would very much appreciate your help.
I am now following steps 40 and 55 to write parallel implementation of the
solution of the system of Maxwell's equations and am having many questions,
but let us start with the following:
After solving the
Dear All,
I am solving a system of Maxwell's equations for geophysical application.
For my application I need to compute two solutions for two different
polarizations of the source field.
In other words, I am solving the system twice with different
right-hand-sides. And for each polarization, I
Hi All,
I am working on the implementation of 3D electromagnetic forward solver for
geophysical applications.
I am following the algorithm as described in Grayver and Kolev, 2015. The
system is solved with preconditioned FGMRES method. However construction of
the preconditioner is not
Thank you for the prompt reply.
I will implement the first option with VectorTools::point_value/gradient()
.
Could you please also point out to me where I can find clarifications to
the following:
1) There are many versions of VectorTools::point_value function. At the
moment I use
void
Hi All,
I have successfully computed the secondary electric field for 3D case as
described in the paper by Grayver, 2014. These field is declared as
BlockVector solution, where the block correspond to real and imaginary
parts.
Now I would like to compute total electric E (secondary plus
Dear Jean-Paul,
I am new to deal.ii and C++ programming and never used Linear operators, so
will have to study documentation. My understanding so far is that when using
LinearOperators I am not actually constructing the matrix, but only know how
such matrix would act on a vector. Am I correct?
oing to do with the result of the matrix addition?
>
> Regards,
> Jean-Paul
>
> On Wednesday, July 26, 2017 at 5:14:11 AM UTC+2, Anna Avdeeva wrote:
>>
>> Dear All,
>>
>> I need to subtract two blocks of a BlockSparseMatrix
>> system_matrix.
>> I a
I have managed to subtract the two blocks by declaring
SparseMatrixEZ tmp1;
Is there a more efficient way?
--
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
Dear All,
I need to subtract two blocks of a BlockSparseMatrix system_matrix.
I am trying to achieve this by
creating SpraseMatrix tmp1, then
tmp1.copy_from(system_matrix.block(0,0));
tmp1.add(-1.0, system_matrix.block(0,1));
However, I am getting the following error message:
"When copying
24 matches
Mail list logo