The differences between the horizontal membership function used in multidimensional fuzzy arithmetic and the inverse membership function used in gradual arithmetic

In the last few years, the number of applications of the multidimensional fuzzy arithmetic (MFA) and the multidimensional interval arithmetic is expanding. Authors of new papers about applications of MFA are often faced with comments from other researchers, especially the gradual arithmetic (GA) proponents, that the horizontal membership function (HMF) used in MFA is the same as the inverse membership function (InvMF) used in GA, and that MFA itself adds nothing new to the fuzzy arithmetic. This view leads to unfair evaluations of scientific papers about MFA applications submitted to scientific journals and to unnecessary disagreements between MFA and GA proponents. The purpose of this paper is to carefully analyze the two types of functions (HMF and InvMF) and to demonstrate their important differences. The basic and decisive difference is the dimensionality of both functions, which is illustrated by examples. It should also be added that HMF has proven its usefulness in solving difficult problems such as: systems of fuzzy equations or fuzzy differential equations, which is confirmed by numerous publications. The paper enable the reader to have a deeper understanding of the multidimensional fuzzy arithmetic.


Introduction
The concept of a horizontal membership function (HMF) is used in the multidimensional fuzzy arithmetic (MFA) Landowski 2013, 2015;Pluciński 2015, 2017), which in first publications was called RDM fuzzy arithmetic (from relative-distance-measure). MFA due to its multidimensional approach to fuzzy calculations is more difficult to understand than 2D types of the fuzzy arithmetic (FA). Therefore, MFA requires a thorough explanation of its essence and details, supported by examples. The scientific team studying MFA has carefully explained how to perform each arithmetic operation in successive papers. Piegat and Landowski (2013) present the implementation of multiplication and division operations, Piegat and Pluciński (2015) present addition, and Piegat and Pluciński (2017) present the division operation, which is slightly more difficult than the other arithmetic operations. Piegat and Landowski (2015) presents the first applications of HMF.
An important advantage of using MFA is its ability to find algebraic solutions of equations, including systems of equations with multiple unknowns. Solving fuzzy equations is very difficult and there are many publications on this subject. Kołodziejczyk et al. (2020) show the advantage of HMF and MFA over one of the methods of solving nonlinear fuzzy equations known in the literature.
Solving differential equations is also a very difficult problem-there are a lot of publications in this area as well. The application of HMF and MFA to this task allowed to obtain precise solutions of differential Andrzej Piegat and Marcin Pluciński are contributed equally to this work. equations, which are very important in the modeling, study and control of dynamic systems. Mazandarani and Zhao (2019), Najariyan and Zhao (2020) and Mazandarani and Xiu (2021) present the results of this type of research.
This paper compares HMF used in MFA and the inverse membership function (InvMF) used to create gradual numbers (GN) used in the gradual interval arithmetic (GIA). GIA can be considered somewhat competitive to MFA. The concept of gradual elements of fuzzy sets was introduced by Dubois and Prade (2008), and the concept of gradual numbers and the fundamentals of GIA were presented by Fortin et al. (2008). Currently, GIA is being developed by, among others (Boukezzoula et al. 2021). An extended notion of GIA-thick gradual intervals, which is equivalent to type-2 fuzzy intervals used in the fuzzy arithmetic, has also been developed (Boukezzoula et al. 2019). In turn, Piegat and Landowski (2021a) present a new MFA based on type-2 HMFs which are extensions of type-1 HMFs. Application examples are also presented in this paper.
The fact, that MFA is more difficult to understand than 2D types of the fuzzy arithmetic has sometimes resulted in erroneous comments about it expressed in some FA papers and in opinions of some reviewers. One of the objections against MFA is the assumption that HMF used in MFA and and inverse membership function (InMF) used in the gradual interval arithmetic (GIA) (Boukezzoula et al. 2021;Dubois and Prade 2008;Fortin et al. 2008) are one and the same. More precisely, that HMF in the theory of gradual numbers is called the assignment function (AF) of a gradual number. The authors of this objection also propose that HMF should be treated as the inverse membership function.
Such opinions have not only theoretical, but also practical implications. Authors of publications on various applications of MFA and HMF inform us about erroneous comments made in reviews of their papers that lead to an unfair evaluation. Erroneous opinions about MFA and HMF are also shared either verbally or by correspondence. Therefore, a strong motivation for the preparation of this paper is to prevent the spread of misinformation and the downgrading of HMF and MFA.
The main contribution of this paper is the detailed analysis of the properties of InvMF used in GIA and HMF used in HMA and the identification of the differences between them. It can be noted at the very beginning that the basic and main difference between InvMF and HMF is the different dimensionality of these two functions resulting from their nature. The higher dimensionality of HMF provides advantages in solving difficult FA problems such as solving systems of equations (Kołodziejczyk et al. 2020), solving fuzzy differential equations (Mazandarani and Zhao 2019;Najariyan and Zhao 2020;Mazandarani and Xiu 2021), implementing operations in type-2 arithmetic (Piegat and Landowski 2021a, b).
In the paper, we want to make it clear that HMF and InvMF are not the same. In Sect. 2, we will explain what HMF is and how it is used in MFA. Sect. 3 will explain the concept of an assignment function, InvMF, a gradual number (GN), a gradual interval (GI) and will show how these terms are used in GIA. Section 4 will explain the differences between HMF and the assignment function. Section 5 will give examples of calculations using HMFs and InMFs and Sect. 6 presents conclusions.
2 Horizontal membership function Figure 1 shows a trapezoidal MF defined by a quadruple ða 1 ; a 2 ; a 3 ; a 4 Þ. This MF presents the knowledge about the value of the variable a that is not known exactly.
The traditional, most commonly used, vertical MF (VMF) lðaÞ is defined by the formula (1). Fig. 1 The trapezoidal membership function is defined by quadruple ða 1 ; a 2 ; a 3 ; a 4 Þ ¼ ð3; 6; 8; 10Þ, where l À ðaÞ means the left border and l þ ðaÞ the right border of MF lðaÞ ¼ for a 2 ½a 1 ; a 2 lðaÞ ¼ l À ðaÞ ¼ a À a 1 a 2 À a 1 for a 2 ½a 2 ; a 3 lðaÞ ¼ 1 for a 2 ½a 3 ; a 4 lðaÞ ¼ l þ ðaÞ ¼ a 4 À a a 4 À a 3 This MF corresponds to the first definition of the fuzzy set (FS) by Zadeh (1965). It defines borders of FS and for each value of the variable a it assigns the degree of membership lðaÞ as the output. On the basis of vertical models of the left l À ðaÞ and right l þ ðaÞ borders, it is possible to calculate inverse functions a À ðlÞ and a þ ðlÞ defined by the formula (2).
a À ðlÞ ¼ a 1 þ lða 2 À a 1 Þ; l 2 ½0; 1 a þ ðlÞ ¼ a 4 À lða 4 À a 3 Þ; l 2 ½0; 1 In inverse functions, the informational output is not the degree of membership l but the value of the variable a, Fig. 2. RDM-variable c a specifies a set of points that lay at an equal, relative distance from the left border a À ðlÞ. Figure 2 shows the exemplary set that is located at the distance c a ¼ 0:5 from the left border. The RDM-variable allows to enter the same (normalized to the range [0,1]) distance of a point located on any a-cut of MF. Using the concept of the RDM-variable it is possible to define HMF (3) for the fuzzy set A.
Formulas (3) and (4) defining HMF allow to specify each point lying inside the fuzzy set A with the coordinates ðl; c a Þ. For example at the level l ¼ 0:5 and c a ¼ 0:5, there is a point with the coordinate a ¼ 6:75, Fig. 2. Since the relation l; c a 2 ½0; 1 is valid, HMF (3) is a mathematical model not only of borders of MF, but also of its entire interior (body) (Fig. 3). The HMF shown in Fig. 3 results from the second definition of the fuzzy set of Zadeh (1975), according to which FS can be composed of its horizontal a-cuts. Modeling of the entire area of FS is possible only due to the fact that in HMF (3) the value of the variable aðl; c a Þ depends not on one but on 2 variables: l and c a . This means that HMF exists in 3D-space, Fig. 4. Figure 3 shows only a projection of this function from 3D to a 2D-space.
HMFs, unlike VMFs, are particularly suitable for performing arithmetic operations, which are simply performed on values of variables a and b that occur in such operations. In the case of HMFs, the objection is sometimes raised that the membership function should be a function having the degree of membership lðaÞ as the output and that the output should not be the value of the variable a itself. In opinion of authors, this objection is rather incorrect,   5,9]. Using HMFs of fuzzy numbers, we can perform arithmetic operations according to the formula (5).
The algebraic result xðl; c a ; c b Þ exists in the 4D-space and it is an universal result, i.e. one that satisfies all possible forms of the mathematical formula in which it occurs. For example, for addition A þ B ¼ X, the result X satisfies the following forms: When performing complex mathematical calculations that are realized in parts, only multidimensional intermediate results should be inserted into the final formula and not the spans of these intermediate results, as it is done today. The span SP x ðlÞ is just a simplified twodimensional information about the full multidimensional result. The span SP x ðlÞ is calculated using the formula (6).
Further information about performing arithmetic operations on FN with HMFs can be found in authors' publications Landowski 2013, 2015;Pluciński 2015, 2017).

