You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: Hdf5ElementForceQuad4.m
    path of file: ./repos/NASTRAN_CoFE/nastran_cofe/hdf5_output_processors
    
    the code of the file until where you have to start completion: %Hdf5ElementForceQuad4 HDF5 data class for QUAD4 force data.

% http://web.mscsoftware.com/doc/nastran/2018.2/release/DataType_v20182.html
% https://web.mscsoftware.com/doc/nastran/2018.2/release/DataType_v20182.xml
% <typedef name="QUAD4_FORCE">
% <field name="MX" type="double" description="Membrane force in x"/>
% <field name="MY" type="double" description="Membrane force in y"/>
% <field name="MXY" type="double" description="Membrane force in xy"/>
% <field name="BMX" type="double" description="Bending moment in x"/>
% <field name="BMY" type="double" description="Bending moment in y"/>
% <field name="BMXY" type="double" description="Bending moment in xy"/>
% <field name="TX" type="double" description="Shear force in x"/>
% <field name="TY" type="double" description="Shear force in y"/>
% </typedef>

% A. Ricciardi
% June 2021

classdef Hdf5ElementForceQuad4 < Hdf5ElementForce
    
    properties
        EID % [uint32] Element identification number
        MX % [double] Membrane force in x
        MY % [double] Membrane force in y
        MXY % [double] Membrane force in xy
        BMX % [double] Bending moment in x
        BMY % [double] Bending moment in y
        BMXY % [double] Bending moment in xy
        TX % [doubl
    
function obj = constructFromElementOutputData(elementOutputData,domainIDs)
            % Function to convert element force output data to HDF5
            %
            % INPUTS
            % elementOutputData [nElements,1 ElementOutputData] element force output data
            % domainIDs [nVectors,1 unit32] HDF5 domain ID numbers
            %
            obj = Hdf5ElementForceQuad4();
            nElements = size(elementOutputData,1);
            nVectors = size(elementOutputData(1).values,2);
            eid = [];
            mx = [];
            my = [];
            mxy = [];
            bmx = [];
            bmy = [];
            bmxy = [];
            tx = [];
            ty = [];
            domain_id = [];
            for i = 1:nElements
                eid = [eid,repmat(elementOutputData(i).elementID,[1,nVectors])];
                
                mx = [mx,elementOutputData(i).values(1,:)];
                my = [my,elementOutputData(i).values(2,:)];
                mxy = [mxy,elementOutputData(i).values(3,:)];
                
                bmx = [bmx,elementOutputData(i).values(4,:)];
                bmy = [bmy,elementOutputData(i).values(5,:)];
                bmxy = [bmxy,elementOutputData(i).values(6,:)];
                tx = [tx,elementOutputData(i).values(7,:)];
                ty = [ty,elementOutputData(i).values(8,:)];
                
                domain_id = [domain_id,domainIDs];
            end