Add step-61, a program that uses the weak Galerkin method #6

sophy1029 · 2018-04-27T16:38:30Z

No description provided.

bangerth · 2018-05-04T17:25:15Z

By the way, if you have time between 2 and 3pm, please come by my office and we can chat about next steps.

sophy1029 · 2018-05-04T17:53:25Z

Yes, I will be there around 2:30pm. I need to go to the measure theory class to get my take home exam first. Thanks.

On Fri, May 4, 2018 at 11:25 AM, Wolfgang Bangerth ***@***.*** > wrote: By the way, if you have time between 2 and 3pm, please come by my office and we can chat about next steps. — You are receiving this because you authored the thread. Reply to this email directly, view it on GitHub <#6 (comment)>, or mute the thread <https://github.com/notifications/unsubscribe-auth/AdMXNwpEVwPS-nDVuRsisO9PzDhn-ue-ks5tvI77gaJpZM4Tqrsh> .

Rename Solver (base) class to SolverBase

Fix a changelog entry.

Allow initialization of Vector with std::initializer_list.

There are no blocks in this object.

Remove '__' from include guards.

Clean up the vector header a bit.

…_number Update PETSc version compatibility

Use a range-based for loop in a couple of places.

CI: archive detailed.log

several doxygen errors and warnings

Unify intrinsics headers

Restrict SymEngine test

Clean up CUDAWrappers::MatrixFree::reinit interface

Arpack Solver documentation fix

…file

…_fix newline bug in grid_in_abaqus fixed

Update module on user interfaces to AD libraries

Edit the results section of step-61.

Enable doxygen to build SD documentation

Fix a couple of typos.

Align a set of equations.

Mapping::get_center()

Step-63

SD scalar operations: Creation of substitution maps