Assignment function used in gradual numbers and the gradual arithmetic
As noted in Sect. 1, it is sometimes alleged that HMF used in MFA is the same as InvMF, also referred to as the assignment function (AF), used in the gradual interval arithmetic (GIA). Therefore, it is necessary to explain exactly what AF is. In publications (Boukezzoula et al. 2021;Dubois and Prade 2008;Fortin et al. 2008) their authors presented the concept of GIA, in which gradual intervals (GIs) are used, which contain gradual numbers (GNs) as their components. Detailed definitions of above terms and various explanations related to them are presented by (Boukezzoula et al. 2021;Dubois and Prade 2008;Fortin et al. 2008). In this paper, the GIA concept will be presented only briefly and we will use the original notations from cited papers. The authors of the GIA concept claim that imprecise value (which in FA is called a fuzzy number) should rather be called a fuzzy interval (FI) (Fortin et al. 2008). According to them, FIs are a subclass of GIs, which are a broader concept, because the left border of GI does not have to be a non-decreasing function and the right border does not have to be a non-increasing function, as is the case of FIs. Hence, every FI is GI but not every GI is FI. Figure 5 shows a trapezoidal and a triangular FI which is also GI. The fuzzy interval (FI = GI) has left and right borders being gradual numbers (GNs) by definition. ''A gradual number allows to look at fuzzy intervals as crisp intervals on the set of gradual real numbers'' (Fortin et al. 2008). Such an interpretation of GIs enables the application of all existing methods of interval arithmetic (IA) in GIA. Boukezzoula et al. (2021) list the following methods of IA: 1. Standard interval arithmetic (SIA). 2. Instantiated interval arithmetic (IIA). 3. Extended interval arithmetic (EIA). 4. Non-standard interval arithmetic (NIA). 5. Generalized Hukuhara interval arithmetic (GHIA). 6. Optimistic interval arithmetic (OIA). 7. Constrained interval arithmetic (CIA). 8. Single level constrained interval arithmetic (SCIA). 9. Requisite constrained interval arithmetic (RCIA).
In the case of FIs shown in Fig. 5, their left and right borders can be described using traditional vertical MFs given by formulas (7) (Fortin et al. 2008): The input information in the formula (7) are values of the variable m and the output information are degrees of membership l À ðmÞ and l þ ðmÞ. From vertical MFs (7) we  (8), where a denotes the degree of membership l (a-cuts).
MFs (8) can be called inverse MFs, because the output information here is the value of the variable m corresponding to the membership degree a. In the case of FIs, each border, left or right, corresponds to only one InvMF (in the case of GIs there may be more) and then the whole function can be defined as the ''assignment function'' (Fortin et al. 2008). ''The assignment function Ar, for each value a 2 ð0; 1 assigns a real number r ¼ ArðaÞ.'' AFs allow to define borders of FIs, which (Fortin et al. 2008) call ''gradual real numbers''. The definition of the gradual number is as follows. ''A gradual real number r is defined by an assignment function Ar from (0, 1] (the unit interval minus 0) to the reals.'' Inverse functions ðl À M Þ À1 and ðl þ M Þ À1 are special cases of GNs. Figure 6 shows 2 examples of GNs.
When performing arithmetic operations on FIs in accordance with the rules of GIA, calculations are performed on borders of FIs, i.e. on GNs. Further on, definitions of addition and multiplication of GNs (Fortin et al. 2008) will be given.

