- bugfix

Author: janbender

examples/pbd.html

@@ -141,7 +141,7 @@ <h3>4. Position correction:</h3>
 				The position correction for a particle is determined using the Lagrange multiplier:
 					$$\begin{align*}
 					\Delta \mathbf x_{i_1} &= \frac{1}{m_{i_1}} \frac{\partial C_i}{\partial \mathbf x_{i_1}} \lambda_i = -\frac{\frac{1}{m_{i_1}}}{\frac{1}{m_{i_1}} + \frac{1}{m_{i_2}}}(\|\mathbf{x}_{i_1} - \mathbf{x}_{i_2}\|-d)\frac{\mathbf{x}_{i_1}- \mathbf{x}_{i_2}}{\|\mathbf{x}_{i_1} - \mathbf{x}_{i_2}\|} \\
-					\Delta \mathbf x_{i_2} &= \frac{1}{m_{i_2}} \frac{\partial C_i}{\partial \mathbf x_{i_2}} \lambda_i = +\frac{\frac{1}{m_{i_2}}}{\frac{1}{m_{i_1}} + \frac{1}{m_{i_2}}}(\|\mathbf{x}_{i_1} - \mathbf{x}_{i_2}\|-d)\frac{\mathbf{x}_{i_1}- \mathbf{x}_{i_2}}{\|\mathbf{x}_{i_1} - \mathbf{x}_{i_2}\|} 
+					\Delta \mathbf x_{i_2} &= \frac{1}{m_{i_2}} \frac{\partial C_i}{\partial \mathbf x_{i_2}} \lambda_i = +\frac{\frac{1}{m_{i_1}}}{\frac{1}{m_{i_1}} + \frac{1}{m_{i_2}}}(\|\mathbf{x}_{i_1} - \mathbf{x}_{i_2}\|-d)\frac{\mathbf{x}_{i_1}- \mathbf{x}_{i_2}}{\|\mathbf{x}_{i_1} - \mathbf{x}_{i_2}\|} 
 					\end{align*}$$
 				and the correction is applied to update $\mathbf x^*$ as
 					$$\mathbf x^* := \mathbf x^* + \Delta \mathbf x.$$
@@ -302,7 +302,7 @@ <h3>References</h3>
 						gradC_x = dx/dl; 
 						gradC_y = dy/dl; 
 						let K = (gradC_x*gradC_x + gradC_y*gradC_y) / this.mass;
-						if (i != 0)
+						if ((this.isActive(index1) && this.isActive(index2)))
 							K = 2*K;
 
 						// comute Lagrange multiplier

examples/xpbd.html

@@ -325,7 +325,7 @@ <h3>References</h3>
 						gradC_x = dx/dl; 
 						gradC_y = dy/dl; 
 						let K = (gradC_x*gradC_x + gradC_y*gradC_y) / this.mass;
-						if (i != 0)
+						if ((this.isActive(index1) && this.isActive(index2)))
 							K = 2*K;
 						K += alpha_tilde;