Fixes: ``` 134: ==1438450==ERROR: AddressSanitizer: stack-use-after-return on address 0x7fa4f4d4bfb0 at pc 0x56248d8610c2 bp 0x7ffeebad19b0 sp 0x7ffeebad19a8 134: READ of size 1 at 0x7fa4f4d4bfb0 thread T0 134: #0 0x56248d8610c1 in ExponentialDecay::ExponentialDecay(double, dealii::PETScWrappers::TimeStepperData const&, bool, bool, int)::'lambda'(double, dealii::[148/1839] ers::MPI::Vector const&, dealii::PETScWrappers::MPI::Vector const&, double, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&)::operator()(double, deali i::PETScWrappers::MPI::Vector const&, dealii::PETScWrappers::MPI::Vector const&, double, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&) const /srv/t estsuite/dealii/tests/petsc/petsc_ts_05.cc:91:15 134: #1 0x7fa538cf812b in std::__1::__function::__value_func<void (double, dealii::PETScWrappers::MPI::Vector const&, dealii::PETScWrappers::MPI::Vector const&, double, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&)>::operator()[abi:v160006](double&&, dealii::PETScWrappers::MPI::Vector const&, dealii::PETScWrappers: :MPI::Vector const&, double&&, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&) const /usr/include/c++/v1/__functional/function.h:510:16 134: #2 0x7fa538cf812b in std::__1::function<void (double, dealii::PETScWrappers::MPI::Vector const&, dealii::PETScWrappers::MPI::Vector const&, double, dealii::PETScWra ppers::MatrixBase&, dealii::PETScWrappers::MatrixBase&)>::operator()(double, dealii::PETScWrappers::MPI::Vector const&, dealii::PETScWrappers::MPI::Vector const&, double, de alii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&) const /usr/include/c++/v1/__functional/function.h:1156:12 134: #3 0x7fa538cf812b in int dealii::PETScWrappers::call_and_possibly_capture_ts_exception<std::__1::function<void (double, dealii::PETScWrappers::MPI::Vector const&, d ealii::PETScWrappers::MPI::Vector const&, double, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&)>, double&, dealii::PETScWrappers::MPI::Vector&, dea lii::PETScWrappers::MPI::Vector&, double&, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&>(std::__1::function<void (double, dealii::PETScWrappers::MP I::Vector const&, dealii::PETScWrappers::MPI::Vector const&, double, dealii::PETScWrappers::MatrixBase&, dealii::PETScWrappers::MatrixBase&)> const&, std::exception_ptr&, st d::__1::function<void ()> const&, double&, dealii::PETScWrappers::MPI::Vector&, dealii::PETScWrappers::MPI::Vector&, double&, dealii::PETScWrappers::MatrixBase&, dealii::PET ScWrappers::MatrixBase&) /srv/testsuite/dealii/include/deal.II/lac/petsc_ts.templates.h:87:9 134: #4 0x7fa538cf78cc in dealii::PETScWrappers::TimeStepper<dealii::PETScWrappers::MPI::Vector, dealii::PETScWrappers::MatrixBase, dealii::PETScWrappers::MatrixBase>::s olve(dealii::PETScWrappers::MPI::Vector&)::'lambda'(_p_TS*, double, _p_Vec*, _p_Vec*, double, _p_Mat*, _p_Mat*, void*)::operator()(_p_TS*, double, _p_Vec*, _p_Vec*, double, _p_Mat*, _p_Mat*, void*) const /srv/testsuite/dealii/include/deal.II/lac/petsc_ts.templates.h:489:26 134: #5 0x7fa538cf73d7 in dealii::PETScWrappers::TimeStepper<dealii::PETScWrappers::MPI::Vector, dealii::PETScWrappers::MatrixBase, dealii::PETScWrappers::MatrixBase>::s olve(dealii::PETScWrappers::MPI::Vector&)::'lambda'(_p_TS*, double, _p_Vec*, _p_Vec*, double, _p_Mat*, _p_Mat*, void*)::__invoke(_p_TS*, double, _p_Vec*, _p_Vec*, double, _p _Mat*, _p_Mat*, void*) /srv/testsuite/dealii/include/deal.II/lac/petsc_ts.templates.h:472:31 134: #6 0x7fa504afbcac in TSComputeIJacobian (/usr/lib64/libpetsc.so.3.19+0x10fbcac) 134: #7 0x7fa504aa1672 (/usr/lib64/libpetsc.so.3.19+0x10a1672) 134: #8 0x7fa504afd307 in SNESTSFormJacobian (/usr/lib64/libpetsc.so.3.19+0x10fd307) 134: #9 0x7fa504a43e66 in SNESComputeJacobian (/usr/lib64/libpetsc.so.3.19+0x1043e66) 134: #10 0x7fa504a35e98 in SNESSolve_NEWTONLS (/usr/lib64/libpetsc.so.3.19+0x1035e98) 134: #11 0x7fa504a4968e in SNESSolve (/usr/lib64/libpetsc.so.3.19+0x104968e) 134: #12 0x7fa504aa25e4 (/usr/lib64/libpetsc.so.3.19+0x10a25e4) 134: #13 0x7fa504aa014a (/usr/lib64/libpetsc.so.3.19+0x10a014a) 134: #14 0x7fa504b03013 in TSStep (/usr/lib64/libpetsc.so.3.19+0x1103013) 134: #15 0x7fa504b03e02 in TSSolve (/usr/lib64/libpetsc.so.3.19+0x1103e02) 134: #16 0x7fa538cd2ccc in dealii::PETScWrappers::TimeStepper<dealii::PETScWrappers::MPI::Vector, dealii::PETScWrappers::MatrixBase, dealii::PETScWrappers::MatrixBase>:: solve(dealii::PETScWrappers::MPI::Vector&) /srv/testsuite/dealii/include/deal.II/lac/petsc_ts.templates.h:880:24 134: #17 0x56248d85849f in ExponentialDecay::run() /srv/testsuite/dealii/tests/petsc/petsc_ts_05.cc:182:32 134: #18 0x56248d855ca6 in main /srv/testsuite/dealii/tests/petsc/petsc_ts_05.cc:242:20 134: #19 0x7fa502a50989 (/usr/lib64/libc.so.6+0x23989) 134: #20 0x7fa502a50a44 in __libc_start_main (/usr/lib64/libc.so.6+0x23a44) 134: #21 0x56248d73e2d0 in _start (/srv/temp/build/tests/petsc/petsc_ts_05.debug/petsc_ts_05.debug+0x742d0) ```