Addition of GNs
The sum of GNsr ands is defined by summing their assignment functions. Forr þs and a 2 ð0; 1 we have: Ar þs ðaÞ ¼ ArðaÞ þ AsðaÞ : An example of adding GNs from Fig. 6 is given by the formula (10).

Multiplication of GNs
The product and the quotient of GNs can be defined in a similar way (with caution when dividing by a gradual numberr such that ArðaÞ ¼ 0 for some a). If two GNsr ands are given (10), then their product Ars is determined by the formula (11).

Gradual interval arithmetic
Knowing the arithmetic of GNs, it is now possible to present the arithmetic of FIs based on GIA principles (FI is a special case of GI). In the GIA arithmetic operations are performed similarly to the interval arithmetic, using for calculations one of existing types of the interval arithmetic. Boukezzoula et al. (2021) describes this principle as follows. ''If borders of conventional interval ½a ¼ ½a À ; a þ are replaced by assignment functions a À ðaÞ and a þ ðaÞ representing two GNs, we get gradual interval ½aðaÞ ¼ ½a À ðaÞ; a þ ðaÞ which borders are used in interval arithmetic operations''. Figure 7 shows two FIs that will be used in the addition example. Fuzzy interval ½aðaÞ is defined by two border GNs represented by two assignment functions a À ðaÞ and a þ ðaÞ given by (12).
The sum of intervals ½aðaÞ þ ½bðaÞ ¼ ½xðaÞ is given by the formula (14) and shown in Fig. 8.
Note, that adding FIs ½aðaÞ and ½bðaÞ using GIA resulted in a fuzzy interval ½xðaÞ containing not one but two assignment functions. Next, we will see if it is the same in calculations using HMFs.

