search

KikeM / romtime

Unsteady Finite Element Reduced Order Model for Time Moving Domains
1 stars 1 forks source link

BUG : Mass Matrix Interpolation Overridden (ROM) #22

Closed KikeM closed 3 years ago

KikeM commented 3 years ago

I closed KikeM/msc-thesis#19 in PR-#34 because I thought I had correctly implemented the MDEIM for the mass matrix too. There was a bug in the code, the ROM interpolation was being computed and then overwritten by the FOM assembly and projection.

def assemble_mass(self, mu, t):

       if self.mdeim_Mh:            
            MN = self.mdeim_Mh.interpolate(mu=mu, t=t, which=self.mdeim_Mh.ROM)        

        # Assemble FOM operator        
        Mh = self.fom.assemble_mass(mu, t)        
        MN = self.to_rom_bilinear(Mh)

        return MN