A computer-implemented method for simulating an assembly of fabric pieces comprising the steps of: - providing a plurality of parallel solvers, a first solver simulating at least one first fabric piece and a second solver simulating at least one second fabric piece, the at least one first and second fabric pieces being modeled by surfaces meshes comprising nodes; and - maintaining the at least one first and second fabric pieces tight by computing a common frontier by the first and second solvers, wherein the common frontier comprises paired nodes of the at least one first and second fabric pieces, and wherein the computation of the common frontier is performed by combining data of the first fabric piece and data of the second fabric piece, the combining comprising reaching a common position target of each of said paired nodes; wherein the combination of data of the first fabric piece and the second fabric piece is performed by mutually exchanging data between the first solver and the second solver; wherein the mutual exchange of data between the first solver and the second solver comprises: - the first solver computes a first local prediction of common frontier data from current common frontier data, and sends the result of the computation to the second solver; - the second solver computes a second local prediction of common frontier data from current common frontier data, and sends the result of the computation to the first solver; - the first solver receives the second local prediction sent by the second solver; - the second solver receives the first local prediction sent by the first solver; - the first solver combines the first local prediction of the common frontier data with the received second local prediction of the common frontier data, the result of the combination being a global negotiated common frontier prediction data for the next time-step of the simulation; - the second solver combines the second local prediction of the common frontier data with the received first local prediction of the common frontier data, the result of the combination being the global negotiated common frontier prediction data for the next time-step of the simulation; - the first solver computes a first corrective force to apply on each of the paired nodes comprised in the common frontier to reach the global negotiated common frontier prediction data at the next time-step of the simulation, the first corrective force being computed from current local common frontier data and the global common frontier prediction data; and - the second solver computes a second corrective force to apply on each of the paired nodes comprised in the common frontier to reach the global negotiated common frontier prediction data at the next time-step of the simulation, the second corrective force being computed from current local common frontier data and the global common frontier prediction data. - providing a plurality of parallel solvers, a first solver simulating at least one first fabric piece and a second solver simulating at least one second fabric piece, the at least one first and second fabric pieces being modeled by surfaces meshes comprising nodes; and - maintaining the at least one first and second fabric pieces tight by computing a common frontier by the first and second solvers, wherein the common frontier comprises paired nodes of the at least one first and second fabric pieces, and wherein the computation of the common frontier is performed by combining data of the first fabric piece and data of the second fabric piece, the combining comprising reaching a common position target of each of said paired nodes; - the first solver computes a first local prediction of common frontier data from current common frontier data, and sends the result of the computation to the second solver; - the second solver computes a second local prediction of common frontier data from current common frontier data, and sends the result of the computation to the first solver; - the first solver receives the second local prediction sent by the second solver; - the second solver receives the first local prediction sent by the first solver; - the first solver combines the first local prediction of the common frontier data with the received second local prediction of the common frontier data, the result of the combination being a global negotiated common frontier prediction data for the next time-step of the simulation; - the second solver combines the second local prediction of the common frontier data with the received first local prediction of the common frontier data, the result of the combination being the global negotiated common frontier prediction data for the next time-step of the simulation; - the first solver computes a first corrective force to apply on each of the paired nodes comprised in the common frontier to reach the global negotiated common frontier prediction data at the next time-step of the simulation, the first corrective force being computed from current local common frontier data and the global common frontier prediction data; and - the second solver computes a second corrective force to apply on each of the paired nodes comprised in the common frontier to reach the global negotiated common frontier prediction data at the next time-step of the simulation, the second corrective force being computed from current local common frontier data and the global common frontier prediction data.