Comparison of the assignment function used in GIA and horizontal membership function used in MFA
There are differences between AF and HMF which will be presented now. The mathematical model of the fuzzy interval used in GIA consists of two assignment functions a À ðaÞ and a þ ðaÞ: ½aðaÞ ¼ ½a À ðaÞ; a þ ðaÞ; a 2 ð0; 1 : In the case of the trapezoidal fuzzy interval from Fig. 7, it takes the form (16).
The model ½aðaÞ is 2-dimensional, see Fig. 7. In the case of MFA, the model of FI is called HMF and is 3-dimensional, (17).
A : a 2 A aða; c a Þ ¼ a À ðaÞ þ c a ½a þ ðaÞ À a À ðaÞ ¼ ð1 þ 2aÞ þ c a ð5 À 4aÞ a; c a 2 ½0; 1 ð17Þ The difference between the two mathematical models is illustrated in Fig. 9. Gradual interval ½aðaÞ represents only borders a À ðaÞ and a þ ðaÞ of the fuzzy interval. For this reason, only borders are used in arithmetic calculations. However, the interior of FI between these borders is not included in calculations. This affects calculation results in many cases. Using HMF, which represents both borders and the interior of FI, we obtain sometimes identical and sometimes different results than in the case of using gradual interval models ½aðaÞ. Examples illustrating this fact can be found in publications of the authors on MFA. Examples will also be given in the next section to illustrate the difference in results of both methods. In gradual arithmetic, borders (AFs) of FI a À ðaÞ and a þ ðaÞ take part in arithmetic operations separately. In the case of HMFs, these borders are combined in a uniform function aða; c a Þ ¼ a À ðaÞ þ c a ½a þ ðaÞ À a À ðaÞ forming a geometrically single whole and only this unity takes part in calculations. This has a significant impact on obtained results.

