Hi Peter,
Thank you for your email. Here is what I have so far:
```c++
// This will be filled in by looping over number of orbitals, with each
orbital containing an array of coefs.
struct Orbital_Data{
\\TODO: need a constructor that allocates memory for all the orbital's info
here?
private:
etc., must be in this
project.
I hope this added more clarity.
Thank you,
Zachary
On Monday, November 30, 2020 at 5:46:59 PM UTC-6 Zachary Streeter wrote:
> Hi everyone,
>
> I am wondering if deal.II has any OpenMP capabilities or if they would be
> easily incorporated?
>
> I know
Ah, I see and that makes perfect sense.
I'v been slightly exposed to TaskFlow and it's interesting y'all are
planning on using it. I will look more into this project!
Thank you,
Zachary
On Monday, November 30, 2020 at 7:18:05 PM UTC-6 Wolfgang Bangerth wrote:
> On 11/30/20 6:05 PM, Zach
ely.
>
> Thank you,
> Zachary
>
> On Monday, December 21, 2020 at 11:27:18 AM UTC-6 Wolfgang Bangerth wrote:
>
>> On 12/21/20 10:16 AM, Zachary Streeter wrote:
>> > Sorry, I should say the relevant elements are the same as the locally
>> owned.
>>
>> Th
21, 2020 at 11:27:18 AM UTC-6 Wolfgang Bangerth wrote:
> On 12/21/20 10:16 AM, Zachary Streeter wrote:
> > Sorry, I should say the relevant elements are the same as the locally
> owned.
>
> That's true for cells but not for degrees of freedom. There may be DoFs
> that
>
So it does seem to be the Cmake version. The problem showed up with Cmake
3.19.1 and updating to 3.19.2 fixed it apparently. Wanted to share if
anyone else comes across this.
On Friday, December 18, 2020 at 11:42:41 AM UTC-6 d.arnd...@gmail.com wrote:
> Zachary,
>
> are you using CMake 3.19
y_pattern(), which really just adds certain entries to the
> sparsity pattern like you do, so if you understand how the set up works in
> step-40, you should be able to understand what you need to do in your case.
>
> Best
> W.
>
>
> On 12/20/20 10:21 AM, Zachary Str
s; ++j )
dsp( i, j ) = stuff_calculated
This gives me a local dsp and I thought the IndexSet add_range() above will
put this local dsp to the correct global index set.
On Sunday, December 20, 2020 at 10:53:33 AM UTC-6 Zachary Streeter wrote:
> Hi there,
>
> I am trying to build a parall
Thanks Wolfgang,
I think you are correct! I just tried to use petsc_full_matrix instead and
it linked no problem!
On Monday, November 9, 2020 at 5:25:30 PM UTC-6 Wolfgang Bangerth wrote:
> On 11/9/20 2:49 PM, Daniel Arndt wrote:
> >
> > Undefined symbols for architecture x86_64:
> >
Ah sorry, yes I have done that already. Sorry for forgetting that
detail...
If it should work assuming I have added my new directory in
source/CMakeLists.txt, maybe my implementation is the problem. Though, I
wouldn't know why since it fails at linking.
On Monday, November 9, 2020 at
Ahh, I don't feel bad not thinking that possibility; brilliant!
On Friday, November 13, 2020 at 7:04:14 PM UTC-6 Wolfgang Bangerth wrote:
> On 11/13/20 6:02 PM, Zachary Streeter wrote:
> > Oh the test sees the number of processors from the .output files name?
>
> Yes. 'make setup
Oh the test sees the number of processors from the .output files name?
On Friday, November 13, 2020 at 5:36:21 PM UTC-6 Wolfgang Bangerth wrote:
> On 11/13/20 1:24 PM, Zachary 42! wrote:
> > I am trying to run some MPI tests and I don’t see on the deal.II website
> how to use ctest for running
z/zstreet/dealii/build/share/deal.II/scripts/run_test.cmake:140
(MESSAGE):
1181: *** abort
1181:
1181:
1/2 Test #1181: base/aligned_vector_01.debug .***Failed 13.52 sec
On Wednesday, October 28, 2020 at 5:34:56 PM UTC-5 Zachary Streeter wrote:
> Hi everyone,
>
> I am getting a st
s,
Zachary
On Thursday, October 29, 2020 at 10:55:44 AM UTC-5 d.arnd...@gmail.com
wrote:
> Zachary,
>
> Did you actually build deal.II, i.e did you run make?
>
> Best,
> Daniel
>
> Am Do., 29. Okt. 2020 um 10:39 Uhr schrieb Zachary Streeter <
> zachary...@gmail.com>
, local_columns_per_process, my_proc);
On Tuesday, January 5, 2021 at 12:12:06 PM UTC-6 Wolfgang Bangerth wrote:
> On 1/5/21 8:21 AM, Zachary Streeter wrote:
> > Let me know if this is okay. This compiled, ran, and produced the same
> error
> > on my end.
>
> Yes, that's the sort of tes
#include
#include
#include
#include
#include
using namespace dealii;
class OneBodyHamiltonianOperator
{
public:
/**
* Declare type for container size.
*/
using size_type = dealii::types::global_dof_index;
OneBodyHamiltonianOperator(const dealii::IndexSet _local_row_set,
My project is in quantum scattering and I would like to have some operators
be distributed PETSc objects. So inside my OneBodyHamiltonianOperator
class (for example), I would like to create a
PETScWrappers::MPI::SparseMatrix and then use SLEPC to solve for the ground
state and excited states.
Hello everyone,
I have a matrix A that has a sparsity pattern. I would like to take advantage
of that sparsity pattern and make a PETSc Sparse Matrix in CRS format from it.
My attempt is to go through the matrix A and have an if statement that saves to
a dynamic sparsity pattern when the
complained.)
On Friday, February 5, 2021 at 10:48:00 PM UTC-6 Wolfgang Bangerth wrote:
> On 2/5/21 6:13 PM, Zachary Streeter wrote:
> > Can you point me to the SparseMatrix::copy_from() that takes a
> FullMatrix? I
> > only see a copy_from() that takes in an
?
That's what I have tried but I get the wrong spectrum. If that's what
you're recommending then I will crack on and debug.
Cheers,
Zachary
On Thursday, February 4, 2021 at 11:50:00 AM UTC-6 Wolfgang Bangerth wrote:
> On 2/4/21 9:35 AM, Zachary Streeter wrote:
> > I have a matrix A
:
> On 2/5/21 4:34 PM, Zachary Streeter wrote:
> > I'm not sure what you mean. I create A and it happens to have some
> sparsity
> > pattern, I don't know if the element I save is zero or non-zero during
> > construction. So I put an if statement and if it's non-zero, t
Can you point me to the SparseMatrix::copy_from() that takes a FullMatrix?
I only see a copy_from() that takes in another SparseMatrix in the
PETScWrappers::MPI::SparseMatrix Class reference.
On Friday, February 5, 2021 at 6:52:01 PM UTC-6 Zachary Streeter wrote:
> Oh very cool, thank
Oh very cool, thank you! I'll take a look at that.
Let's say my matrix gets huge, how would you tackle this? I'm expecting it
will, hence why I want CRS.
On Friday, February 5, 2021 at 6:22:32 PM UTC-6 Wolfgang Bangerth wrote:
> On 2/5/21 5:14 PM, Zachary Streeter wrote:
> > O
Oh really? I'm going to check it out right now! Debugging parallel
programs is still a bane of my existence so definitely would like to hear
your take on it!
Cheers,
Zachary
On Friday, February 19, 2021 at 1:29:43 PM UTC-6 Wolfgang Bangerth wrote:
> On 2/19/21 12:23 PM, Zachary Stree
Hi everyone,
I have created a PETSc::MPI::SparseMatrix that has a global matrix distributed
across all the processes. I then pass it to SLEPc to solve for its spectrum.
I get the correct spectrum so I am pretty sure everything is correct. My
question is regarding strange behavior with more
25 matches
Mail list logo