[deal.II] PETSc Vector set deprecated

2021-03-28 Thread Zachary 42!
Hello, I am building PETSc vectors with the set operation and I see it’s deprecated. The documentation suggests to use “import( )” but I don’t see how to use it. How can I use import( ) as a replacement for set( ) as the documentation suggests? Please provide an example. Thank you,

Re: [deal.II] PETSC SparsMatrix initialization error

2021-01-08 Thread Zachary 42!
ccurred while calling a SLEPc function -Zachary > On Jan 8, 2021, at 5:25 PM, Wolfgang Bangerth wrote: > > On 1/8/21 3:58 PM, Zachary 42! wrote: >> Yes, the error comes when I try to use the Slepc solver (the matrix builds >> without errors). Running the att

Re: [deal.II] PETSC SparsMatrix initialization error

2021-01-08 Thread Zachary 42!
No, that's not what you want. You need to partition the columns in the same way as you partition the vectors you want to multiply the matrix by. We generally partition the vectors in the same way as the rows of the matrix, and so then you also have to choose the same partitioning for the columns

Re: [deal.II] PETSC SparsMatrix initialization error

2021-01-08 Thread Zachary 42!
_reinit.cc Description: Binary data > On Jan 7, 2021, at 10:43 PM, Wolfgang Bangerth wrote: > > On 1/6/21 8:07 AM, Zachary 42! wrote: >> Here is my code again with hopefully better comments. Looking at the loop >> structure for the dynamic_sparsity_pattern should make it clear

[deal.II] Re: Slepc Solver Dimension error

2021-01-07 Thread Zachary 42!
> On Jan 7, 2021, at 7:15 PM, Zachary 42! wrote: > > Hi everyone, > > I get a “sum of local lengths 32 does not equal global length 128” PETSC > error in the following test file. It dies at the eigensolver.solve step. I > tried to use the test and tutorial programs but to

[deal.II] Slepc Solver Dimension error

2021-01-07 Thread Zachary 42!
Hi everyone, I get a “sum of local lengths 32 does not equal global length 128” PETSC error in the following test file. It dies at the eigensolver.solve step. I tried to use the test and tutorial programs but to no avail. Still wrapping my head around the logic so please excuse me. Thank

[deal.II] AlignedVector set at index

2021-01-06 Thread Zachary 42!
Hi everyone, Ideally I would like to create an AlignedVector with a begin and end and then push_back only to those indices within those bounds. Is there a way to accomplish this? I assume the push_back function works similar to std::vector and starts the index count at zero, I would like to

Re: [deal.II] PETSC SparsMatrix initialization error

2021-01-06 Thread Zachary 42!
Sorry, I am not sure why my code’s format changed when I copied and pasted, trying again... #include #include #include #include #include using namespace dealii; class OneBodyHamiltonianOperator { public: /** * Declare type for container size. */ using size_type =

Re: [deal.II] PETSC SparsMatrix initialization error

2021-01-06 Thread Zachary 42!
W., Ah yes you are correct in your thinking. My test case was every process did have the same number of rows but I would need to communicate that with a gather to account for the number of rows not divisible by the number of given processors. I will use the 4th reinit function to avoid

[deal.II] PETSC SparsMatrix initialization error

2021-01-04 Thread Zachary 42!
Hi everyone, I am trying to debug this strange behavior. I am trying to build a PETSC sparse parallel matrix using 4 processors. This gives me 32 local number of rows (so 128 global number of rows). But when I pass the local_num_of_rows variable into the reinit function, this is the PETSC

[deal.II] Distribute local sparsity pattern

2020-12-20 Thread Zachary 42!
Hi there, I am trying to build a parallel PETSc sparse matrix by building a local sparsity pattern and subsequently distributing the pattern to all processes so I then have the full sparsity pattern for the PETSc sparse matrix (I read this needs to be done but please correct me if I am wrong).

[deal.II] New Build Error? Usable debug flags failed

2020-12-18 Thread Zachary 42!
Hello, I recently tried to build deal.II locally because the cluster is down and I have the error below. How can I attack this? I appreciate any guidance in advance! -- Include /Users/zachary/Documents/Coding/dealii/cmake/setup_finalize.cmake -- Performing Test

[deal.II] Calculating local PETSC matrix

2020-12-08 Thread Zachary 42!
Hi everyone, I want to create a PETSC matrix across multiple nodes using MPI. This should be simply finding the appropriate local index within my loops and build up the matrix. What are some of the functionalities I should look into for this? There are multiple “partition” so I thought it

[deal.II] OpenMP

2020-11-30 Thread Zachary 42!
Hi everyone, I am wondering if deal.II has any OpenMP capabilities or if they would be easily incorporated? I know the deal.II project uses TBB for multithreading but I would like to run multithreaded on computers that may not have TBB. Cheers, Zachary -- The deal.II project is located at

[deal.II] Creating a SOA with nice OOP like interface

2020-11-21 Thread Zachary 42!
Hi everyone, I saw on lecture 22 “Some data structure design considerations” Prof. Wolfgang described a struct of arrays (SOA) with an accessor interface. This keeps nice encapsulation while also providing nice cache use. I would like to understand this design better for my personal

[deal.II] Compiling f90 files

2020-11-20 Thread Zachary 42!
Hi everyone, I want to use some fortran modules in my test build of deal.II. What would be the best way to create this structure? Maybe create a fortran directory in source and create a obj_fortran OBJECT but I am not sure if the DEAL_II_ADD_LIBRARY macro will work with .f90 files. This

[deal.II] Building on Cluster

2020-11-17 Thread Zachary 42!
Hi everyone, I am having an issue with building deal.II on my groups cluster. I am trying to build with Intel MPI compilers like mpiicc and mpiicpc. Here is the output where the error is: - Include /global/home/users/lstreeter/dealii/cmake/checks/check_01_cxx_features.cmake -- Performing

[deal.II] CTest with MPI

2020-11-13 Thread Zachary 42!
Hi everyone, I am trying to run some MPI tests and I don’t see on the deal.II website how to use ctest for running MPI programs. I just saw how the output file needs to indicate the number of processors used in the test so numdiff works correctly. How do you use ctest for running a MPI test

[deal.II] Adding new obj target

2020-11-09 Thread Zachary 42!
Hi folks, I am trying to add another include/source directory and subsequent obj to dealii for potential future contributions to dealii. Though I am having linking errors because my new obj needs to link agains the lac_obj (need Trilinos and other interfaces). I thought the ADD_DEPENDICIES

[deal.II] Running new tests with argument

2020-11-04 Thread Zachary 42!
Hi folks, I am trying to add a new test that requires reading in a file so I’m passing it in as an argument. I am not sure how to use CTest and the other test setup in deal.II to pass in a file as an argument. Cheers, Zachary -- The deal.II project is located at http://www.dealii.org/ For