Hi there,
Is there anyone working on xfem? I browsed the online tutorial, but
unfortunately couldn't find any implementation of XFEM in deal II. I have
written XFEM, Fast marching, level set codes with Matlab, but not familiar
with deal ii library yet. anyone can give me some hints? I read some
examples recently, but still have some points unclear
1. By justifying the signs of level set function in each element, I can
easily identify the standard element and enriched element. But how can I
change the settings of dof's in these elements. For example, in the
enriched elements, additional dof's exist due to the enrichment. Also, the
local_dof_indices are changed in this process.
2. We can achieve loops over each active cell to calculate local stiffness
matrix using DoFHandler<dim>:: active_cell_iterator. My question is can I
do the loop operation over each nodes and also find the corresponding
supporting elements around it? Any connection functions from nodes to
elements? for example, 3x3 Q1 element, totally 4x4 nodes, I would like to
do the following loops
for i=1:total_num_nodes
support_element_number = some_function( i );
n_support_element = size( support_element_number );
for j=1:n_support_element
... operations in current j-th element
...
end
end
3. If I assembled two sets of global matrix, say, [ A ], [ B ] and rhs {f},
respectively
[ A ]*{u} + [ B ]*{a} = {f}
[ B ]_T *{a} = {g}
how can I combine them as a big matrix K ? for example,
K = [ A B ]
[ BT 0 ]
F = {f}
{g}
Thank you very much.
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