Examples illustrating differences in results
obtained with the use of gradual interval and horizontal membership function models First example will consist in solving the equation A þ X ¼ B, where fuzzy intervals A and B represent the knowledge about the uncertain values of the variables a and b. The intervals ½aðaÞ and ½bðaÞ are shown in Fig. 7. These intervals have the following form: ½aðaÞ ¼ ½a À ðaÞ; a þ ðaÞ ¼ ½1 þ 2a; 6 À 2a and ½bðaÞ ¼ ½b À ðaÞ; b þ ðaÞ ¼ ½5 þ 2a; 9 À 2a. Interval insides include AFs of left and right borders of intervals. First, the equation will be solved using GIA, and SIA of Warmus-Sunaga-Moore will be used as the interval arithmetic. It should be noted that in the case of GIA it is assumed that the solution of the equation is 2D gradual interval ½xðaÞ ¼ ½x À ðaÞ; x þ ðaÞ, i.e. the same mathematical form as in the case of intervals ½a À ðaÞ; a þ ðaÞ and ½b À ðaÞ; b þ ðaÞ.
The solution of Eq. (19) is x À ðaÞ ¼ 4 and x þ ðaÞ ¼ 3, which is the improper interval ½xðaÞ ¼ ½x À ðaÞ; x þ ðaÞ ¼ ½4; 3. This 'solution' means that there is no real value of x satisfying equation (18), because the value of x would have to satisfy the condition x ! 4 and x 3, what is impossible. However, let's check, using the common sense method of Instance Testing, whether there is actually any solution to the equation under study. At the level a ¼ 0 the values a ¼ 1 and b ¼ 9 are possible, Fig. 8 Result of the addition of two fuzzy intervals from Fig. 7 using the gradual interval arithmetic based on assignment functions resulting in x ¼ b À a ¼ 8. Hence, at least one instance solution is possible: ða; b; x ¼ b À aÞ ¼ ð1; 9; 8Þ. We will see further that there are not only this one, but an infinite number of other possible solutions of the equation under study. For this purpose, MFA and HMFs will be used. The formulation of the task in terms of MFA is given by the formula (20).
A : aða; c a Þ 2 A; B : bða; The solution of the problem is the set ðA; B; X ¼ B À AÞ of triples ða; b; x ¼ b À aÞ defined by the formula (21).
The variable xða; c a ; c b Þ is a result variable that exists in 4D-space and cannot be seen directly. However, it is possible to determine and see the span SP x of this variable defined by the formula (22).
where: xða; c a ; c b Þ ¼ 4 À c a ð5 À 4aÞ þ c b ð4 À 4aÞ, see formula (20). Minimum xða; c a ; c b Þ is obtained for c a ¼ 1, c b ¼ 0 and maximum for c a ¼ 0, c b ¼ 1. Finally, the span of the result variable is given by the formula (23) and its visualization is shown in Fig. 10.
Formula (21) makes it possible to generate all possible states of the system realizing the dependence a þ x ¼ b for each level of membership l ¼ a. For a ¼ 0 possible states are e.g. ða; b; x ¼ b À aÞ: ð6; 5; À1Þ, (3, 6, 3), (1, 5, 4), (1,9,8). In general, there are an infinite number of possible states ða; b; x ¼ b À aÞ satisfying the relationship A þ X ¼ B. Why was GIA based on gradual intervals unable to detect any of these states? This is due to the misconception that the correct span of the result variable x can always be detected by 2D computations without finding a multidimensional solution, as MFA does. Kołodziejczyk et al. (2020) present an example of the use of HMF and MFA to solve system of equations of the order higher than 1. The paper shows that MFA made it possible to find algebraic solutions in tasks where other kinds of fuzzy arithmetic have failed.
In the second example, a real object moving in uniformly decelerated motion will be considered. Let us assume a deceleration a ¼ À0:5 m s 2 and an initial velocity V 0 ¼ 22 m s . Our task will be to determine what distance the vehicle will travel after about 50s of movement.
Time of movement is defined imprecisely as 'about 50s' and we will define it using the triangular fuzzy number shown in Fig. 11. We can describe the time using HMF as: The distance can be calculated from the formula: After substituting time t, we have: dðl; c t Þ ¼ 0:5að40 þ 10l þ c t ð20 À 20lÞÞ 2 þ V 0 ð40 þ 10l þ c t ð20 À 20lÞÞ : The solution dðl; c t Þ is a function of two variables: l and c t . Now, we can determine the span: ; c t 2 ½0; 1 : The plot of the span is presented in Fig. 12. We can perform similar calculations using GIA. The time defined using InvMF has the form: t ¼ ½40 þ 10l; 60 À 10l : We can now determine formulas for the left and right borders of the resulting membership function, which is shown in Fig. 13.
From the Fig. 13, we can see that the resulting fuzzy number determined using GIA is not consistent with the definition of a fuzzy number, unlike the result determined using MFA (Fig. 12). This is because MFA takes into account not only borders but also interiors of the fuzzy numbers in its calculations.

Conclusions
The paper presents a discussion with the view found among some fuzzy arithmetic researchers proclaiming the identity of the inverse membership function used in GIA and the horizontal membership function used in MFA. The similarity of InvMF and HMF also indirectly suggests the similarity (or even identity) of GIA and HMA. The paper shows that there are fundamental differences between InvMF (which is part of the gradual number used in GIA) and HMF. The gradual number is a mathematical object defined in 2D-space, and the InvMFs that are its elements (GN borders) are also 2D objects. In gradual arithmetic, arithmetic calculations are performed using boundaries of GNs. On the other hand, HMF is a mathematical object defined in 3D-space, including not only the borders but also the interior of the membership function. This has a great impact on the computational results that are obtained.
AFs, as models of borders of MF, are used in GIA to calculate the results separately, while HMF is a single overall model of MF which takes part in the calculations as one whole, without dividing it into parts. In general, this often leads to different computational results provided by GIs-based and HMFs-based fuzzy arithmetic types. First of all, the results of calculations in MFA are objects with dimensionality 3 or higher, while GIA always provides results with dimensionality 2. This fact causes MFA to be able to provide, for example, algebraic solutions to equations that verify them after substitution. In the case of GIA, this is not always true, as shown by the simple example presented in the paper. The spans of multidimensional solutions provided by MFA are always consistent with the definition of a fuzzy number, while the solutions provided by GIA are sometimes inconsistent, as shown in the example and what is reported in the literature (Boukezzoula et al. 2021).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons. org/licenses/by/4.0/.