TFG Albert Font [Domain Decomposition]

AlbertTFG files/.~lock.k_longer.csv#

@@ -0,0 +1 @@
+,LEGIONALBERT/fontg,LegionAlbert,27.05.2025 16:20,file:///C:/Users/fontg/AppData/Roaming/LibreOffice/4;

AlbertTFG files/CoarsePlotSolution.m

@@ -0,0 +1,224 @@
+classdef CoarsePlotSolution < handle
+    %UNTITLED Summary of this class goes here
+    %   Detailed explanation goes here
+
+    properties
+        mesh
+        u
+        originalBCApplier
+        outputFileName
+        r
+        centroids
+    end
+
+    methods
+        function obj = CoarsePlotSolution(u, mesh, bcApplier,outputFileName, r, centroids)
+
+            if isempty(r)
+                obj.plot(u, mesh, bcApplier,outputFileName);
+            else
+                obj.init(u, mesh, bcApplier,outputFileName,r, centroids);
+                obj.makeHole(bcApplier,outputFileName);
+
+            end
+        end
+
+        % function makeHole(~, x, mesh, bcApplier,outputFileName,r,holeCoords)
+        %     x = obj.createXHole(x, mesh, r, holeCoords);
+        %     mesh = createHole(mesh, r, holeCoords);
+        %     obj.plot(x, mesh, bcApplier,outputFileName);
+        %o
+        % end
+
+
+
+    end
+
+    methods (Access = private) %%% treure Static
+
+        function makeHole(obj, bcApplier,outputFileName)
+
+            %uH = obj.createUHole();
+            uMeshFun = obj.createUnfittedMeshFun();
+            obj.plot(uMeshFun,outputFileName);
+
+        end
+
+        function init(obj, u, mesh, bcApplier,outputFileName,r,centroids)
+            obj.mesh        = mesh;
+            obj.u           = u;
+            obj.originalBCApplier   = bcApplier;
+            obj.outputFileName      = outputFileName;
+            obj.r                   = r;
+            obj.centroids   = centroids;
+
+        end
+
+        function U = createUHole(obj)
+            U = obj.u;
+
+            for i = obj.mesh.nnodes:-1:1
+                if sqrt( (obj.mesh.coord(i,1) )^2 + (obj.mesh.coord(i,2))^2 ) <= obj.r
+                    U(2*i-1:2*i) = [];
+                end
+
+            end
+
+        end
+
+        function uMeshFun = createUnfittedMeshFun(obj)
+            % mR    = obj.createReferenceMesh();
+
+            gPar.type         = 'Circle';
+
+            % Initialize the sum as the zero function
+            f_sum = @(x) 0;
+            u = obj.u;
+            mesh = obj.mesh;
+            f = @(x) 0;
+            s.x0 = obj.centroids(:,1);
+            s.y0 = obj.centroids(:,2);
+            s.r  = obj.r;
+            s.type = 'Circles';
+            g = GeometricalFunction(s);
+            ls = g.computeLevelSetFunction(obj.mesh);
+
+            lsCircle          = ls.fValues;
+            lsCircleInclusion = -lsCircle;
+
+            sUm.backgroundMesh = obj.mesh;
+            sUm.boundaryMesh   = obj.mesh.createBoundaryMesh;
+            uMesh              = UnfittedMesh(sUm);
+            uMesh.compute(lsCircleInclusion);
+            uMeshFun = uMesh.obtainFunctionAtUnfittedMesh(u);
+%             s.fValues = reshape(uMeshFun.innerMeshFunction.fValues,2,[])';
+
+        end
+
+        function newMesh = createReferenceMesh(obj)
+            xmin = -1;
+            xmax = 1;
+            ymin = -1;
+            ymax = 1;
+
+            n1 = numel(find(obj.mesh.coord(:,1)==xmin));
+            n2 = n1;%numel(find(obj.mesh.coord(:,2)==ymin));
+
+
+            %UnitMesh better
+            x1      = linspace(xmin,xmax,n1);
+            x2      = linspace(ymin,ymax,n2);
+            [xv,yv] = meshgrid(x1,x2);
+            [F,V]   = mesh2tri(xv,yv,zeros(size(xv)),'x');
+            s.coord  = V(:,1:2);
+            s.connec = F;
+            newMesh = Mesh.create(s);
+        end
+
+        function levelSet = createLevelSetFunction(obj,mR)
+            sLS.type        = 'CircleInclusion';
+            sLS.xCoorCenter = 0;
+            sLS.yCoorCenter = 0;
+            sLS.radius      = obj.r;
+            g               = GeometricalFunction(sLS);
+            lsFun           = g.computeLevelSetFunction(mR);
+            levelSet        = lsFun.fValues;
+        end
+
+        function uMesh = computeUnfittedMesh(obj, mR,levelSet)
+            sUm.backgroundMesh = mR;
+            sUm.boundaryMesh   = mR.createBoundaryMesh();
+            uMesh              = UnfittedMesh(sUm);
+            uMesh.compute(levelSet);
+        end
+
+        function plot(~, uMeshFun,outputFileName)
+            row = 0;
+            col = 0;
+            iter = 0;
+            flag = 0;
+
+%             if ~isempty(bcApplier)
+%                 x = bcApplier.reducedToFullVectorDirichlet(x);
+%             end
+%             if nargin <7
+%                 flag =0;
+%             end
+%             %             xFull = bc.reducedToFullVector(x);
+%             if size(x,2)==1
+%                 s.fValues = reshape(x,2,[])';
+%             else
+%                 s.fValues = x;
+%             end
+%             %
+% 
+%             s.mesh = mesh;
+%             s.fValues(:,end+1) = 0;
+%             s.ndimf = 2;
+%             s.order = 'P1';
+%             xF = LagrangianFunction(s);
+%             %             xF.plot();
+%             if flag == 0
+%                 xF.print(['domain',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+%             elseif flag == 1
+%                 xF.print(['DomainResidual',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+%             elseif flag == 2
+%                 xF.print(['Residual',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+%             elseif flag == 3
+%                 xF.print(['domainFine',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+%             elseif flag == 4
+%                 xF.print(['domainNeuman',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+%             end
+
+
+
+            %uMeshFun.innerMeshFunction.print(['inner',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+            %uMeshFun.innerCutMeshFunction.print(['innerCut',num2str(row),num2str(col),'_',num2str(iter)],'Paraview')
+
+            fvalues = [uMeshFun.innerMeshFunction.fValues;
+                         uMeshFun.innerCutMeshFunction.fValues];
+
+            s.coord = [uMeshFun.innerMeshFunction.mesh.coord;
+                         uMeshFun.innerCutMeshFunction.mesh.coord];
+
+            s.connec = [uMeshFun.innerMeshFunction.mesh.connec;
+                         uMeshFun.innerCutMeshFunction.mesh.connec  + max(uMeshFun.innerMeshFunction.mesh.connec(:))];
+
+            mh = Mesh.create(s);
+
+            ss.mesh = mh;
+            ss.fValues = fvalues;
+            ss.order = 'P1';
+            ss.ndimf = size(fvalues,2)
+
+            u = LagrangianFunction(ss);
+
+            u.print([outputFileName],'Paraview')
+
+            % fileName = [num2str(row),num2str(col),'_',num2str(iter)];
+            % 
+            % s = dir(pwd);
+            % s = struct2table(s);
+            % idx = startsWith(s.name, fileName);
+            % s = s(idx,:);
+            % oldFileName = s.name;
+            % newFileName = replace(oldFileName, fileName, outputFileName);
+            % fclose('all');
+            % movefile(oldFileName{1}, newFileName{1});
+            % 
+            % fileName = ['innerCut',num2str(row),num2str(col),'_',num2str(iter)];
+            % 
+            % s = dir(pwd);
+            % s = struct2table(s);
+            % idx = startsWith(s.name, fileName);
+            % s = s(idx,:);
+            % oldFileName = s.name;
+            % newFileName = replace(oldFileName, fileName, [outputFileName, '_cut']);
+            % fclose('all');
+            % movefile(oldFileName{1}, newFileName{1});
+
+        end
+
+
+    end
+end