Re: [deal.II] Use Gmsh to import external grid calculation elasticity problem
Okay, I re-run the program I wrote vs told me that ,what that mean? *2>step-8.obj: error LNK2019: Unresolved external symbol "public: void __cdecl Step8::ElasticProblem<2>::run(void)" (?run@?$ElasticProblem@$01@Step8@@QEAAXXZ) , The symbol is quoted in the function main 2>D:\dealii-9.2.0\dealii-9.2.0\examples\step-8a\Debug\step-8.exe: fatal error LNK1120: 1 unresolved external command* 在2020年11月4日星期三 UTC+8 上午11:09:17 写道： > On 11/3/20 7:33 PM, Nick Wang wrote: > > > > After I studied and understood the tutorials of step-8 and step-49. I > tried > > to use step-49 to import external grids to to calculate the force > condition > > and stress analysis, > > So I modified the statement of step8 grid generation, > > while failed... > > Did I make a mistake or are there other ways to finish the force > analysis > > from external grids? > > Help us understand what the error message is that you get. You already > have it > on your screen, but we would need to download your program, put it > somewhere, > compile it, maybe run it. It would take someone a couple of minutes to do > that > -- time we wouldn't have for answering someone else's problems. > > Best > WB > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- 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 you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/4fa565f4-1ab1-4919-8818-3b58fe6169c4n%40googlegroups.com.
Re: [deal.II] Running new tests with argument
Hi Zachary, Are you adding a test to deal.II itself? If so all tests are compiled with the macro SOURCE_DIR defined as the absolute path to the test's source directory. This is how we load grids in, e.g., tests/grid/grid_in.cc. If you are adding tests using our test framework but not deal.II: I am pretty sure that this macro is defined in the same way. Best, David From: dealii@googlegroups.com on behalf of Zachary 42! Sent: Wednesday, November 4, 2020 6:38 PM To: dealii@googlegroups.com Subject: [deal.II] Running new tests with argument 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 mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/F4D11B95-0211-40DD-A97C-326B6A358261%40gmail.com. -- 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 you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/BN6PR03MB3345FEEE06EDB1D15F90A344EDEE0%40BN6PR03MB3345.namprd03.prod.outlook.com.
[deal.II] Running new tests with argument
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 mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/F4D11B95-0211-40DD-A97C-326B6A358261%40gmail.com.
[deal.II] Post-processing velocities with Trilinos vectors during simulation
Hi everyone!
I'm working on post-processing velocities with Trilinos solution vectors
during the simulation on Lethe.
Calculating average velocities and pressures (, , and ) works
well using Trilinos vectors with no ghost cells and .equ(...) function.
The calculation of average Reynolds stress (, and and
average shear stress (), where u' = u - , is quite not easy.
It seems that I can't do what I want in parallel without working with loop.
The problem is that loops seem to take way too much time on Trilion vectors
and the local_range() function doesn't work with BlockVector. (Am I wrong?)
I have tried some ways to do it.
1. Trying working with Deal.ii vectors : did not work in parallel
2. Doing summations and/or multiplications of the vectors
local_evaluation_point to get u'u', v'v' and w'w' and replacing the fourth
element of the evaluation point with a loop for u'v' : too much time.
3. Doing loops on the Trilinos vectors : too much time
Does anyone can give me an advice to find a way to get one or two vectors
with my precious time-averaged Reynolds stress and time-averaged shear
stress?
I did not mention the exact way I'm trying to calculate the average of all
of that, but I'm posting my code where I'm trying to do a loop on the
solution vector.
The variables average_velocities, total_time and inv_range_time are
calculated in an other function is my class AverageVelocities.
template
VectorType
AverageVelocities::calculate_reynolds_stress(
const VectorType &local_evaluation_point,
const std::shared_ptr &simulation_control,
const DofsType & locally_owned_dofs,
const MPI_Comm & mpi_communicator)
{
if (simulation_control->get_step_number() == 0)
{
// Reinitializing vectors with zeros at t = 0
sum_reynolds_stress_dt.reinit(locally_owned_dofs,
mpi_communicator);
reynolds_stress.reinit(locally_owned_dofs,
mpi_communicator);
}
else if (abs(total_time) < 1e-6 || total_time > 1e-6)
{
VectorType reynolds_stress_dt(locally_owned_dofs,
mpi_communicator);
// ***Won't work with BlockVectors
if constexpr (std::is_same_v)
{
for (unsigned int i = local_evaluation_point.local_range().first;
i < local_evaluation_point.local_range().second; i++)
{
if ((i + 4) % 4 == 0)
{
// Calculating (u'u')*dt, (v'v')*dt (w'w')*dt and (u'v')*dt
reynolds_stress_dt[i] =
(local_evaluation_point[i] - average_velocities[i]) *
(local_evaluation_point[i] - average_velocities[i]) * dt;
reynolds_stress_dt[i + 1] =
(local_evaluation_point[i + 1] - average_velocities[i + 1]) *
(local_evaluation_point[i + 1] - average_velocities[i + 1]) *
dt;
reynolds_stress_dt[i + 2] =
(local_evaluation_point[i + 2] - average_velocities[i + 2]) *
(local_evaluation_point[i + 2] - average_velocities[i + 2]) *
dt;
reynolds_stress_dt[i + 3] =
(local_evaluation_point[i] - average_velocities[i]) *
(local_evaluation_point[i + 1] - average_velocities[i + 1]) *
dt;
// Summation of all reynolds stress during simulation
sum_reynolds_stress_dt += reynolds_stress_dt;
// Calculating time-averaged reynolds stress if output needed
if (simulation_control->is_output_iteration())
reynolds_stress.equ(inv_range_time, sum_reynolds_stress_dt);
}
}
}
}
return reynolds_stress;
}
Thank you in advance!
--
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 you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/7939af05-e474-4b4e-bf5f-8d6a9b4786c1n%40googlegroups.com.
Re: [deal.II] Boundary DOFs using DOF Handler
Thank you, Professor Bangerth. This is really helpful. Animesh On Wednesday, November 4, 2020 at 11:45:46 PM UTC+5:30 Wolfgang Bangerth wrote: > On 11/4/20 12:20 AM, Animesh Rastogi IIT Gandhinagar wrote: > > > > I have a quick question about the dof_handler in dealii. I have a > rectangular > > geometry and would like to access the numbering of the degree of > freedoms > > which are on the right hand side boundary of the rectangle (+X face). I > know > > the boundary_id for that face. Essentially, I want to use that number to > > extract the x component of the displacement at that node using the > > displacement vector. > > You can get all DoF indices on a particular boundary via > DoFTools::extract_boundary_dofs(). You can correlate this > DoFTools::extract_dofs() with a component mask that corresponds to your > x-displacement. With the information produced by these two you should be > able > to pick out which DoFs you want to work on. > > Best > W. > > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > www: http://www.math.colostate.edu/~bangerth/ > > -- 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 you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/e3e30b1c-9d88-452d-b9fb-a1f6cc3cd697n%40googlegroups.com.
Re: [deal.II] Boundary DOFs using DOF Handler
On 11/4/20 12:20 AM, Animesh Rastogi IIT Gandhinagar wrote: I have a quick question about the dof_handler in dealii. I have a rectangular geometry and would like to access the numbering of the degree of freedoms which are on the right hand side boundary of the rectangle (+X face). I know the boundary_id for that face. Essentially, I want to use that number to extract the x component of the displacement at that node using the displacement vector. You can get all DoF indices on a particular boundary via DoFTools::extract_boundary_dofs(). You can correlate this DoFTools::extract_dofs() with a component mask that corresponds to your x-displacement. With the information produced by these two you should be able to pick out which DoFs you want to work on. Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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 you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/dbf0b8de-1a64-49ee-b324-528a7330cc58%40colostate.edu.
[deal.II] deal.II Newsletter #141
Hello everyone! This is deal.II newsletter #141. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #11144: Label functions as such (and not as classes). (proposed by drwells) https://github.com/dealii/dealii/pull/11144 #11143: step-7: mark some ints as doubles (proposed by tjhei) https://github.com/dealii/dealii/pull/11143 #11141: Add internal::TriangulationImplementation::ImplementationMixedMesh (proposed by peterrum) https://github.com/dealii/dealii/pull/11141 #11140: Store reference_cell_types in Triangulation (proposed by peterrum) https://github.com/dealii/dealii/pull/11140 #11139: Avoid the use of 'illegal' in error messages. (proposed by bangerth; merged) https://github.com/dealii/dealii/pull/11139 #11137: Introduce Utilities::MPI::mpi_processes_within_communicator (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11137 #11136: Use global_active_cell_index in MatrixFree::initialize_indices (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11136 #11135: Extend docu on global_active_cell_index (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11135 #11134: Expose mesh loop bug for anisotropic grids. (proposed by luca-heltai) https://github.com/dealii/dealii/pull/11134 #11133: DoFHandler housekeeping. (proposed by marcfehling) https://github.com/dealii/dealii/pull/11133 #11132: [WIP] Enable element activation distributed::SolutionTransfer (proposed by Rombur) https://github.com/dealii/dealii/pull/11132 #11131: Transfinite interpolation: add a test (proposed by tamiko) https://github.com/dealii/dealii/pull/11131 #11130: set number_cache.n_global_dofs to the right size in MultiGrid scenario (proposed by mathsen; merged) https://github.com/dealii/dealii/pull/11130 #11129: fix type in Physics rotation_matrix_2d and 3d (proposed by tjhei; merged) https://github.com/dealii/dealii/pull/11129 #11128: ADHelpers: Fix indexing for symm tensor extractor (proposed by jppelteret) https://github.com/dealii/dealii/pull/11128 #11127: Transfinite interpolation: inverse quadratic approximation in pull-back (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/11127 #11126: Transfinite interpolation: Fix refresh of Jacobian (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/11126 #11125: Update ghost particles without rebuilding them (proposed by blaisb; merged) https://github.com/dealii/dealii/pull/11125 #11123: Hide hp_capability_enabled parameter. (proposed by marcfehling; merged) https://github.com/dealii/dealii/pull/11123 #11122: Check MSVC's constexpr support. (proposed by drwells) https://github.com/dealii/dealii/pull/11122 #11121: Implement face_to_cell_vertices for pyramids and wedges. (proposed by drwells; merged) https://github.com/dealii/dealii/pull/11121 #9: update to clang-format-11 restart with fresh .clang-format (proposed by tamiko) https://github.com/dealii/dealii/pull/9 #8: Fix Wedge::n_faces(). (proposed by drwells; merged) https://github.com/dealii/dealii/pull/8 #7: Fix Wedge::n_faces(). (proposed by drwells) https://github.com/dealii/dealii/pull/7 #6: Fix multiple definitions for CUDAWrappes::internal::used_objects (proposed by masterleinad; merged) https://github.com/dealii/dealii/pull/6 #4: Switch SymmetricTensor -> Tensor for better performance (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/4 #3: Do not capture by reference in range-based for over IndexSet (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/3 #11104: Exchange ghost particles optimization by not rebuilding vertex_to_neighbor_subdomain (proposed by blaisb; merged) https://github.com/dealii/dealii/pull/11104 #11103: Pass n_q_points_1D to Simplex::QGauss (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11103 #11098: Introduce VectorDataExchange::Full (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11098 #11092: MatrixFree: add DoFInfo::dof_indices_contiguous_sm (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11092 #11090: Extend GridOut::write_svg() for simplices (proposed by eliasstudiert; merged) https://github.com/dealii/dealii/pull/11090 #11085: Inverse quadratic approximation for MappingQ::real_to_unit_cell (proposed by kronbichler; merged) https://github.com/dealii/dealii/pull/11085 #11074: Introduce MPI-3 shared-memory capabilities in L:d:V (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/11074 #10490: Give each cell a globally unique (level) ID (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/10490 #10439: Introduce internal::TriangulationImplementation::Policy (proposed by peterrum; merged) https://github.com/dealii/dealii/pull/10439 ## And
