You can subscribe to this list here.
2010 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}

_{Oct}
(1) 
_{Nov}
(1) 
_{Dec}


2011 
_{Jan}

_{Feb}

_{Mar}
(1) 
_{Apr}

_{May}
(1) 
_{Jun}
(9) 
_{Jul}
(2) 
_{Aug}
(6) 
_{Sep}
(11) 
_{Oct}
(15) 
_{Nov}
(4) 
_{Dec}
(9) 
2012 
_{Jan}
(7) 
_{Feb}
(14) 
_{Mar}
(11) 
_{Apr}

_{May}

_{Jun}
(3) 
_{Jul}

_{Aug}

_{Sep}

_{Oct}

_{Nov}
(3) 
_{Dec}

2013 
_{Jan}
(8) 
_{Feb}

_{Mar}
(1) 
_{Apr}
(3) 
_{May}

_{Jun}

_{Jul}
(3) 
_{Aug}
(2) 
_{Sep}

_{Oct}

_{Nov}

_{Dec}
(4) 
2014 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}
(2) 
_{Jun}
(2) 
_{Jul}

_{Aug}

_{Sep}
(5) 
_{Oct}
(3) 
_{Nov}
(1) 
_{Dec}
(4) 
2015 
_{Jan}
(5) 
_{Feb}
(10) 
_{Mar}
(12) 
_{Apr}
(14) 
_{May}
(26) 
_{Jun}
(7) 
_{Jul}

_{Aug}

_{Sep}
(4) 
_{Oct}

_{Nov}

_{Dec}

2016 
_{Jan}
(1) 
_{Feb}
(9) 
_{Mar}
(3) 
_{Apr}
(3) 
_{May}
(12) 
_{Jun}
(2) 
_{Jul}
(13) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 

1

2

3

4

5

6

7

8

9
(2) 
10

11

12

13

14

15

16
(2) 
17

18

19

20

21

22

23

24

25

26

27

28
(1) 
29

30

31
(2) 




From: Rafael Vázquez <vazquez@im...>  20120131 18:23:55

Dear Alexis, if all the data is correct (boundary conditions, source function...), the reason is probably that you are not matching the interfaces in the correct way. Have you checked the computed solution in paraview? Is it continuous along the interfaces? Rafa 
From: alexis papagiannopoulos <alexis.papagiannopoulos@gm...>  20120131 12:18:28

