Solved. Thanks all. Evan
On Mon, Nov 22, 2010 at 12:42 AM, Markus Bürg <[email protected]> wrote: > 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 > >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