``` ==1080297==ERROR: AddressSanitizer: attempting free on address which was not malloc()-ed: 0x60200003da70 in thread T0 #0 0x55bcdd907b7d in operator delete(void*) (/srv/temp/testsuite-IQZ1b8kK/build/tests/scalapack/scalapack_06b.debug/scalapack_06b.debug+0x17cb7d) #1 0x7fc52a4d2047 in void std::__1::__libcpp_operator_delete[abi:v160006]<void*>(void*) /usr/include/c++/v1/new:276:3 #2 0x7fc52a4d2047 in void std::__1::__do_deallocate_handle_size[abi:v160006]<>(void*, unsigned long) /usr/include/c++/v1/new:300:10 #3 0x7fc52a4d2047 in std::__1::__libcpp_deallocate[abi:v160006](void*, unsigned long, unsigned long) /usr/include/c++/v1/new:316:14 #4 0x7fc52a4d2047 in std::__1::allocator<int>::deallocate[abi:v160006](int*, unsigned long) /usr/include/c++/v1/__memory/allocator.h:131:13 #5 0x7fc52a4d2047 in std::__1::allocator_traits<std::__1::allocator<int>>::deallocate[abi:v160006](std::__1::allocator<int>&, int*, unsigned long) /usr/include/c++/v1/__memory/allocator_traits.h:288:13 #6 0x7fc52a4d2047 in std::__1::__split_buffer<int, std::__1::allocator<int>&>::~__split_buffer() /usr/include/c++/v1/__split_buffer:362:9 #7 0x7fc52a4d2047 in std::__1::vector<int, std::__1::allocator<int>>::__append(unsigned long) /usr/include/c++/v1/vector:1049:5 #8 0x7fc52a8282d2 in std::__1::vector<int, std::__1::allocator<int>>::resize(unsigned long) /usr/include/c++/v1/vector:1910:15 #9 0x7fc52a8282d2 in dealii::ScaLAPACKMatrix<double>::eigenpairs_symmetric(bool, std::__1::pair<unsigned int, unsigned int> const&, std::__1::pair<double, double> const&) /srv/temp/testsuite-IQZ1b8kK/dealii/source/lac/scalapack.cc:1684:17 #10 0x7fc52a826813 in dealii::ScaLAPACKMatrix<double>::eigenpairs_symmetric_by_index(std::__1::pair<unsigned int, unsigned int> const&, bool) /srv/temp/testsuite-IQZ1b8kK/dealii/source/lac/scalapack.cc:1446:12 #11 0x55bcdd911907 in void test<double>(unsigned int, unsigned int, double) /srv/temp/testsuite-IQZ1b8kK/dealii/tests/scalapack/scalapack_06b.cc:134:21 #12 0x55bcdd90ebc8 in main /srv/temp/testsuite-IQZ1b8kK/dealii/tests/scalapack/scalapack_06b.cc:207:9 #13 0x7fc4f5250989 (/usr/lib64/libc.so.6+0x23989) #14 0x7fc4f5250a44 in __libc_start_main (/usr/lib64/libc.so.6+0x23a44) #15 0x55bcdd7f73f0 in _start (/srv/temp/testsuite-IQZ1b8kK/build/tests/scalapack/scalapack_06b.debug/scalapack_06b.debug+0x6c3f0) ```

sophy1029 force-pushed the step-61 branch from 6156ce3 to 12eebbb Compare August 29, 2018 20:03

sophy1029 force-pushed the step-61 branch 3 times, most recently from 4d201ee to ad3df66 Compare January 30, 2019 17:13

bangerth and others added 24 commits March 28, 2019 09:00

Allow initialization of Vector with std::initializer_list.

9680c30

Add changelog.

f0b996c

Add a test.

292c1c9

Fix a changelog entry.

cebc4fc

Merge pull request dealii#7864 from jppelteret/rename_solver_class

ab39c1b

Rename Solver (base) class to SolverBase

Merge pull request dealii#7872 from bangerth/fix-changelog-entry

0e60f73

Fix a changelog entry.

Merge pull request dealii#7871 from bangerth/initializer_list

0df3420

Allow initialization of Vector with std::initializer_list.