Dear Geopde users, I am currently working on a multipatch imposed onto an ellipses domain ,which is constructed by two symmetrical patches(semiellipses along x).the geometry file is the following: # HEADER 2 2 3 1 2 1 5 2 0.0000 0.0000 0.0000 0.5000 0.5000 1.0000 1.0000 1.0000 0.0000 0.0000 1.0000 1.0000 0.7500 0.3750 0.0000 0.3750 0.7500 1.0000 0.7071 0.0000 0.7071 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.3536 0.5000 0.3536 0.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 0.7071 1.0000 0.7071 1.0000 2 1 5 2 0.0000 0.0000 0.0000 0.5000 0.5000 1.0000 1.0000 1.0000 0.0000 0.0000 1.0000 1.0000 1.0000 0.7071 0.0000 0.7071 1.0000 0.7500 0.3750 0.0000 0.3750 0.7500 0.0000 0.3536 0.5000 0.3536 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.7071 1.0000 0.7071 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 INTERFACE 1 1 1 2 1 1 INTERFACE 2 1 4 2 3 1 INTERFACE 3 1 2 2 2 1 SUBDOMAIN 1 1 2 BOUNDARY 1 1 1 3 BOUNDARY 2 1 2 4 and the function file used: function [error_h1,error_h2,dofs] = ... ellipses_multi(i) clear problem_data problem_data.geo_name = 'ellipses.txt'; % Type of boundary conditions for each side of the domain problem_data.nmnn_sides = [ ]; problem_data.drchlt_sides = [1 2]; % Physical parameters problem_data.c_diff = @(x, y) ones(size(x)); % Source and boundary terms problem_data.f = @(x, y) exp(x).*((x.^2 + y.^2  1).*sin(x.*y)  2*y.*cos(x.*y)); problem_data.g = @(x, y, ind) test_ellipses_bc_g_nmnn (x, y, ind); problem_data.h = @(x, y, ind) exp(x) .* sin (x.*y); % Exact solution (optional) problem_data.uex = @(x, y) exp(x) .* sin (x.*y); problem_data.graduex = @(x, y) cat (1, ... reshape ( exp(x).*(sin(x.*y) + y.*cos(x.*y)), [1, size(x)]), ... reshape (exp(x).*x.*cos(x.*y), [1, size(x)])); data_names = fieldnames (problem_data); for iopt = 1:numel (data_names) eval ([data_names{iopt} '= problem_data.(data_names{iopt});']); end % Construct geometry structure, and information for interfaces and boundaries [geometry, boundaries, interfaces] = mp_geo_load (geo_name); npatch = numel (geometry); geometry=mp_hrefinement(geometry,[i i]); msh = cell (1, npatch); sp = cell (1, npatch); for iptc = 1:npatch % Define the refined mesh, with tensor product structure knots{iptc}=geometry(iptc).nurbs.knots; % Compute the quadrature rule rule = msh_gauss_nodes (geometry(iptc).nurbs.order); [qn, qw] = msh_set_quad_nodes (knots{iptc}, rule); msh{iptc} = msh_2d (knots{iptc}, qn, qw, geometry(iptc)); % Evaluate the discrete space basis functions in the quadrature points sp{iptc} = sp_nurbs_2d(geometry(iptc).nurbs, msh{iptc}); dofs=sp{iptc}.ndof; end % Create a correspondence between patches on the interfaces [gnum, ndof] = mp_interface_2d (interfaces, sp); % Compute and assemble the matrices rhs = zeros (ndof, 1); ncounter = 0; for iptc = 1:npatch [rs, cs, vs] = op_gradu_gradv_tp (sp{iptc}, sp{iptc}, msh{iptc}, c_diff); rows(ncounter+(1:numel (rs))) = gnum{iptc}(rs); cols(ncounter+(1:numel (rs))) = gnum{iptc}(cs); vals(ncounter+(1:numel (rs))) = vs; ncounter = ncounter + numel (rs); rhs_loc = op_f_v_tp (sp{iptc}, msh{iptc}, f); rhs(gnum{iptc}) = rhs(gnum{iptc}) + rhs_loc; end stiff_mat = sparse (rows, cols, vals, ndof, ndof); for iref = nmnn_sides for bnd_side = 1:boundaries(iref).nsides iptc = boundaries(iref).patches(bnd_side); iside = boundaries(iref).faces(bnd_side); msh_side = msh_eval_boundary_side (msh{iptc}, iside); sp_side = sp_eval_boundary_side (sp{iptc}, msh_side); x = squeeze (msh_side.geo_map(1,:,:)); y = squeeze (msh_side.geo_map(2,:,:)); gval = reshape (g (x, y, iref), msh_side.nqn, msh_side.nel); rhs_nmnn = op_f_v (sp_side, msh_side, gval); global_dofs = gnum{iptc}(sp_side.dofs); rhs(global_dofs) = rhs(global_dofs) + rhs_nmnn; end end % Apply Dirichlet boundary conditions u = zeros (ndof, 1); [u_drchlt, drchlt_dofs] = mp_sp_drchlt_l2_proj (sp, msh, h, gnum, boundaries, drchlt_sides); u(drchlt_dofs) = u_drchlt; int_dofs = setdiff (1:ndof, drchlt_dofs); rhs(int_dofs) = rhs(int_dofs)  stiff_mat(int_dofs, drchlt_dofs)*u_drchlt; % Solve the linear system u(int_dofs) = stiff_mat(int_dofs, int_dofs) \ rhs(int_dofs); % 4) POSTPROCESSING % 4.1) EXPORT TO PARAVIEW output_file = 'ellipses'; vtk_pts = {linspace(0, 1, 20), linspace(0, 1, 20)}; fprintf ('The result is saved in the file %s.pvd \n \n', output_file); mp_sp_to_vtk (u, sp, geometry, gnum, vtk_pts, output_file, 'u') % 4.2) COMPARISON WITH THE EXACT SOLUTION npatch = numel (geometry); for iptc = 1:npatch [error_h1(iptc), error_l2(iptc)] = sp_h1_error (sp{iptc}, msh{iptc}, ... u(gnum{iptc}), problem_data.uex, problem_data.graduex); end error_l2 = sqrt (sum (error_l2 .* error_l2)) error_h1 = sqrt (sum (error_h1 .* error_h1)) end However when i test the function for several levels of hrefinement(in both patches simultaneously) i get a relatively stable error. Any suggestions why that could be happpening? Thank you, Aleix Alexis Papagiannopoulos 
From: Dominik Mokriš <dominik.tramvaj@gm...>  20120128 22:11:05

It seems to be working now, however >pkg uninstall geopdes_base nurbs you suggested, wasn't enough to uninstall. I don't remember, what exactly happened but the octave kept telling, nurbs and geopdes are still installed. I have made a note somewhere, I hope to find it and tell you, how I resolved this. So the conclusion: Users of Ubuntu 10.04, use >pkg install local nurbs1.3.5.tar.gz >pkg install local geopdes_base2.0.1.tar.gz instead of usual >pkg install ( ... ) In default package repositories I've found only older version of the nurbs package, I think it's 1.0.2, so I don't recommend getting the nurbs package this way. Thanks for your help and sorry for delay with an answer. Dominik On Sat, Dec 24, 2011 at 3:18 PM, c. <carlo.defalco@...> wrote: > > > On 24 Dec 2011, at 14:36, Dominik Mokriš wrote: > > > Hello all, > > > > recently I've started using Ubuntu 10.04 LTS and got into trouble with the > > GeoPDEs installation. Here's how: > > > > What is the trouble? I got an error similar to these warnings on > > Ubuntu10.10. > > It seems you are trying to install the packages systemwide rather than locally, > are you running Octave as root? This is not recommended as it might conflict with packages > installed by your system package manager [*]. > > try uninstalling all packages then run "pkg install" again with the option "local", e.g. > > pkg uninstall geopdes_base nurbs > pkg install local nurbs1.3.5.tar.gz > pkg install local geopdes_base2.0.1.tar.gz > > > Thanks for any suggestions and have a nice Christmas. > > Dominik Mokriš, Charles University in Prague > > HTH, > Carlo > > [*] If you want to install packages systemwide on Ubuntu I'd recommend using synaptic or aptget instead, e.g.: > > sudo aptget install octavenurbs > > unfortunately geopdes is not distributed this way yet ... 
From: Rafael Vazquez <vazquez@im...>  20120116 16:40:35

