Hello Evan,
is the h-convergence as it is supposed to be?
Best Regards,
Markus
Am 22.11.10 05:52, schrieb pleramorphyllactic:
Hi all,
I'm having an issue with my p-convergence. Basically, if I use the
native mesh generator I have no problems, using:
GridGenerator::hyper_cube(triangulation, -5, 5);
triangulation.refine_global(n_refinements);
However, when I try to input a very simple square mesh
std::ifstream input2_file("square.inp");
GridIn<dim> grid_in;
grid_in.attach_triangulation(triangulation);
grid_in.read_ucd (input2_file);
I lose the p-convergence of my solutions; which seems to indicate an
orientation problem to me (but the mesh "looks fine" in paraview).
Does anybody have a suggestion? It seems that read.ucd should already
be creating triangulation compatibility, so I am a little lost at the
moment.
Thanks again,
Evan
Here's the square.inp file
121 140 0 0 0
1 -5.0 5.0 0
2 -5.0 -5.0 0
3 -5.0 4.0 0
4 -5.0 3.0 0
5 -5.0 2.0 0
6 -5.0 1.0 0
7 -5.0 0 0
8 -5.0 -1.0 0
9 -5.0 -2.0 0
10 -5.0 -3.0 0
11 -5.0 -4.0 0
12 5.0 -5.0 0
13 -4.0 -5.0 0
14 -3.0 -5.0 0
15 -2.0 -5.0 0
16 -1.0 -5.0 0
17 0 -5.0 0
18 1.0 -5.0 0
19 2.0 -5.0 0
20 3.0 -5.0 0
21 4.0 -5.0 0
22 5.0 5.0 0
23 5.0 -4.0 0
24 5.0 -3.0 0
25 5.0 -2.0 0
26 5.0 -1.0 0
27 5.0 0 0
28 5.0 1.0 0
29 5.0 2.0 0
30 5.0 3.0 0
31 5.0 4.0 0
32 4.0 5.0 0
33 3.0 5.0 0
34 2.0 5.0 0
35 1.0 5.0 0
36 0 5.0 0
37 -1.0 5.0 0
38 -2.0 5.0 0
39 -3.0 5.0 0
40 -4.0 5.0 0
41 -3.99975235074 3.99924883023 0
42 -3.99989359546 2.9995966829 0
43 -4.00013521822 1.99962042975 0
44 -4.0001508689 0.999285214472 0
45 -4.00041147345 -0.000304069995063 0
46 -3.99980818189 -1.00029449897 0
47 -3.99996693887 -1.9998312425 0
48 -4.000040745 -2.99988833531 0
49 -3.99987043447 -3.99980731004 0
50 3.99977526303 -3.99925880274 0
51 3.99995289676 -2.99961415031 0
52 4.00025861047 -1.99964076392 0
53 4.00035979732 -0.999312624167 0
54 4.00055829848 0.000269678516263 0
55 3.99998703832 1.00025187577 0
56 4.00002489984 1.99980329952 0
57 4.00005681287 2.9998805897 0
58 3.99987462604 3.99980531835 0
59 -3.00008463094 -3.99983034595 0
60 -1.9999232651 -4.00000198665 0
61 -0.999394142418 -4.00014448777 0
62 0.0018794941171 -4.00052150483 0
63 1.00104717485 -4.00061075955 0
64 2.00017748184 -3.99980543411 0
65 2.99987719957 -3.99967876582 0
66 3.0000869903 3.99982928682 0
67 1.99988435329 4.00000095343 0
68 0.999439151929 3.99996852278 0
69 -0.00178415570107 4.00040405921 0
70 -1.00093949378 4.00050637683 0
71 -2.00008548861 3.99974038843 0
72 -2.99981517431 3.99964736308 0
73 -2.99984276228 -3.00013671807 0
74 -1.99941052588 -3.00011339788 0
75 -0.998697122202 -3.00028756076 0
76 0.00166110513774 -3.00047941055 0
77 1.00101614937 -3.00035234917 0
78 2.00029833069 -2.99990143369 0
79 3.00022212828 -2.99969706146 0
80 2.9998310508 3.00011696448 0
81 1.99961793044 2.99991047454 0
82 0.99893758772 3.00012770259 0
83 -0.00142487299093 3.00031379179 0
84 -1.00080405875 3.00022158013 0
85 -2.00010417364 2.9998006105 0
86 -3.00007713367 2.99964591748 0
87 3.00030060993 -1.99923262533 0
88 3.00076625804 -0.999512209558 0
89 3.00062506033 0.000332370132171 0
90 2.99990222364 1.00022112094 0
91 2.99988030663 2.00001186702 0
92 -3.00002810167 1.99915492967 0
93 -3.00048882165 0.999414151808 0
94 -3.00026628885 -0.000458794438716 0
95 -2.99967527287 -1.00033719449 0
96 -2.99964659429 -2.00013795347 0
97 1.99957960894 1.99994320692 0
98 0.999096100565 2.00010863604 0
99 -0.000825745281158 2.00007227662 0
100 -1.0005440083 1.99987241053 0
101 -2.00036986188 1.99935271211 0
102 -1.99928330302 -2.00011998997 0
103 -0.998804038433 -2.00026332116 0
104 0.0011126718878 -2.00023835675 0
105 1.00081734026 -2.00002925918 0
106 2.00063566551 -1.99948223508 0
107 -2.00043392468 0.999573208968 0
108 -2.00000263765 -0.000296906707415 0
109 -1.99957682024 -1.00022376128 0
110 2.00074961191 -0.999733533326 0
111 2.00034435558 0.000128681597118 0
112 1.99986629714 1.00008210862 0
113 -1.00038181027 0.999814113821 0
114 -0.999876784609 -0.000182761960104 0
115 -0.999353901652 -1.000206728 0
116 1.00070179374 -0.999986081415 0
117 1.00022710044 1.76932347369e-05 0
118 0.999667214503 1.00003830359 0
119 0.00067099396712 -1.0001346169 0
120 -0.000329278307687 0.999954348114 0
121 0.000173007524032 -8.63342286404e-05 0
1 1 quad 3 41 40 1
2 1 quad 4 42 41 3
3 1 quad 5 43 42 4
4 1 quad 6 44 43 5
5 1 quad 7 45 44 6
6 1 quad 8 46 45 7
7 1 quad 9 47 46 8
8 1 quad 10 48 47 9
9 1 quad 11 49 48 10
10 1 quad 2 13 49 11
11 1 quad 23 50 21 12
12 1 quad 24 51 50 23
13 1 quad 25 52 51 24
14 1 quad 26 53 52 25
15 1 quad 27 54 53 26
16 1 quad 28 55 54 27
17 1 quad 29 56 55 28
18 1 quad 30 57 56 29
19 1 quad 31 58 57 30
20 1 quad 22 32 58 31
21 1 quad 14 59 49 13
22 1 quad 15 60 59 14
23 1 quad 16 61 60 15
24 1 quad 17 62 61 16
25 1 quad 18 63 62 17
26 1 quad 19 64 63 18
27 1 quad 20 65 64 19
28 1 quad 21 50 65 20
29 1 quad 33 66 58 32
30 1 quad 34 67 66 33
31 1 quad 35 68 67 34
32 1 quad 36 69 68 35
33 1 quad 37 70 69 36
34 1 quad 38 71 70 37
35 1 quad 39 72 71 38
36 1 quad 40 41 72 39
37 1 quad 59 73 48 49
38 1 quad 60 74 73 59
39 1 quad 61 75 74 60
40 1 quad 62 76 75 61
41 1 quad 63 77 76 62
42 1 quad 64 78 77 63
43 1 quad 65 79 78 64
44 1 quad 50 51 79 65
45 1 quad 66 80 57 58
46 1 quad 67 81 80 66
47 1 quad 68 82 81 67
48 1 quad 69 83 82 68
49 1 quad 70 84 83 69
50 1 quad 71 85 84 70
51 1 quad 72 86 85 71
52 1 quad 41 42 86 72
53 1 quad 52 87 79 51
54 1 quad 53 88 87 52
55 1 quad 54 89 88 53
56 1 quad 55 90 89 54
57 1 quad 56 91 90 55
58 1 quad 57 80 91 56
59 1 quad 43 92 86 42
60 1 quad 44 93 92 43
61 1 quad 45 94 93 44
62 1 quad 46 95 94 45
63 1 quad 47 96 95 46
64 1 quad 48 73 96 47
65 1 quad 81 97 91 80
66 1 quad 82 98 97 81
67 1 quad 83 99 98 82
68 1 quad 84 100 99 83
69 1 quad 85 101 100 84
70 1 quad 86 92 101 85
71 1 quad 74 102 96 73
72 1 quad 75 103 102 74
73 1 quad 76 104 103 75
74 1 quad 77 105 104 76
75 1 quad 78 106 105 77
76 1 quad 79 87 106 78
77 1 quad 93 107 101 92
78 1 quad 94 108 107 93
79 1 quad 95 109 108 94
80 1 quad 96 102 109 95
81 1 quad 88 110 106 87
82 1 quad 89 111 110 88
83 1 quad 90 112 111 89
84 1 quad 91 97 112 90
85 1 quad 107 113 100 101
86 1 quad 108 114 113 107
87 1 quad 109 115 114 108
88 1 quad 102 103 115 109
89 1 quad 110 116 105 106
90 1 quad 111 117 116 110
91 1 quad 112 118 117 111
92 1 quad 97 98 118 112
93 1 quad 104 119 115 103
94 1 quad 105 116 119 104
95 1 quad 99 120 118 98
96 1 quad 100 113 120 99
97 1 quad 117 121 119 116
98 1 quad 118 120 121 117
99 1 quad 115 119 121 114
100 1 quad 120 113 114 121
1 1 line 22 32
2 1 line 32 33
3 1 line 33 34
4 1 line 34 35
5 1 line 35 36
6 1 line 36 37
7 1 line 37 38
8 1 line 38 39
9 1 line 39 40
10 1 line 40 1
11 2 line 12 23
12 2 line 23 24
13 2 line 24 25
14 2 line 25 26
15 2 line 26 27
16 2 line 27 28
17 2 line 28 29
18 2 line 29 30
19 2 line 30 31
20 2 line 31 22
21 3 line 2 13
22 3 line 13 14
23 3 line 14 15
24 3 line 15 16
25 3 line 16 17
26 3 line 17 18
27 3 line 18 19
28 3 line 19 20
29 3 line 20 21
30 3 line 21 12
31 4 line 1 3
32 4 line 3 4
33 4 line 4 5
34 4 line 5 6
35 4 line 6 7
36 4 line 7 8
37 4 line 8 9
38 4 line 9 10
39 4 line 10 11
40 4 line 11 2
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