Use a range-based for loop.

cb09939

Fix a comment.

7f12c2c

There are no blocks in this object.

Clean up some headers and code samples.

330bd8c

Remove '__' from include guards.

c9539d4

Merge pull request dealii#7874 from drwells/underscore-guards

ffb94e8

Remove '__' from include guards.

Merge pull request dealii#7873 from drwells/vector-header

0c65f3a

Clean up the vector header a bit.

Update PETSc version compatibility

8d35766

CI: archive detailed.log

3ba0690

fix path

cd518ac

Merge pull request dealii#7877 from masterleinad/update_petsc_version…

85bfd86

…_number Update PETSc version compatibility

Merge pull request dealii#7878 from bangerth/for

704f110

Use a range-based for loop in a couple of places.

Merge pull request dealii#7875 from tjhei/jenkins_archive_detailedlog

2a76046

CI: archive detailed.log

several doxygen errors and warnings

ec8ee20

Unify intrinsics headers

3a1f07a

Merge pull request dealii#7879 from tjhei/doxygen_fixes2

e309199

several doxygen errors and warnings

Merge pull request dealii#7881 from masterleinad/unify_intrinsic_headers

b68e137

Unify intrinsics headers

set LAPACK_INCLUDE_DIRS and LAPACK_USER_INCLUDE_DIRS

fb6d449

jppelteret and others added 26 commits April 30, 2019 09:35

Merge pull request dealii#7981 from masterleinad/restrict_symengine_test

7a7bbd0

Restrict SymEngine test

Arpack Solver documentation fix

1a78ff6

newline bug in grid_in_abaqus fixed

f844461

grid_in_abaqus test added

2d04811

Merge pull request dealii#7966 from masterleinad/cleanup_cuda_mf_reinit

8ff21b1

Clean up CUDAWrappers::MatrixFree::reinit interface

Merge pull request dealii#7982 from rezarastak/patch-1

4aff835

Arpack Solver documentation fix

fix typo in 2d_short_handwritten.inp

5e7b354

Edit the results section of step-61.

98a2c16

Use 'align' instead of 'eqnarray'.

d956724

reverting the effects of the last commit on the last line of the inp …

a780a7d

…file

Build documentation with SymEngine enabled

e039666

Add placeholder for documentation on symbolic differentiation

f15930f

Update module on user interfaces to AD libraries

e09d3b4

Merge pull request dealii#7983 from rezarastak/grid_in_abaqus_newline…

b6f4b71

…_fix newline bug in grid_in_abaqus fixed

Merge pull request dealii#7985 from jppelteret/ad-doc_02

e1d8afb

Update module on user interfaces to AD libraries

Merge pull request dealii#7954 from bangerth/61-4

09eb6bf

Edit the results section of step-61.

Merge pull request dealii#7984 from jppelteret/sd-doc_01

c9b35c9

Enable doxygen to build SD documentation

Align a set of equations.

3cd2401

Fix a couple of typos.

44a7a39

Merge pull request dealii#7987 from bangerth/fix-typos

2035b7f

Fix a couple of typos.

Merge pull request dealii#7986 from bangerth/61-5

92af75a

Align a set of equations.

Merge pull request dealii#7913 from luca-heltai/mapping-center

89057df

Mapping::get_center()

Merge pull request dealii#7738 from tcclevenger/step-63

406a15d

Step-63

Merge pull request dealii#7974 from jppelteret/sd-map_utilities_03

7d2c55c

SD scalar operations: Creation of substitution maps

step-61 announcement.

f12d8b1

Minor update to a changelog entry.

2b35c1d

bangerth merged commit 2b35c1d into bangerth:master May 13, 2019

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Add step-61, a program that uses the weak Galerkin method #6

Add step-61, a program that uses the weak Galerkin method #6

sophy1029 commented Apr 27, 2018

bangerth commented May 4, 2018

sophy1029 commented May 4, 2018 via email

Add step-61, a program that uses the weak Galerkin method #6

Add step-61, a program that uses the weak Galerkin method #6

Conversation

sophy1029 commented Apr 27, 2018

bangerth commented May 4, 2018

sophy1029 commented May 4, 2018 via email