Dear Alexis, the .mat files are not supported in multipatch because the information from the nurbs structure is not enough, since you also need the boundaries and interfaces. But since this information does not change with the refinement, you can load the refined geometry with the function geo_load (without mp) for each single patch. This is exactly what is done in the examples in geopdes_elasticity/multipatch. Regards, Rafa Il 16/01/2012 14:11, alexis papagiannopoulos ha scritto: > Dear geopde users, > > I am working on multipatch and i want graph some convergence tests. To do > this i would like to try hrefinement simultaneously.So i used te > following code: > > function geometry=mp_hrefinement(mp_geo,knts ) > > for j=1:numel(mp_geo) > nurbs=mp_geo(j).nurbs; > [rknots,zeta,nknots]=kntrefine(nurbs.knots,knts,nurbs.order1,[0 0]); > nurbs=nrbkntins(nurbs,nknots); > end > geometry=mp_geo_load(mp_geo.nurbs); > > > end > > > but as mp_geo_load only accepts only txt files it is not possible. Is there > anyway to convert txt files to mat files and vice versa?Any suggestions are > gratefully accepted. > > Alexis Papagiannopoulos. >  > RSA(R) Conference 2012 > Mar 27  Feb 2 > Save $400 by Jan. 27 > Register now! > http://p.sf.net/sfu/rsasfdev2dev2 > _______________________________________________ > Geopdesusers mailing list > Geopdesusers@... > https://lists.sourceforge.net/lists/listinfo/geopdesusers > 
From: alexis papagiannopoulos <alexis.papagiannopoulos@gm...>  20120116 13:11:06

Dear geopde users, I am working on multipatch and i want graph some convergence tests. To do this i would like to try hrefinement simultaneously.So i used te following code: function geometry=mp_hrefinement(mp_geo,knts ) for j=1:numel(mp_geo) nurbs=mp_geo(j).nurbs; [rknots,zeta,nknots]=kntrefine(nurbs.knots,knts,nurbs.order1,[0 0]); nurbs=nrbkntins(nurbs,nknots); end geometry=mp_geo_load(mp_geo.nurbs); end but as mp_geo_load only accepts only txt files it is not possible. Is there anyway to convert txt files to mat files and vice versa?Any suggestions are gratefully accepted. Alexis Papagiannopoulos. 
From: Rafael Vazquez <vazquez@im...>  20120109 11:56:34

Hi Alexis, in the mailing list we didn't get the images, but by your description of the geometry, probably this is similar to what Dominik was doing for the circle, in September. You can read those emails and the article I recommended to him, to see if it is helpful. In general, if your geometry has a side that is degenerated into a point, you don't need to do anything special in the txt file, but you will need to add some lines in the code to force the continuity of the solution (see the article). Regards, Rafa Il 09/01/2012 12:27, alexis papagiannopoulos ha scritto: > Dear Geopde users, > > I am currently working on a multipatch problem, regarding an ellipseslike > domain consisting of 2 patches(patch1.jpg and patch2.jpg). each of the > patches consists of 2 sides. How does the txt configuration work in this > case , and generally if the patch consists of less than 4 sides? > > > >  > Ridiculously easy VDI. With Citrix VDIinaBox, you don't need a complex > infrastructure or vast IT resources to deliver seamless, secure access to > virtual desktops. With this allinone solution, easily deploy virtual > desktops for less than the cost of PCs and save 60% on VDI infrastructure > costs. Try it free! http://p.sf.net/sfu/CitrixVDIinabox > > > _______________________________________________ > Geopdesusers mailing list > Geopdesusers@... > https://lists.sourceforge.net/lists/listinfo/geopdesusers > 
From: alexis papagiannopoulos <alexis.papagiannopoulos@gm...>  20120109 11:28:03

Dear Geopde users, I am currently working on a multipatch problem, regarding an ellipseslike domain consisting of 2 patches(patch1.jpg and patch2.jpg). each of the patches consists of 2 sides. How does the txt configuration work in this case , and generally if the patch consists of less than 4 sides? 